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20.  ABSTRACT  (Continue  on  reverae  aide  It  neceaaary  and  Identify  by  block  number) 

This  research  program  investigated  a  new  tool  for  the  analysis  of  elec- 
tromangetic  coupling  and  shielding  problems:  the  finite-difference, 
time-domain  (FD-TD)  solution  of  Maxwell's  equations.  The  objective  of 
the  program  was  to  evaluate  the  suitability  of  the  FD-TD  method  to 
determine  the  amount  of  electromagnetic  coupling  through  an  aperture  into 
an  enclosed  conducting  container  and  the  interaction  and  coupling  of  the 
penetrating  fields  with  internal  electronics.  Two  specific  container  \ 
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models  were  used  for  the  evaluation.  The  first,  a  conducting  cylinder 
with  one  open  end;  the  other,  the  guidance  section  of  a  missile.  Each 
of  these  two  configurations  was  modeled  to  calculate  the  electromagnetic 
field  ocupled  into  the  structure.  . 

The  following  specific  questi{hts"were  addressed  during  this  program 

1.  Can  the  FD-TD  method  accurately  model  electromagnetic  coupling 
into  a  conducting  structure  for  arbitrary  angles  of  incidence  and 
arbitrary  wave  polarization? 

2.  Can  the  FD-TD  method  accurately  model  a  complex  structure  with 
both  irregular-shaped  apertures  and  dielectric  or  permeable  materials 
within  the  interior? 

3.  Can  the  FD-TD  method  predict  the  voltage,  current,  or  power  that 
the  coupled  field  may  indice  on  wires  or  cable  bundles  inside  of  a  metal 
structure? 

4.  What  are  the  size  limits  of  an  object  that  can  be  modeled  using 
the  FD-TD  method? 

5.  What  is  the  capability  of  the  FD-TD  method  to  interface  with 
other  analysis  techniques,  such  as  "Method  of  Moments? 

Overall,  this  program  showed  that  the  FD-TD  method  can  be  success¬ 
fully  applied  to  electromagnetic  coupling  problems  involving  conducting 
structures  with  hole  and  sleeve-type  apertures  illuminated  by  a  plane 
wave  having  an  arbitrary  polarization  and  angle  of  incidence.  Further, 
the  FD-TD  method  can  be  applied  to  complex  cavity-like  structures  having 
internal  metal  and  dielectric  materials,  as  well  as  connecting  wires. 
Accuracy  of  the  FD-TD  results  was  very  good  relative  to  the  uncertainties 
of  available  experimental  and  numerical-theory  approaches.  Convergence 
of  the  leectromagnetlc  fields  to  the  sinusoidal  steady  state  occurred 
within  about  3  cycles  of  the  incident  wave  when  a  slight  value  of  iso¬ 
tropic  loss  was  assigned  to  the  interior  of  the  structures  modeled.  This 
resulted  in  program  central  processor  times  of  less  than  5  minutes  for 
FD-TD  lattices  containing  as  many  as  1.5  million  unknowns  time-stepped  to 
n  =  800  using  the  Control  Data  STAR-100  and  Cyber  203  computers. 

maX  This  program  also  established  the  feasibility  of  a  hybrid  MOM/FD-TD 
analysis  technique  based  upon  the  use  of  Schelkunoff ' s  equivalent 
electric  current  theorem  at  apertures  of  a  cavity.  This  hybrid  technique 
was  shown  to  give  consistent  results  for  apertures  strongly  coupled  to 
internal  wires  or  other  metal  and  dielectric  structures. 

This  program  has  shown  that  the  pure  FD-TD  and  hybrid  MOM/FD-TD 
methods  have  great  promise  for  applications  involving  complex  conducting 
and  dielectric  structures  illuminated  by  plane  waves  at  arbitrary  angles 
of  incidence  and  polarization.  Further,  great  promise  is  shown  for  those 
structures  that  are  simultaneously  electrically  large  compared  to  a  wave¬ 
length  and  penetrated  by  locally  complex  cavity-backed  apertures  having 
dimensions  comparable  to  a  wavelength.  Such  structures  requiring 
resolution  of  both  large  and  small  details  at  the  same  time  have  not  been 
well  treated  by  any  one  previous  analytical  or  numerical  approach. 
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PREFACE 


I IT  Research  Institute  (IITRI)  is  pleased  to  submit  this  Final  Report  on 
"Evaluation  of  Time  Domain  Electromagnetic  Coupling  Techniques"  to  Rome  Air 
Development  Center  (RADC/RBCT).  The  report  covers  work  performed  by  IITRI 
under  Air  Force  Contract  No.  F30602-79-C-0039,  designated  as  IITRI  Project 
No.  E6461.  The  report  is  in  two  volumes.  Volume  1  covers  details  of  the 
technical  work,  including  relevant  theory  and  numerical  results.  Volume  2 
provides  listings  of  the  Fortran  computer  programs  used  to  obtain  the  results 
of  Volume  1 . 

The  principal  investigator  on  this  program  was  Dr.  Allen  Taflove,  with 
support  provided  by  Dr.  Korada  Umashankar  of  IITRI  on  Task  5;  Mr.  Neil 
Robertson  of  IITRI  on  data  reduction  and  contour  map  plotting;  and  Prof. 
Donald  Wilton  and  Dr.  Allen  Glisson  of  the  University  of  Mississippi  on  com¬ 
putation  of  frequency-domain  data  for  checking  of  results  and  for  Task  5. 

The  project  duration  was  7  December  1978  to  29  February  1980. 

Respectfully  submitted, 

I  IT  RESEARCH  INSTITUTE 


Allen  Taflove,  Ph.D. 
Research  Engineer 
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VOLUME  2  -  LISTINGS  OF  FORTRAN  COMPUTER  PROGRAMS 


1.0  INTRODUCTION 

This  volume  documents  the  FD-TD  computer  programs  written  during  the 
present  research  effort.  Included  are  listings  of  the  programs  for  the  fol¬ 
lowing  problems: 

A.  Penetration  of  a  19.0  cm  diameter,  68.5  cm  long,  open-ended 
aluminum  cylinder  by  a  300  MHz  plane  wave  at  broadside 
incidence  and  TE  polarization  (Section  3.1  of  Volume  1); 

B.  Same  as  Problem  A,  but  for  TM  polarization  (Section  3.2  of  Volume  1) 

C.  Penetration  of  a  12.8  cm  diameter,  28  cm  long,  missile  guidance 
section  by  a  300  MHz  plane  wave  at  axial  incidence,  for  the  case 
of  the  interior  dielectric  components  modeled  (Section  4.1  of 
Volume  1 ) ; 

D.  Same  as  Problem  C,  but  for  the  case  of  the  interior  dielectric 
and  metal  components  and  wires  modeled  (Section  4.2  of  Volume  1); 

E.  Hybrid  method  of  moments/FD-TD  analysis  of  penetration  of  a  19.0 
cm  diameter,  68.5  cm  long,  open-ended  aluminum  cylinder  by  a 

300  MHz  plane  wave  at  axial  incidence  (Section  7.4.1  of  Volume  1); 

F.  Same  as  Problem  E,  but  for  TM  polarization  of  the  incident  wave 
at  45°  incidence  relative  to  the  cylinder  axis  (Section  7.4.2  of 
Volume  1 ) ; 

G.  Hybrid  method  of  moments/FD-TD  analysis  of  loaded  missile  guidance 
section  of  Problem  D  (Section  7.4.3  of  Volume  1). 


1 


The  computer  program  listings  for  Problems  A  and  B  were  run  without 
requiring  any  data  cards,  since  the  descriptions  of  the  cylinder  geometries 
employed  are  contained  within  the  programs.  Program  listings  for  Problems 
C,  D,  E,  F,  and  G  required  data-cards  decks  of  sizes  720  cards,  720  cards, 

72  cards,  72  cards,  and  528  cards,  respectively,  to  specify  the  interaction 
geometry.  The  data-card  format  for  these  programs  is  summarized  in  Section 
3.1  of  this  volume.  Further,  program  listings  for  Problems  E,  F,  and  G 
required  data-card  decks  of  sizes  144  cards,  288  cards,  and  225  cards,  re¬ 
spectively,  to  specify  the  aperture  excitation  for  the  hybrid  method  of 
moments/FD-TD  approach.  The  data-card  format  for  this  specification  is  sum¬ 
marized  in  Section  3.2  of  this  volume. 

All  FD-TD  computer  programs  were  written  using  STAR  Fortran  Version  2.1 
for  processing  by  the  Control  Data  STAR-100  or  Cyber  203  computer  systems 
under  the  1.2  operating  system.  This  Fortran  version  contains  certain  exten¬ 
sions  to  ANSI  standard  Fortran  [1]  that  permit  usage  of  the  vector  processing 
capabilities  of  these  computers.  The  reader  is  referred  to  the  STAR  Fortran 
Manual  for  detailed  discussion  of  these  features  [2]. 
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2.0  FD-TD  PROGRAM  LISTINGS 

In  each  of  the  programs  of  this  section,  FREQ  =  3.0  E  +  8  denotes  the 
operating  frequency,  f  =  300  MHz;  DX  denotes  the  lattice  cell  size,  6  (in 
meters);  MPR  denotes  the  total  number  of  media  within  the  model;  DATA  EPS 
and  DATA  SIG  give  the  assumed  relative  dielectric  constant  and  conductivity 
(mhos/m)  of  each  medium;  and  NMAX  gives  the  number  of  the  last  time  step  of 
the  algorithm. 

2.1  Problem  A  --  Task  1,  Case  1 
(Section  3.1  of  Volume  1) 

The  following  11  pages  list  the  computer  program  for  the  159  x  64  x  24 
cell  —  800  time  step  run  of  Problem  A.  The  problem  solved  is  penetration  of 
a  19.0  cm  diameter,  68.5  cm  long,  open-ended  aluminum  cylinder  by  a  300  MHz 
plane  wave  at  broadside  incidence  and  TE  polarization. 


r>  r»io  oio  r> 


C 

C 


C 


C 

C 


PROGRAM  F0T0( INPUT* OUTPUT #TAPE60« INPUT# TAPE8-TAPE8  ) 


RUN  TASK1 —  STEADY  300  MHZ  PLANE  WAVE  IRRADIATION  OF  A 
19.0  CM  DIAMETER*  68.5  CM  LONG*  OPEN-ENDED 

_ ALUMINUM  CYLINDER _ 

CASE  I-  TE  POLARIZATION  OF  THE  INCIDENT  WAVE 

_ BROADSIDE  INCIDENCE _ 

INCIDENT  WAVE  COMPONENTS  EZ  AND  HX 

159  X  64  X  24  CELL  CUBIC  SPACE  LATTICE  IS  USED _ 

UNIT  CELL  DIAMETER  *  DX  •  0.5  CM  «  WAVELENGTH/200 
EVEN  SYMMETRY  ABOUT  LATTICE  PLANE  Z  »  24PQX  IS  ASSUMED 
SOFT  TEM  WAVE  SOURCE  CONDITION  IS  USED  AT  PLANE  Y  •  3.0*DX 

SOFT  LATTICE  TRUNCATIONS  ARE  USED _ 

PROGRAM  IS  OPTIMIZED  FOR  THE  CDC  STAR-100 


REAL  MUR*MUZ 

DIMENSION  A  ( 2422500  )  * _ A1  (  3840  I  »  A2 1 3840 1 »  AAL (  38*0  >» 

1  A AR ( 38  4  0 ) »  DD( 160) #AL <24 )#AR<24) *CA<9), 

2  _ C8(  9  )  »  EPS  (  3)»SIG(3)» _ KE  XB  (  19  )  *KEXC  ( 19  ) 

COMMON  A 

DESCRIPTOR  D1 *  02*  03*  D4* NF D> BV> BW, B X* B Y> AIX* A2X* A1Y* A2Y* A1Z* 
1  A2Z 

BIT  BV*BBBC384Q1*BW*B8W(3840)»BX*BBX(3840)»BY»BBY(159) 
ASSIGN  A1X»A1(162;3678) 

ASSIGN  A2X*A2(162»3678) _ 

ASSIGN  AIY,A1(I61}3679) 

ASSIGN  A2Y*A2(161;3679) _ 

ASSIGN  A1 Z *  A1 < 1 . 3839  ) 

ASSIGN  A2Z*  A2 ( 1 »  3839  ) _ 

ASSIGN  BW*BBW(l;3840) 

ASSIGN  BX»BBX< 1 »  3840  ) _ 

ASSIGN  BY,BBY(l;159) 

ASSIGN  D4>QD(l;159l _ 

T1  -  SECOND(CP) 

PRINT  150*  T I _ 

150  FORMAT ( F20. 5 ) 


C  . I.  PROBLEM  PARAMETERS 

_ FREQ  »  3 .0  E*8 _ 

DX  »  0.005 

_ MPR  «  3 _ 

DATA  EPS/I. 0*  1.0*  1.0/ 

_ DATA  SIG/O.O*  3.7E*7»  0.01/ 

NM AX  ■  800 


C _ 

C  . CYLINDER  GEOMETRY  DESCRIPTION . 

_ DATA  KEX 3/19, 17*15»13»12*11»10*9»9»8»8*7»7»7»6» 6* 6*6*6/ 

DATA  KEXC/2 5* 19* 17* 15 *13* 12* 11* 10*9*9* 8*8 *7* 7* 7*6* 6*6*6/ 

C _ _  _ 

C  . IJ.  BASIC  AND  DERIVED  CONSTANTS . 


PI  »  3.141592&5 _ 

M'iZ  4.0  *  PI  *  1.0E-7 
;Z  «  8 .854E-12 _ 


A 


o  r>|o 


OT  *  OX  /  6.0E+8 

_ NHALF  -  0.5  /  FREQ  /  OT 

R  ■  DT  /  2.0  /  EPSZ 

_ R A  »  DT  *»2  /  PX»*2  /  MUZ  /  EPSZ 

RB  *  DT  /  OX  /  MUZ 

_ RC  »  1 .0  E«4  /  RB _ 

RD  •  2.0  *  PI  *  FREO  *  OT 
I  PUN  ■  0 
C 

Aid.  3840 )  *  0. 

A2  < 1 ; 3840 )  *  0. 

CAd;9)  *  0. 

CB(1J9)  *  0. 

B W  *  Q8VMKOd»160;BW  ) 

BX  *  Q8VMKZ (159*160; BX ) 

DO  2  1*1, MPR 

EAF  «  R  *  SIG(I)  I  EPS(I) 

CAd)  »  (1.0-EAF)  /  d.0«-EAF) 

2  CB(I)  *  RA  /  EPS(I)  /  (l.OfEAF) 
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C 


A 


C 

c 


c 
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. III.  LOAD  VECTOR  A. 

. ZERO  INITIAL  FIELDS 

DO  3  1*1,65 

I PEL  «  (1-1)  ♦  36900 _ 

A ( IDEL+l ; 36900)  *  0. 
A(2398501;24000)  *  0. 


. TYPE  OF  MEOIUM . 

. ANISOTROPIC  LOSSY  AIR 

DO  A  1*2,64 _ 

IDEL  *  JI-l)  *  36900 

A C IDEL* 1 ; 3840 )  *  3. _ _ 

AI IDEL*8161;3840)  *  3. 

A( IDEL+16321; 38401  «  1. _ 


. . broadside  circular  cylinder 

DO  14  1*15,33 

IDEL  *  (I-l)  *  36900 _ _ _ _ _ 

KB  *  KEXB ( 1—14 ) 

KC  «  KEXCU-14) _ 

00  7  K*KB, KC 

A ( IDE L»1 2 ♦(K-ll»160; 1371  *  2. _ 


A(I0EL»8171»(KB-l)»160;i3g)  *  2. 
DO  10  K  »KB, 24 

A( IPEL*8308*(K— 1)*160)  *  2. 

KC  »  KE  XC ( 1-14  >  ~  1 _ 

DO  12  K*KB,KC 

AdDEL*16331*<K—l)*160;i 38)  *  2. 
KD  *  KC  ♦  1 

DO  13  K«KD, 24 _ 

A C  IDE L*1 6331* ( K—l ) *160;  137)  •  3. 


! 

i 


13  A  (  IDEl  *16468* (  K-l  >♦  1  6 0  1  »  2. _ 

C 

_ IA  »  67  -  I _ 

IDELA  *  CIA-1)  *  369C0 

_ A (  IQELA*l53840)  «  A ( IOEL* 1 ; 3840 ) _ 

AC  IDELA*1632l53840)  *  A( IDEL*16321 5 3840) 

_ IB  »  66  -  I _ _ _ 

IOELB  -  C  IB-1  )  *  3690C 

14  AC  I 0ELB*8161 5  3  84 0)  »  AC  IDEL*8161| 3840) 

C 

_ T 2  »  SECONQ(CP) _ 

PRINT  150, T2 

C _ _ _ 

c  . IV.  TIME-STEPPING  LOOP . 

00  200  N «1 # NMA X 

TERM  «  SIN (FLOAT (N)*RD)  *  RB 

MC  ALL  »  3  ♦  IFIXCFL0ATCN)/6.0) 

IF (MCALL.GT.21 )MCALL*21 


r 

C  7. ...TRANSVERSE  PLANE  NO.  1 

C  . EX,  EZ  TRUNCATIONS.... 


A  C  400 1 ; 4  000 )  *  ACi;4C00) 

_ A(l;4000)  »  A (40901 5  4000 ) _ 

A ( 20321 5 4000 )  *  A( 16321 54000) 

_ AJC  163  21  >4000 )  «  A  (57  221 ,4000 ) _ 

C 

C . . EY  ITERATION . . 

ASSIGN  01#. OYN. 3679 

_ ASSIGN  02, .OYN. 3679 _ 

ASSIGN  03, • OYN. 3679 

_ A  ( 1 21 61 )  »  0.5  »  (AC  16161 )»AC 16162) ) _ 

A ( 12162; 157 )  -  0.333  *  ( A C 161 61 5 15 7 > *  A ( 1 6162 ; 157) 

1 _ *A (16163 , 157) ) 

A ( 12319)  *  0.5  *  (AC  16318 )*A< 16319)  ) 

_ A  C 161 61 5 1 59  )  -  A (  123  21? 159 ) _ 

01  *  C A ( 3  )  *  A(12321  53679  ) 

_ 02  *  A(24641; 3679)  -  A(  24 48153679)  ♦  A ( 33011 5 3679 ) 

1  -  A ( 3301 2  5  36  79 ) 

_ 03  «  CBC  3)  ♦  02 _ _ _ 

A  C 123 21 5  3679 )  «  D1  *  03 

_ FREE _ 

C 

C . . TRANSVERSE  PLANES  2-64 . . 

00  82  J  Y* 1 , MC  ALL 

_ JOEL  -  (JY-1)  *  110700  _ 

C 

C . . EX  ITERATION . . 

ASSIGN  Oi, .OYN. 3678 

_ __  ASSJGN  D2,.0YN.  3678_ _ 

ASSIGN  03, .OYN. 3678 

A  S  S 1 5  N  NJFD  ,  .  OYN. 3678 _ 

ASSIGN  BV,BBB( 153678 ) 


A(40902*H)  «  0.5  *  (  A(44902*HI  ^(44903^)) 


Hi  u'ltHUU  ClvI'r'U.fl  t-I-i  Ui  CIlI'FUifl  H-1 1 


1  «A(44904*M;156) ) 

A(41059+H)  »  0.5  *  ( A(45058»H)»A(45059»H) 


A(44902*M»158)  «  AC 41062*M Jl58) 


. MAIN  EX  LOOPS . 

NFD  »  At 37062»Hi3678 ) 


DO  22  JJ*1#MPR 
•  NFD. 


A1X  *  08VCTRLCCAC JJ >,BV;A1X> 
A2X  «  Q8VCTRL(C 8( J4 ) t BV : A2 X ) 


01  •  A1X  *  AC  41062*M ; 36  78 ) 

«  A(69912+H;3678)  -  AC33012* H;3678)  ♦  A(65542» M .3678 ) 


1  -  AC65702+M ;3678 ) 

03  ■  A2X  ♦  02 


30  AC41062+MJ3678)  ■  01  ♦  03 
FREE 


. EY  ITERATION. 


# . OYN. 
ASSIGN  02# . DYN. 3679 
ASSIGN  03#. OYN. 3679 
ASSIGN  NFD#. DYN. 3679 


ASSIGN  B V# BBB ( 1 # 3679 ) 


00  40  MA-1»3 
M  -  JOEL  ♦  36900* (HA— 1 ) 


. SOFT  LATTICE  TRUNCATION . 


A ( 49061 ♦H )  ■  0.5  *  ( A(53061*M)*AC 53062+M) ) 

AC 49062+H# 1 57 )  »  0.333  *  ( A( 53061 ♦« ;1 57 ) *A ( 53062»M ; 1 57 ) 


♦A(53063*N;157> I 

A ( 49219+M )  »  0.5  ♦  ( A ( 53218*H) *A C 53219*H) > 


A (53061+Mj 159 )  -  A(49221*H;159) 


. MAIN  EY  LOOPS . 

NFO  *  AC  45221 *M  #3679 ) 


DO  32  4J-1#MPR 
BV  *  NFQ.EQ.JJ 


»  Q8VCTRLCCAC JJ)»BV?A1Y) 
32  A2Y  •  Q8VCTRLCCBC JJ)»BVJA2Y> 


01  «  A1Y  *  A(49221*M?3679) 

02  »  A(61541*M»36  79)  -  AC  61 381*M; 3679  I 


03  «  I2Y  *  02 


40  A(49221*HJ3679)  ■  01  ♦  03 
FREE 


♦  A(69911*MJ3679) 


-  A ( 6991 2* M  J  3679 ) 
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i 


c 


c 


c 

c 


£ 


c 

£ 


c 

£. 


c 


c 


. . SOFT  LATTICE  TRUNCATIONS . . 

. LEFT  TRUNCATION . 

AL<1)  »  0.5  ♦  <A<69541+H)+A<69542+H)) _ 

AL 12; 22 )  •  0.333  *  < A<6954i+M ; 22 )♦ A (69542+H ; 22 ) 

1 _ »A<69543  +  H;22) ) _ 

AL<24)  *  0.333  *  < A ( 69563+N ) ♦ 2. 0  +  A ( 69564+ M ) ) 

AAL ( 1 ; 38 AO )  »  Q8VXPN0<AL(1;24),BW:AAL<H3840)) _ 

A ( 69541+H  »  24 )  -  Q8VCHPRS<A<65542+M;3840)*BW>A<69541+H>24)) 


. RIGHT  TRUNCATION..... 

ARID  -  0.5  »  <A<69566+H) +A<69567»H)) _ 

AR (2;22)  «  0.333  *  < A I69566+H J 22 )♦ A I69567+H; 22 ) 

1 _ »A  < 69568  »rt:  22  ) ) _ 

AR <24 )  *  0.333  *  < A < 69588+M ) +2.0+A < 69589+H > ) 

AARC1 ;3840 )  •  08 VXPND ( AR ( 1 ; 24 ) , B X ; A AR C 1 ; 3840 )  ) _ 

A ( 69566+H »  24 )  -  Q8VCHPRS< AI65540+H; 3840)* BX| A<69566+H;24 ) ) 


. MAIN  HY  LOOPS . 

D1  »  A<65542+H;3838) _ 

02  «  A<57222+M*3838)  -  A< 57221+H; 3838 )  ♦  A (40902+M ; 3838 ) 

1 _ ~  A<41062+H;3838) 

A<65542+H;3838)  -01+02 


A<65541+M;3840)  «  00 VCTRL < AAL < 1 ; 3840) * BW? AI65541+M 13840 ) ) 
_ A  <  fc.554.ltf1  j  3840 )  -  08  VCT  RL  <  AAR  <  1 ;  384Q)  *  BX;  A  ( 65541+H  ;  3840  )  ) 

. . HY  HOR.  PLANE  ENVELOPE  COMPUTATION . 

04  -  VABS<A<69221+M*159);04) 

_ BY  -  P4«GT«At6959.1+j;lj.91 _ 

70  A ( 69591 +H ; 1 59 )  »  Q8VCTRLC D4, B Y; A< 69591 +MJ 1 59 ) ) 

_ Oil _ 


. HZ  ITERATION 

ASSIGN  01* .DYN .3678 
ASSIGN  02* .0YN.3678 


00  80  HA-1*  3 
N  -  J DEL  ♦  36900* <MA-1 ) 


. SOFT  LATTICE  TRUNCATIONS . 

. . LEFT  TRUNCATION . . 

AL  < 2)  -  0.5  ♦  <A~<36852+H)+A<36853+M)) 

AL  <  3  >  21 )  -  0.333  ♦  < A( 3 68 52+H ; 21 ) »A < 36853 +H; 21 ) _ 

I  +  A ( 36854  +  N ; 2 1 )  ) 

AL  <24 )  -  0.5  ♦  <A<36873+H)+A<36874+H)  ) _ 

AAL ( 1 >3840 )  -  08VXPNO<AL<1>24)*BW*AAL<1;3840)) 

A< 36851+H*24)  -  Q8VCHPRS< AI32852+H; 3840) » BW; A< 3685 1+H ; 24  ) ) 


. . RIGHT  TRUNCATION . . 

AR<2)  -  0.5  *  <A<36877+H) +A(36878+H)) 

AR  <3»  21 )  »  0.333  ♦  < A < 368 77+H ; 21 )♦ A ( 36878+H i 21 ) 

L  +A<  36879 +H»  21 ) ) 
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_ ARC24)  «  0.5  *  CAC36898»H)»AC36899»H)) _ 

AAR  C 1 ; 3840 )  •  Q8VXPN0 ( AR< 1 ; 24 ) , BX ; A AR ( 1 ; 3 040 I) 

_ AC36876»H;24)  -  Q8VCHPR S ( A ( 32 850+H) 38401 , BX ) AC36876*H)24 >) 

C 

C  . MAIN  HZ  LOOPS..... 

01  •  AC33012+M;3678) 

_ 02  «  A(41062*H;3678)  -  AC 4162+H; 3678)  ♦  AC12321+H; 3678) 

1  -  AC 12322+M; 3678) 

_ AC33012*H;3678)  »  01  ♦  02 _ 

C 

_ AC328511H13840I  -  QB VCT RL ( AAL. fl ) 3840 ) , BWt A( 32851*H )3840 1) 

80  AC32851+MJ3840)  *  Q8VCTRI C  AAR  C 1 J 3840) * BXJ  AC  3285 i*Hj 3840 )  ) 

_ FREE _ 

82  CONTINUE 


_ IFCHCALL.LE. 20160  TO  94 

C 

C  . TRANSVERSE  PLANE  NO.  65 

M  «  63  ♦  36900 

C _ 

C  .....EX,  EZ  TRUNCATIONS 


_ A  C  40901 »M ;  4000 )  »  A( 36901 ♦H;4000) _ 

A(  36901 *M» 4000 )  -  A C  45061+M; 4000 ) 

_ AC4S061*H;4000)  •  A(  4001*  Mi  4000) _ 

C 

_ ACS7221*H;4000)  »  AC 53221»H;4000) _ 

AC53221*M;4000)  -  A  ( 61381  *M;  4000  > 

_ AC6l381»N;4000)  «  A ( 20321 *H ;4000 > _ 

C 

C . . HX  ITERATION . . 

ASSIGN  D1,.0YN.3639 

_ ASSIGN  D2,  . OYN» 3839 _ . 

01  «  AC2448i*MJ3839) 

_ 02  «  AC12321*M?3839)  -  A( 1 2161» W ; 3 839 )  ♦  AC 20321+H? 3839 > 

1  -  AC  57221*N ; 3839 ) 

_ AC24481»H;3839)  »  01  ♦  D2 _ 

D4  «  VABSC AC28161+M;159);D4> 

_ 8 Y  »  D4.GT.AC28481»H;159) _ 

A  C  284  81 *M  » 1 59 )  »  Q8VCTRLC D4,BY; AC28481+H) 159) ) 

_ FREE _ 

C 

C . . HZ  ITERATION . . 

ASSIGN  01,. DYN. 3678 

_ ASSIGN  02, .DYN. 3678 _ 

At (2 1  «  0.5  *  CAC36852*MI*AC36853*H)) 

AL(3;?1)  «  0.333  ♦  C AC36852»H;21)*AC36853*H;21) _ 

I  ~  ♦AC36854»M;21l I 

_ ALC24)  «  0.5  *  CAC36873*H)*AC  36874  »H)  I _ 

AALC1  ,3840)  -  Q8VXPN0C ALC 1 ; 24 ),BV; AAL C 1; 38401) 

A  C  36  8  51 *M  » 24)  »  Q8VCHPRSCAC32852*H) 3840) , 8W I  AC  3685 !♦ Ml  24 ) ) 
AF  C  2  }  *  0.5  ♦  CAC  368  77*H)*AC3  68  78*11)) 

Ai  <3)21)  »  0.333  ♦  C AC36877*H)21 )»A C36878»H;21 ) _ 

1~  *AC36879*H;2i)) 


10 


AR  <  24  >  ■  0.5  »  (  Af36898»M)»At36899»M)  ) _ 

A ARf 1 ; 3840 )  »  Q8VXPN0f ARfl;24),BX;AARfl;3840)> 

A C  36876+M » 24 )  «  QBVCHPRSf A( 32850+M; 3840) , BX t At 36876+H ;24 ) I 
01  *  Af33012+N*3678) 

02  »  Af41062*M?3678)  -  AC 4162+H { 3678)  *  A I12321+H; 3678) 

1  -  A  f  1 2  322  »  36781 

Af 33012+H;3678)  »  01  ♦  02 _ 

At  328  51*M  J  3840  )  -  Q8VCTRL  UAL  <1  f  3840)  ,  BWJ  A  (  32851*H  f  3840 1  ) 
Af32851*M«3840)  -  08 VCT RL f A AR f li 3840) , BX i A f 32851  +  M ; 3840 )  ) 
FREE 


. VERT.  PLANE  ENVELOPE  COMPUTATION 

94  ASSIGN  01,  .DYN.3840 _ _ 

ASSIGN  02* . DYN. 3840 

ASSIGN  B V*  BB8 f  1 » 3840  ) _ 

M  »  32  *  36900 

HZ  «  65  ♦  36900 _ 


. EX  ENVELOPE . . 

01  •  0.5  ♦  (A(4001+M;3840)*A(40901*N;3840)) 

D2  -  VABSf  01; 02 ) _ 

BV  -  02.GT.Af MZ*1J3840> 

Af  HZ  +  1 ;  3840 )  »  Q8VCTRLf02,BV:AfHZ»l;3840) ) _ 

. EY  ENVELOPE . . 

02  -  VABS(A(12161+M}3840) ;D2) 

BV  »  02.GT. Af  MZ+4001 1 3840  I _ 

AfHZ* 4001*3840)  *  08 VCTRL C 02. BV; Af HZ*4001 J 3840) I 


. EZ  ENVELOPE . 

D1  »  0.5  ♦  fAf20321+H;3840)*Af57221*H;384Q)) _ 

D2  -  VAB S  f  01 ; D2  ) 

BV  ■  02.GT.AfMZ*8001;3B40) _ 

AfMZ+8001t3840)  »  08 VCTRL  f  02*  BV; Af  MZ^BOOl > 3840 ) > 


. HX  ENVELOPE....  • 

D2  »  VABSf A(24481*M}3840) :D2) _ 

BV  -  D2.GT • Af MZ+12001 ; 3840) 

AfHZ+12001;3840)  »  Q 8 VCTRL  fD2,BV)Af HZ ♦12001*  3840) ) 

. HY  ENVELOPE . . 

01  ■  0.5  ♦  fAf 2864l4M;3840)+Af65541tH?3840) ) 

D2  -  VA8SfDl;02) _ 

BV  *  D2.GT.AfHZ*16001j3840> 

Af HZ*16001:3B40)  «  OB VCTRL f 02 . B V ; A f MZ+16001 t 3840 )) 

. HZ  ENVELOPE . . 

02  «  VABStAf32851+M;3640) » 02 ) 

BV  »  02. GT.AfHZ*2 0001*3840) _ 

AfHZ* 20001*3840)  «  Q 8 VCTRL  fD2,BV»  Af  MZ* 20001 *3840) ) 
FREE 

T3  »  SECOND f CP ) 

PRINT  150,  T3 


c 

C . . FIELD  ENVELOPE  PRINTOUT  ROUTINE 

DO  100  L  *NH ALF  , NM AX, NHALF 
I F ( N. EO.L ) GO  TO  101 
100  CONTINUE 

IF (N«  EQ.NNAXIGO  TO  101 
GO  TO  199 


C  . AT  HORIZONTAL  SYMMETRY  PLANE . 

101  IFIN.EQ.NMAX)IPUN»1 
PRINT  102*  N 

102  F0RMATC1H1,52X*27HEZ  ENVELOPE  FOR  TIME  STEP  -»I5* 

1  /Z»62X»15HPLANE  Z  »  24*DX,  / /,2X,  l  HJ,  /  ) 

CALL  ENVH!24321,RC»IPUN)  _ _ 


PRINT  103,  N _ 

103  F0RMAT!1H1,52X,27HHX  ENVELOPE  FOR  TIME  STEP  ■•IS, 

1 _ //»62X,15HPLANE  Z  ■  24»DX,  / /,  2  X,  1H  J»  / ) 

CALL  ENVH!28481,3.77E*6,IPUN» 


PRINT  104,  N 

104  F0RMAT!1H1,52X,27HHY  ENVELOPE  FOR  TIME  STEP  »>I5, 

1  //,62X,15HPL ANE  Z  *  24*DX» //* 2X* 1HJ, / ) 

_ CALL  ENVH!32691,3.77E*6,IPUN) _ 

C 

C . . AT  VERTICAL  SYMMETRY  PLANE . . 

PRINT  105,  N 

105  FORMAT! 1H1, 52 X,27HEX  ENVELOPE  FOR  TIME  STEP  »,IS, 

1  //,62X,17HPLANE  Y  »  32.5*DX , / t ,2X, 1HK, / ) 

CALL  ENVVU,RC,IPUN> _ 


PRINT  106,  N _ _ _ 

106  FORMAT! 1H1»5ZX»27HEY  ENVELOPE  FOR  TIME  STEP  »,I5, 

1 _ //,62X,17HPLANE  Y  ■  32. 5*DX  ,  /  /  ,2X,  1HK,  /  ) 

CALL  ENVV(4001,RC#IPUN) 


PRINT  107,  N 

107  FORMAT! 1H1, 52 X»27HEZ  ENVELOPE  FOR  TIME  STEP  »,I5, 

1  // »62X» 17HPL  ANE  Y  *  32.5ADX,/ Z,2X,1HK,/ I 

CALL  ENVV!8001,RC,IPUN) _ 

PRINT  108,  N _ 

108  FORMAT! 1H1, 52X, 27HHX  ENVELOPE  FOR  TIME  STEP  »,I5» 

1 _ //,62X,17HPLANE  Y  ■  32.5»0X,//,2X,1HK,M 

CALL  E !VV(12001,3.77E+6,IPUN> 


PRINT  109,  N 

109  FORMAT! 1H1, 52X, 27HHY  ENVELOPE  FOR  TIME  STEP  »,I3, 
I  // ,62X,17HPLANE  Y  «  32.5*0X,//,2X,1HK, / ) 

_  CALL  ENVV!16001,3»77E*6,IPUN) _ 

_ P F INT  110,  N _ 
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110  F0RMAT(1H1»52X»27HHZ  ENVELOPE  FOR  TINE  STEP  »,I5* 

1 _ //,62X,17HPL  ANE  Y  »  32. 5*0X  ,  t  /  ,2X,  1HK,  /  ) 

CALL  ENVVI20001*  3  .77E*6»IPUN) 


199  CONTINUE 

200  CONTINUE 

T A  ■  SECOND(CP) 
PRINT  1 50,  T4 
STCP 

ENQ _ _ 


SUBROUTINE  ENVH!LOCA, SCALE, IPUNCH) 
DIMENSION  A (2422500) , IP  110240 ),NN( 160> 
COMMON  A 

z - 

00  9  1*1,160 
9  NN ( I )  •  r 
C 

DO  1  LY*2*64 

LOC  *  LOCA  ♦  (LY-11436900 
LdCI  -  1  ♦  (L Y-l ) *160 
IP(LOCI;160)  •  SCALE  *  A(L0C;160) 

1  A(L0CJ160)  *  0. 

C 

IF(IPUNCH.EQ.0)6Q  TO  4 

_ DO  2  LY-15,52 _ 

LOCI  «  1  ♦  (LY—1 ) *160 

_ LOCH  »  LOCI  »  158 _ 

2  WRITE (8»3)  ! IP(LL),LL*LOCI, LOCII) 

3  FORMAT! 1615 ) _ 

C 

4  DO  7  LX«1,6 _ 

L X A  *  5  ♦  !LX-1)*25 

_ LXB  «  LXA  »  24 _ 

DO  5  LY«2,64 

_ LYY  ■  66  -  LY _ 

LOCI  «  LXA  ♦  (LYY-114160 
_ LOCII  »  LOCI  »  24 _ 

5  PRINT  6,  LYY,  ( IP!LL )»LL-l OCI, LOCII) 

6  F0RMAT(1X,I2,2X,25I5) _ 

7  PRINT  8,  (NNCLL  )  > LL*LX A ,L XB ) 

8  FORMAT!//, 4X, 2515,////) _ 

C 

_ RETURN _ 

END 


SUBROUTINE  ENVV! L OCZ , SC  ALE, IPUNCH ) 

DIMENSION  A(2422500) ,NN (160 ),IPR(25),IPP(160) 
COMMON  A 

z - — - — - 


DO  8  1*1,160 


LHN  »  LQCZ  »  65*36900 


IPf IPUNCH.EQ.O)GQ  TO  3 _ 

00  1  LZ«l,24 

LOC  »  LHN  ♦  ( L Z-l  ) *1 60 _ 

IPP<1;159)  *  SCALE  ♦  ACL0CJ159) 

1  WRITE  <8,2 )  (IPP(LL),LL«1,159) 

2  FORMAT! 1615) 

3  DO  6  LX»1 »6 

L X A  «  5  ♦  ( LX-1 ) *25 _ 

LXB  *  LX A  ♦  24 

DO  4  LZ» 1, 24 _ 

LZZ  -  25  -  LZ 

LOC  *  LHN  ♦  ( LZZ-1)  *160  »  LXA-1 
IPRC1525)  «  SCALE  *  A!L0C;25> 

4  PRINT  5,  LZZ*  < IPR!LL ),LL  »1,25) 

5  FORMAT! IX,I2»2X,25I5) 

6  PRINT  7,  <NN!LL),LL«LXA,LXB) 

7  FORMAT!//, 4X, 2515,////) 

A ( L MN ; 3840 )  *  0. 

RETURN _ 


2.2  Problem  B  --  Task  1,  Case  2 
(Section  3.2  of  Volume  1) 

The  following  11  pages  list  the  computer  program  for  the  159  x  64  x  24 
cell  --  800  time  step  run  of  Problem  B.  The  problem  solved  is  penetration 
of  a  19.0  cm  diameter,  68.5  cm  long,  open-ended  aluminum  cylinder  by  a  300  MHz 
plane  wave  at  broadside  incidence  and  TM  polarization. 
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I 


o  |c->  o  |o 


VEL3Clli4'J 


PUN  TASK  1 --  STEADY  ->nO  MHZ  PLANE  WAVE  IRRADIATION  OF  A 
9.0  CM  niAMETFR.  (-8.5  CM  LONG.  OP 


aluminum  cylinder 

ASE  TI-  TN  POLARIZATION  OF  jHt  INCIDENT  WAV 


PRO  ADS  1 OE  IN'C  IOF.MCE 
N 


ISP  X  <S4  X  24  CELL  CUBIC  SPACE  LATTICE  IS  USED 
NIT  CELL  DIAMETER  =  DX  =  0. 


EVEN  SYMMETRY  ABOUT  LATTICE  PLANE  ?  =  24.S*0X  IS  ASSUMED 


Bwcmnn 


soft  lattice  TRUNCATIONS  are  used 

PROGRAM  IS  OPTIMIZED  FOR  THE  CD 


REAl  MUR .M 


DIMENSION  A (242P500) , CEXP ( 24S760 ) » CEYB ( 245760 ) , CEZB ( 24B760 ) 


mm 


»  CP (9) .FPS(S) .SI6<3> ♦ 

XH. 


K F. X 0  ( 1 9 )  *KEXC  (19) 


CA(1I9)  *  0. 


00  1  KK=1*?4 


DO  2  I*1*MPR 


CA ( I )  s  (1.0-EAF)  /  (1.0*EAF) 
EPS (I)  /  M  .0*£AF 


.ZERO  INITIAL  FIELDS, 


Mill 


IOEL  =  (1-1)  *  36900 


A <2398501 1P4000)  =  0. 


.TYPE  OF  MFHIUM. 


■rjuuimiminiM-i 


tt-Vii 


A  ( IOEL* 1 I 3840 )  *  1. 


4  A(IDEL*163?1I3840)  =  3 


.BROADSIDE  CIRCULAR  CYLINDER, 


nin 


IDEL  «  <I-1)  *  36900 


klaVill'i 


KB  a  KEXBU-14) 


C . . WM1IH  me  r  r  i  Ut.jua..  .  . . 

ASSIGN  A1..DYN.3840 

_ ab<;k.n  A2..nvN.3BA0 - 

ASSIGN  NFD, .DYN, 3840 

_ ft.  I  - 

A?  =  0. 

_ DO  23  I  =-L»6-4 - 

IDEL  *  (1-1)  *  36900 

_ IDEM  -s—t-U-U  4-364-0 - 

C 

_ NEfl  a  .Al.LDFL*  1 1.384,0.) - 

DO  20  JJS 1 » MPR 

_ BV  a  NF-CUEQ.J.J - 

A 1  =  OBVCTRL (CA (JJ)  .BVIA1 ) 

20  A2  =  OBVCTRL(Cb(JJ)  .RV1A2) _ 

A ( I DEL* 1  *  3840 )  =  A1 

_ CEXB( IOEM+1 1364Q)  =  a2 _ 

C 

_ NFD  =  A(I0EL*6161I384Q) _ 

DO  21  JJ=1,MPR 

_ BV  =  MFD.EQ.JJ _ 

A  1  s  08VCTRL <CA <JJ> ,PVIA1 > 

21  A?  =  Q8VCTRL (CP ( JJ) . PVI A2) _ 

A<IDEL»8161 13840)  *  *1 

_ CEYB(IDEM*1I3B40)  s  a2 _ 

C 

_ NFD  =  A(IDEL*l632H3P40) _ 

DO  ?2  JJ=1«MPR 

_ PV  s  NFD.FQ.JJ _ 

A 1  =  08VCTPL (CA (JJ) ,PV|A1 ) 

??  A 2  =  OPVCTRL(CB(JJ) ,RV1A2) _ 

A(!0EL+163?lt 3840 )  =  A1 

23  CEZBnOEWll384Q)  s  A 2 _ 

FREE 

_ A  ( 8 1 6 1  I  3840 )  =  0. _ 

C 

_ T2  =  SECOND (CP) _ 

PRINT  150, T2 

C _ 

C  . IV.  TIME-STEPPING  LOOP 

_ DO  200  Nal.NMAX _ 

TERM  s  SIN(FLOAT(N)*PD) 

MCALL  =  3  ♦  IFIX(FLOAT(N)/6.0) 

I F  ( MCA’  L.GT.21)MCALl=21 


C  . TPANSVERSF  PLANE  NO.  1 

C _ ....  HZ  truncations.... 

A  (4001*4000)  =  A  (816]  14000) 

_ A(8161  14000)  -  A(40R01>4000) 

A (20J21 *4000)  =  A (28641 *4000) 
_ A  (28641  *4000  ~  A  (57221  *4000  ) 


ASSIGN  D1..0YN.3679 


A  ( 12161 )  =  0.5  *  (AM6161)*A(1M6?)  ) 


l  ♦AM61&3I157)) 


A  ( 1 61 61  *  1  69)  =  A  ( 1 2??1  ?  159) 
M  73? 1  1 3679 I 


D?  =  A {24641  13679)  -  A { ?44fc 1 1 3679 )  ♦  A < 330 1 1 1 3679 ) 


(MilMknCJ 


A ( 12321 1 3679)  =  01  -  D2 


DO  62  JY=1«MCALL 

-li  e  ii 


JDEM  s  (JY-1 )  *  11520 


.....HX  ITERATION, 
N. 


ASSIGN  D2..DYN.3676 


DO  30  MA= 1 . 3 


*  (Ma-1 l 


A(4090?*M)  =  0.5  *  ( A ( 44902*M ) ♦A(44903+M) ) 


. MAIN  HX  LOOPS . 

D 1  a _ A  (4106?*m«367 


D2  s  A (69912  +  MJ3678)  -  A ( 330 1 2*M I  3678 )  ♦  A ( 6554?*M » 3678 ) 

S 


30  A  < 41 06?*M  j 3676 )  s  Dl  -  D? 


ASSIGN  D1..0YN.3679 
N 


r:Knin«*Mfiu 


M  s  JDEL  ♦  36900* ( MA- 1 ) 


.SOFT  LATTIOE  TRUNCATION, 


A ( 4906?*M ; 157)  =  0.373  *  ( A < 5306 1 *M 11 57 ) ♦ A < S3062*M 1 1 57 ) 


A (49219*M)  =  0.5  *  ( A <5321B*M) *a (53P19*M) ) 
is  iUQ??i.Mtmwi 
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I  |  |ooi  |  oi  mo 


DO  60  MAs] »  3 


M  a  JDFL  ♦  36900* (M a-1 ) 
MM  r  JOEM  ♦  36A0*<M/s-l) 


OOPS . 


01  =  A  ( ?4481 *M  t 3839 )  *  A(1*M»3839) 

0?  =  A  ( 1?3?1»M13B39)  -  A  (  1?161*M<3P39)  ♦  A  <  ?03P1  I  3P39  > 


1  -  A<S7?21*MI3839) 

03  *  CFXR ( 1*MM<  3839)  *  D? 


A(?4«P1*M13239)  =  D)  -  D3 


00  70  MAs 


M  =  JOEL  ♦  36900* (MA-1) 
MM  a  JDEM  ♦  3840* (MA 


AL  ( 24 )  =  0.333  *  (A (A9563*M) *2.0*A (69564*M) ) 


. RIGHT  TRUNCATION . 

«  (A(69^*6*MW 


AR  ( ?  I ?2 )  =  0.333  *  (A  (69566+MJ 2?)  +A (69567*Ml?2) 

( 


AR ( 24 )  =  0.333  *  (A (*9588*M) ♦?.0*A (695894M) ) 


DO  61  KK*1 « 24 
Mf 


MR!  a  M  ♦  MQ(KK) 

>  a  A ( 65S42  *MP  I ) 


61  A (69565*MLE )  *  A(65699+MRI) 


. MAIN  EY  LOOPS . 

01  =  A<6554?*M13A3A]  *  A { 45062*M » 383B 


DP  *  A (57222+MJ3838)  -  A (57221 ♦M13838)  ♦  A (4090?*Ml 3838) 

-  A ( 4 1 06P*M 


D3  =  CE YB  ( 3 84 2 ♦MM  <  3 A3fl }  *  02 
♦MI3836)  a 


MLE  a  M  ♦  MG(KK) 

541 ♦MLt) _ *  AL(Kk 


69  A(65700*MLE)  =  AR(KK) 


. EY  HOR.  PLANE  ENVELOPE  COMPUTATION.. 

BS  (A  (69? 21  ♦MM  59 )  »D4) 


=  D4.GT.A(69591*m»159) 

6959 1 +M  t 1 59 )  a  08V0TR1  ( DA . BY t A < *989 1 *M t 1 69 >  1 


FREE 


iEZ  ITERATION, 


r~wi  in  ?m  ii  vn  uimMi 


ASSIGN  02..0YN.3678 

N. 


M  =  JOEL  ♦  36900* (MA-1 ) 


C  . LEFT  TRUNCATION . 

_ AL  ( 2 )  =  0.5  »  (A  (3t>P^2*M)  *A  (36653+M)  ) _ 

AL(T;?1  )  =  0.333  *  <A<36852*M»21)*A<36R53*MI?1) 

_ 1 _ ♦  A  ( 36a54*H  )  2  H  ) 

AL(?4)  =  0.5  *  (A <3A«73*M) *A (3h874*M) ) 


C  . RIGHT  TRUNCATION . 

_ AR(2)  =  0.5  »  (A (36877+H) *A (36878+M)) _ 

AR (3121 )  =  0.333  *  (A(36877»MI?1)*A(36878*MI?1) 

_ ] _ ±A  (3-6BLS*HJ.2.H ) 

AR ( 24 )  =  0.5  *  (A (36A98+M) +A (36899+M) ) 


00  71  KK=1.?4 

_ MLf—S  H  ♦  KK _ 

MR  I  =  M  ♦  HQ(KK) 

_ A(3685n*HlE)  =  AL32^Z»al) 

71  A<36875*MLF)  =  A(33009*MRI) 


C _ 

C  . RAIN  EZ  LOOPS . 

_ D1  a  A<3301?»M<367fn  »  A  t  16482+M  I  3678  > _ 

D?  a  A  <4108?*h;3678)  -  A  (4162+M367B)  ♦  A ( 12321 +M 13678) 

_ l _ -  A  ( 1232.2  ».M  1 3678) 

D3  a  rFZB(162*HMI3h7A)  *  02 


Ain3012-*R  1 3678 )  J=  jQJ.  ---03 _ 

HO-  BA  KK=  l_t,2-4 _ 

MLE  a  R  ♦  Nfl(KK) 

_ F.i  =  al  <kki _ 

80  A(33010*MLE)  =  AR(KK) 

_ F.EEE _ 

82  CONTINUE 

F  _ LE.1-MCALL.  L  E  .  2  0 )  Gfl.  TO. .  94 _ 

C 

C _ .»...TRANS)IFRS£  El  AAiF-MO>-ft5.  .  _ 

M  =  63  *  36900 

C - 

C  . . •  .  • HX *  HZ  TRUNCATIONS . 

_ A  (.4 Q9.0  1  I4QQ.0 1  s L  A  U690 140  0  0  ) - 

A ( 36901 ♦M I  4  00  0 )  a  A f 4506 1 *M J 4 0 0 0 ) 

_ A-1A_5Q61  tM  1 4.00  Q  )  =  Ai4CL01^8t400Q.) - 

C 

A  <S72?1  tMUnnm  a  A  t  ^  32  ?1_»M14QQQ) - 

A  <53221  ♦MI40OO)  =  A  (  4  1  36  1  *M  |  4  00  0  ) 

A ( 6 1 38 1 l  4  p  0  0 )  a  A ( ?  0  3? 1 *M  )  4  0  0  0  ) _ 

_ ..  .  .  »E  X  ITERATION..... 

ASSIGN  01*. OYN . 3839 

ASSIGN  02  «  «  DYN . 3839 _ 

ASSIGN  03  » . DYN . 3839 

0)  =  A  (2^461  *413839)  »  C.A(l) _ 

02  a  A ( 1 2321 ♦HI  38  39 )  -  A ( 1 2 1 6 1 +M l 3839 )  ♦  A (20321 ♦MI3839) 

i _ «  A(57221*Mt3B39) 

03  a  C8(l)  *  02 
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_ A  (24481  *NI3839)  =  Dl  -  D3 _ 

DA  b  VABS(A(28161*Mtl59) |D4) 

_ BY  =  D4.GT.A (28461+MI159) _ 

A ( 28481 ♦Mil  59 )  =  QBvrTRL(DA»aYIA(?84ftl»M|159) ) 

_ FREE _ 

C 

C . . F.7  ITERATTON . . 

ASSIGN  D 1 ♦ » DYN .3678 

_ ASSIGN  D2«. DYN. 3678 _ 

ASSIGN  D3.. DYN. 3678 

_ AL(2)  =  0.5  »  (  A  (36m52*M)  *A  (36853+H)  ) _ 

AL(3I?1)  =  0.333  *  (A (36852+HI21) .A ( 36853*M I  ?1 ) 

_ 1 _ ♦A(36854*M|?l)  ) _ 

AL<?4)  =  0.5  *  (A  (36«73*H) ♦A(36874»M)  ) 

_ AP(2)  =  O.S  «  (  A  (36877*H) *A (36878  +  M) ) _ 

AR  ( 31  ?1 )  =  0.333  *  («  <3f.ft77*MI21)  4.A(36878*M|?D 

1 _ *A  ( 36879+M  I  ?1 )  ) _ 

AR  ( 24  )  =  0.5  *  (A  (3<.»9B*M)  ♦A(3fc89Q*M)  ) 

_ DO  91  KK  =  1,24 _ 

MLE  =  M  ♦  KK 

_ MR I  =  M  ♦  MQ(KK) _ 

A (36850+MLE )  *  A ( 32 A62*MR I ) 

9 1  A(36875-MLE)  »  A  (3300-9 R I ) _ 

Dl  =  A (3301?*MI3676)  *  CA(3) 

_ 0?  =  A (41Q62*MI3678)  -  A  ( 4 1  b2*M  1 3678 )  <•  A  ( 1 2321  ±M  I  3878) 

1  -  A  (12322*hI3678) 

_ 03  =  CP (3 )  «  02 _ 

A (3301?*MI3678)  =01-03 

_ S12..9.2  K*  =  l..i24 _ 

MLE  =  M  ♦  HQ (KK ) 

_ A  l32fl51.*MLEl  =..AL(KK) _ 

92  A (33010 ♦MLE )  *  AR(Kkj 

_ FREE _ 

C 

C _ .....VERT,  PLANF  ENVELOPE  COMPUTATION..... _ 

94  ASSIGN  Dlt. DYN. 3840 

_ ASSIGN  D2,. DYN. 3840 _ 

M  s  3?  *  38900 

_ MZ  =  65  »  36900 _ 

C 

C . . HX  ENVEl  OPF . . 

D?  =  VABS(A(4090UH|3840)  ID2) 

_ BV  =  D2.GT.A  (HZ*1  13340) _ 

A (MZ* 1 1  3840 )  =  08VCT9L(D2»BVIA(MZ*H3840)  ) 


C 


c 


. HY  ENVELOPE . 

Dl  =  0.5  »  fA(1216)*Ml3P40)*A(49061*MI3840)) 

02  s  VABS (Dl ID2) 

BV  =  D2.GT .A (HZ  +  400) :3840) _ 

A (MZ44001 13840)  =  GBVCTRL(D2*HV|A(MZ*4001I3840) ) 


P?  »  VABS( 


flAQ)  10?.) 


HZ  FNVELOPF 
A(57221*Mn 
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RV  =  D2.GT.A (MZ.8001 t3840) 

_ A  (M7.POOH38AP)  =  PRvCTRL  (D2.HVIA  (MZ  +  800H3840)  ) _ 

C 

C _ .....LX  ENVELOP?. ..... _ 

D1  =  0.5  *  (A (24481*m»3840) ♦A(fel381.Ml3840) ) 

D2_  =  VABS(DHP2) _ 

BV  a  D2.GT.A (MZ*l?0ni 13840) 

_ A(MZ*1200H38»0)  =  QflVCTPL(Pg.8VtA(MZ*l200H3B40n 

c  \ 

c . . EY  ENVELQor . . 

D?  a  VABS(A(85541*M»1840) »D2) 

BV  a  P?.GT.A(MZ*160fM  <3840) _ 

A (MZ.16001 *3840)  =  ORVCTRL (02.8VIA (MZ*16001 13840) ) 


C  . FZ  ENVELOP? . 

_ QI_s._iU5_»  -(A_(.32a51.*Mf3ft.4.Ql.tA  {fr9?51±hOM.(L)  ) _ 

D2  a  VARS(D1|02) 

BV  a  D?«GT .A (MZ+20001  13840) _ 

A(MZ*?0001 *3840)  a  OR VCTRL ( 02.BV *  A ( MZ*2000 1 » 3840 ) ) 
_ EHEE _ 


T3  a  SFCONO  (CP ) 

_ PPINT  150.  T3 _ 

C 

C . . F  IELD  ENVFI  OPE -PRINTOUT.  ROUTINE . 


DO  100  L«NHAl  F .NMAX.MHALF 

_ I£(M>Ffl^L)GO  TO  1 Q 1 _ 

100  CONTINUE 

_ _ IE(M.FQ*N8A-X1-GQ  TO  .1  (U _ 

GO  TO  199 


C  . AT  HORIZONTAL  SYMMETRY  PLANE..... 

)  0  1  IF  (N.FQ«NMAX1_IP-UN=1 _ 

PPINT  102.  N 

10?  FORMAT < 1H1 .52X.27HH7  ENVELOPE  E QR_  TIME  STEP  a. 15. 
1  //.62X.15HPI.ANE  Z  *  24*DX,//»2X*  1MJ,/) 

_ CAU.  FNVH  (2A.321  «  1  0000.0.  IPliDl) _ ; _ 

C 

_ PRINT  1-03.— N - 

103  FORMAT < 1H1 ,52X,?7HE X  ENVELOPE  FOR  TIME  STEP  *»I5. 

_ 1 _ //.1.2X. -15HPLANE  Z  =-2 4«QX  . // .  2X-«  1 HJ  .  / ) - 

CALL  FNVH (28481 .RC. IPUN) 

C - 

PPINT  104,  N 

104  FORMATMH1  .52X.27HFV  ENVELOPE  FOR  TTME  -STEP  a.TS. 

1  //,6?X , 15HP(  ANE  Z  =  24*DX,//,2X*1HJ, /) 

_ CAM  FNVH (32891  .RC. TPUN) _ 

C 

C . . AT  VFRTTCAl  SYMMETRY  PI  ANF . . 

PRINT  105,  N 

05  FORMAT  UH l.» h2.)U 21 HHX  ENVF1  QPE  FOR  TIMF.SIEP  4.15. 
1  /  /  »8?X  »  1 7HP|.  ANE  Y  =  32 .5*DX  »//»2X  ♦  1HK ,/) 

_ CALL  ENV  V  ( 1  ,  1-0  0  0  0  .  Q.»  I  PUN) - 

C 
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_ PR  I  MI — 1  flh  *.  hi - 

106  FORMAT  <  1H1  »5?X.?7HHV  ENVELOPE  FOR  TIME  STEP  =  ,IS, 

1 _ //.62X . 17PPLANF  Y  =  32  »  5-tO  X  «  X  ,  1  Hft  ,  / 1 

CALL  FNVV (4001 . 1000O.0* IPUM) 


PPINT  107,  N 

107  FORMAT ( 1H1 ,5?X,27HFI7  ENVELOPE  FOR  TIME  STEP  =  ,Ib, 
1  //»62X, 17HP|  ANE  Y  =  3?.S*nx,//»?X,lHK,/) 

CALL  FNVV(R 001, 10000,0, IPUM _ _ 


PRINT  10P.  _N _ _ 

1  OR  FORMAT  ( 1H1  ,S?X,?7HEx  ENVFLOPt  FOR  TIME  STEP  =,|tj, 

_ 1 _ //,6?X,  17HPLANE  Y  =  3g.S»OX,//,2X,lHK,/) 

CALL  ENVV(1?001,PC,TPUN) 


PRINT  109,  N 

109  FORMAT  ( IH1 ,5?X,?7HFy  ENVELOPE  FOR  TIME  STEP  s,  IS, 
1  //,f ?X,  17HP|  />NE  Y  =  32.6*DX*//*2X»1HK,/) 

CALL  FNVV ( 16001 ,RC. T»UN) _ 


-PRINT  110,  N _ _ _ 

110  F0PMAT(1H1,S?X,?7HF7  ENVELOPE  FOR  TIME  STEP  s, is, 

_ 1 _ //»6?X,  17HPL  ftNE  Y  =  3?»6<>DX  ,//,  2X  , lMK , /) 

CALL  FMVV (?0001 ,RC, IPUN) 


19Q  CONTINUE 

200  CONTINUE _ _ _ _ 

T 4  a  SECOND (CP) 

PRINT  150,  T4 _ _ 

STOP 

_ £NJB _ 

Sl'RPOUT INF  FNVH  (LOCA  ,SCALE»  IPUNCH) 

_ DIMENSION  4(24^00)  ,IP(10?40)  .nN(ISQ) 

COMMON  A 

DO  9  1=1,160 _ 

9  NN ( I )  a  I 


DO  1  LYa?,64 

_ LOC  a  |  OCA  ♦  (LY»1 )»  i6900 _ 

LOCI  a  1  ♦  (I.Y-1)*1S0 

_ IPtLOCniEO)  =  SCALP  *  A(LOC*lMQ) 

l  A (LOC? 160 )  a  0, 


IF ( IPUNCH.FO.OIGO  TO  4 

_ DO  2  LY=lb,S? _ 

LOCI  =  1  ♦  ( L  Y- 1 ) * 1  so 

_ LOCII  =  LOCI  t  1*8 _ 

?  »<R I TF  ( R , 3 )  (IP(LL)  .l.L=LOCI, LOCII) 


<=  PRINT  *.  LYY.  <  IP  (t.U  «LL  =  LOCI  »LOCI  I  > 

6  FORMAT  MX.I?.2A..251_^ _ 

7  PRINT  R.  (NN(LL) .IL=LXA,LXR> 

R  FORMAT  LLL « 4X  «  25.1  5  </.//-/  ) _ 


_ _ KE.IURN _ 

SUPROUTINf-  FNVV  (L0C7.SCALF-  .  IPOnCH) 

_ DIMENSION  A (?«■?? 500 )  ,NN ( 160) » I  PR (?5)  . IRP (160) 

COMMON  A 


DO  R  1=1.160 

R_  NN  (  I )  =  I _ 

L.MM  =  lOCZ  ♦  65*36R:V) 


IF  (  IPIINCH.FO.O  )  00  TO  3 

DO  -1  _  L7  =.  1_*  P-A _ 

LOC  =  LMN  ♦  (17-1)<M*>0 
IRP(Hl5R)  =  SCALE  «  A  (LOCi  1  6>  V  > 
1  *P  I  TE  (  R  »  ? )  (IPP  <LL  >  .LL  =  1  .  1  ^  V#  > 

?  FORMAT (16151 _ 

3  00  6  lXal.fr _ 

LX A  =  F  ♦  (LX-1)*?S 

LXh  =  _LAA _ i _ Z/t _ 

00  4  L7=l.?4 

_ LZZ  =  t^LZ _ 

LOC  =  LMN  ♦  (L/Z-1)*160  ♦  LXA-1 
_ IPR(1<?5)  =  SChLE  »  A(L0C;?5) 

4  PRINT  S.  L77.  (IPR(LL) «LL  =  l»<?b) 

5.  FORMAT  (1X.I?.?X«?5IF) _ 

6  PRINT  7,  (NN (LL ) .LL=l X A,LXP> 

7_  FORMAT t //.4X,?S IS, ////) _ 


A ( LMN I 3R4Q 
RFTUPN 
END _ 


0 


2.3  Problem  C  --  Task  2,  Case  1 
(Section  4.1  of  Volume  1)- 

The  following  10  pages  list  the  computer  program  for  the  24  x  100  x  24 
cell  --  1800  time  step  run  of  Problem  C.  The  problem  solved  is  penetration 
of  a  12.8  cm  diameter,  28  cm  long,  missile  guidance  section  by  a  300  MHz 
plane  wave  at  axial  incidence,  for  the  case  of  the  interior  dielectric  com¬ 
ponents  modeled. 
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PROGRAM  FTTH  (INPUT, OUTRUT»TAPE60=INPUT,TAPFftsTAPEB) 

C 

c  pun  TACK?—  STEADY  30 0  Mh7  PLANE  WAVE  IRRADIATION  OF  A 

r  IP. 8  CM  DIAPFTEw  MISSILE  GUIDANCE  SECTION 

C  CASE  T_  INTERIOR  OIFLFCTPIC  COmPONFNTS  M00ELF.0 

c  axial  incidence 

C  INCIDENT  wave  COMPONENTS  F7  ano  HX 

c  ?4  X  100  X  ?4  CEIL  CIJRIC  SPACE  LATTICE  IS  USED 

C  UNIT  CELL  niAKETrP  =  DX  =  0.3?  CM  =  pAVt'LENOTH/300 

r  FVFM  CYMK'FTWY  apOWT  LATTICE  PLANES  x  =  ?4.5*0X  AND 

C  7  =  PA.pooX  TS  ASSUMED 

C  SOFT  TFM  WAVF  soupcf  condition  IS  used  AT  plane  Y  *  3.0*DX 
C  SOFT  LATTICE  TRUNCATIONS  APE  USED 

C  PPOOPAM  IS  OPT  I M 1 7F.  0  FOB  TWF  CDC  STAR-100 

C 

RF  AL  VHP.  N(i7 

niNFNSTON  A  (SPG636)  »CF XR (60000 )  ,  C^YP  <*0000  )  •CFZ«(60000)  . 

I  7  (8B36)  .  AAA  (600)  »AA(?S),PD(?4)*DE(?4),CA(Q), 

?  CR(C) .EPS (6) ,SIG<6) 

COMMON  A.rFxP.CFYP.CEZP 

nFSCPTPTOP  D1 .OP.O 3,04,08. NFU»A1 .A?»mV.BW»BX,Ry 
PIT  RV.RRP (F00) .RW.RRW (600) ,1X.RRX (600) ,HY.PRY (?4) 

A c S  T OF1  PV  .PRO  ( 1  5*00  ) 

ASSIGN  PW.HPW (1 *600) 

ASSIGN  RX .pEX ( I  *  bO 0 ) 

ASSIGN  PY,RRY ( 1 *?4) 

A ASTGW  DA  ,00  <  1  *2A  ) 

ASSIGN  OS ,DF (1*74) 

T i  =  sFCONO(rP) 

PPINT  180,  T 1 
180  FOPMAT  (FP0.8) 


C 

c  . I.  NPOflLKM  PAP ARE  TEES 

FPF  0  =  3,nF*A 
OX  =  0.01/3.0 


MPP  =  6 

DATA  EPS/  1.0,  1.0,  1.0,  8.5,  4,5,  5.3  / 

DATA  SIG/  0.0,  3.7F,7,  0.025,  0.0024*  O.OOOfl,  0.0  / 

mm  a  x  =  lft  nr. 

MPA  X  a  1P00 
C 

c  . IT.  BASIC  ANO  DERIVED  CONSTANTS..,.. 

PT  =  3.1418QP6S 

M|i7  =  4.0  API  *  1.0F-7 

►  pc7  -  p.<-84F-l? 

OT  a  OX  /  6 , OE  *  ft 

rinf;  F  a  0.8  /  Fppo  /  OT 

P  a  OT  /  ?.0  /  FPS7 

PA  =  OT**?  /  OX**?  /  MU7  /  FPSZ 

wh  =  i'T  /  Ox  /  MUZ 

^C  a  1.0F+A  /  RR 

40  a  P.n  *  PJ  *  F»FO  *  OT 

IPIIM  a  0 
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DO  O 


C  A  ( 1  *  9 )  =  0. 

CR  ( 1  * «» )  =  0. 

mw  =  oavmkc. n  »?s*«w) 

MX  =  OPVMK7  (  *  *^  X  ) 

00  ?  7=1, MRP 

FAp  r  P  *  SIG(I)  /  FPS(I) 

CA(i)  =  (i.o-fafi  /  (l.r^peF) 

?  CP(I)  =  RA  /  FPS(T)  /  ( 1 ,  0  ♦  F  A  F  ) 

. I T I •  LOAD  VECTOR  o . 

. 7FPO  IMTTIAL  FI  FIPS . 

Z(H593ft)  =  o. 

A(l*5036)  =  0. 

A (593601 *  5936 )  =  0. 

0  .....TYPF  OF  MEOIIIP . 

9FAO(60.4,ENO=R01 ,EPR=?01 )  (7 (1 ) ,1=1 ,600) 

4  FORMAT  (75F1.0) 

PFAP  ( f  0  » 4  »F  ND=?0 1  «FRR=?0I  )  (Z(I)  ♦  I  =  260  1  *  3?0  0  ) 

RFAn(AP.4,FK'0=Rni  ,FPR=?01  )  (Z(l), 1=1301,1900) 
00  5  J=?,6 
JPF  L  =  (J-l)  *  5936 

5  A  ( JOFl.  ♦  1  *  cc<36 )  =  7(1*5936) 

C 

00  MS  J=7,14 

RFAP(F.n.4,tNO  =  ?Ol.FRP  =  ?01  >  (7  <  J  )  «  1  =  1  ,600  ) 

PFAO  (60  .4.FMO  =  ?01  ,FRPa?01 )  (  Z  ( I  1  ,  !=?601 

RFAD<ftn,4,FKO=?01,F9Q=?n] )  (7(1) ,1=1301,1900) 

JPFL  =  (J-l )  *  5936 
55  A  (JDEI  ♦ 1*6936)  =  7(1*5936) 

C 

00  6  JA=1S,71«7 

RFAO(f>P,4,FNr=?01  ,FPP=P01  )  (Z  (  I  )  ,  1  =  1 ,600) 

RE  AO  (f.n,4,FNP  =  ?01  ,F  PR  =  ?ni  )  (7(7)  »J=?60l»320fl) 

RFAD(ftP,4,FNP=?01 ,FPR=?0] )  (7(1), T=1 301, 1900) 

JOFL  =  (JA-1)  *  5936 
A  ( JOEL  ♦  1 1  5936  )  =  7(11 59.34  ) 

RFAD(60»4  , FNO=?0 1  , FRR=?01  )  (Z  ( I )  ,  T  =  1 ,600) 

PFAO  (60»4,FN0=?01  ,F°R=?01  )  (7(1)  ,  I  =  ?6  0  1 , 3?0  0  ) 

00  6  JP=1,6 

JOFl.  =  (  J  A  ♦  JR  - 1  )  «  5935 

6  A  (JDEI.*  1  *5935)  =  7(1*5936) 

DO  7  J=°3, 100 

JOFL  =  (J-l)  *  5936 

7  A  (JOFl. ,1*6936)  a  7(1*  5936  ) 

C 

RFAn(60,A,ENO=?01 ,FOR=?ni )  (7 ( I ) , 7=1 ,600) 
RFAn(f0,4,Fwn  =  201,t'RP  =  ?01  )  (7(1)  ,  1  =  2601  »3?00  ) 
«FAO(60,4 ,FNO=?01 ,ERR=?01 )  < 7 ( I ) , 1=1 301 , 1 900 ) 

DO  H  J=75 ,77 
JrFL  =  (J-l )  *  5936 
f  A ( JDFl ♦  1  * 6°36 )  =  7(1*5936) 
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">  ") 


r 


c 


r 


C. 


r 


C. 


C 


we  AO  <60.4*EW0=201  ,ERR=?01 )  (Z(1)*I=1301»1900) 

J  =  76 

jnp|_  =  (J-1  )  *  5936 
A  (JDFt  =  Zdl^W) 

OFA0  <6p.4,FK.n  =  ?01  ,ERP=?01 )  < 2 ( I > *  I  =  1 ♦ 60 0 ) 

PFA.n<6O,4.FNP  =  ?01*EPR  =  ?Ol>  (Z(I)  *1=2601*3200) 
J  =  79 

JDFL  =  (J-1)  *  5936 
A  (  JOF  L  ♦  1  5  t'c>  76  )  =  7(1*6936) 


PFAn(60,4.FNn=?01 ,F9P=?01 )  <7 (I) *1=1*600) 

RFAn(60,4,F^n=?01.EHP=?01)  (Z(I) *1=2601*3200) 

■if  AO  (60,4  ,FMr.  =  ?01  ,EP»=201 )  <Z(I)  *1  =  1301*1900) 

no  Q  J=6p,H4 
JPFL  =  (J-1 )  *  593ft 
Q  A  (J0EL*1  ?u>c36)  =  7(1*6976) 


RFAD(60*4*E  Ai0  =  20 1 «EPP=201 )  ( Z ( I ) *  I  =  1 30 1 » 1 900 ) 

J  =  96 

JOFL  =  (J-1)  *  697ft 

a ( JOEL* 1 s  697ft )  =  7(156936) 


9F  AO  (60.4  ,FA'0  =  /r01  ,F  PP=?01  )  (Z  (  I  )  .  1  =  1  *600) 

PFAn (60,4,FFn=?01 ,FPR=201 )  (Z ( I ) *1=2601*3200) 
PEAP(60.4.FNO=?01,E9R=?01)  (Z ( I )♦ 1=1301 . 1900) 
00  in  J=P6,«7 
JOFL  =  (J-1)  *  6q7ft 
]  0  A ( JOF|  ♦ 1 5  c  93ft )  =  7(116936) 


NFAn(60,4.EHn  =  ?01 *FPR  =  201 )  <  Z  C  r  )  ♦  1  =  1*600) 

PFAn(ftn,4.EM0=?01*EPR=?01)  (Z(I) *1=2601*3200) 

3FAO(6n,4.FMr.  =  ?ni  » ERP=20 1  )  (Z(  I)  *1  =  1301*1900) 

no  11  J=6  F  *91 
JOFL  =  (J-1)  *  59^6 
11  A  ( JOFl  ♦  1  J6C-.,*  )  =  Z(  1  *6936) 


PFAO (60,4,FfO  =  201 *F  PR=20 1 )  ( Z ( I ) » I  =  1  * 60 0 ) 

kF  AO(6O.4.FF-r=201*EPP=?01  )  (7  (1)  *1  =  2601*3200) 

J  =  9? 

jr>F|_  =  (j-1)  «  6c'3ft 

A ( JOFI  ♦ 1  *  6°  36 )  =  7(1*6936) 


. VFOIA  COEFFICIENTS . 

ACSION  Al».OYt^,600 
A9STF;*-  A?«  .i^Y^  »600 
A  9  6  I G*  NF(  .  .OYf-  .600 
A  1  =  0. 

A?  =  0. 

no  ?3  j=?.ioo 

JOFL  -  (J-1)  *  S93ft 

JOF m  -  ( J-] )  *  600 
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o  o  r>  d 


c. 


NFP  =  A  (Jf>F|.  -»1  *^00) 

00  ?0  JJ=1.MPK 
HV  =  PF0.FO.JJ 
A  1  =  OSVCTPL  (CA  <JJ)  ,F«V  t  A  1  > 

?o  a?  =  opvctrl  «r*  (  jj>  ,pv:  a?) 

A ( JOFL  *1 !f Aft)  =  A  1 
FFXP  <  JPF  1  JM'f.)  =  A? 

c 

NFL)  =  A  <  JO^L*  1301  *400) 

00  ?1  JJ=1  »*PR 
GV  =  MFP.F.r.  JJ 
A  1  =  DRVCTK  (CM  JJ)  *PVIA1  ) 

?)  A?  =  OPVCTPl.  (0M  (  JJ)  *OVJ  A?) 

A  (JOFL  *1301 (600)  =  A  1 
CFYR ( J0FP  + 1 J600 )  =  A? 

C 

MFO  =  AljrFL  +  ’hOUMlP) 

00  ??  jjsl.opw 
3V  =  MF0.FO.JJ 
A]  =  ORVCTPL (CA < JJ) «RV»A1 ) 

PP  A?  =  npVfTP|  (Ch  (J,j)  ,PV;  A?) 

A  (  J0£|.  ♦?*<■:  1  Jf  0  0)  =  A  1 

?3  CF7P( JOFm+1 SFOO)  s  A? 

FWF  F 

r 

TP  =  SFCONfMCP) 

PPTNT  ISO,  T? 

. IV.  TI^F-STFPPT NC»  l.fOp . 

00  POO  N=).k>max 
TFPM  =  SIN  (FLOAT  (N)  *P0)  *  Oh 
MOALL  =  3  ♦  TFIX (FLOAT (N)/f.O) 

IF  (PCALl.GT.33)  PC  A|_L  =  3  3 

. transvfpsf  pi. anf  mo,  i . 

. FX»  fc  7  TPIIMCATIONR.  .  .  .  . 

A  (hPfytPPb)  =  A  (  1  *4?S) 

A  (  I  ;  6?G  )  =  A  ( A  l>6?  *  4  PS  ) 

aoppmgpm  =  AtpoonfipR) 

A(?4ni*f?t)  =  A(914P*4PS) 

C  . FY  ITERATION. . .  .  . 

ASSIGN  0] . .0YN.S74 
ASSIGN  0?..0YN..S74 
ASSIGN  D3..DYN.S74 

A(]9?*)  =  O.G  *  (A (PS5] ) +A (?SS?) ) 

A ( ] 9P7  tPP)  =  0.333  «  (  A  (?SS1  *?;- ) +A  (PSSplP?)  ♦MPSS.3t?P)  ) 

A ( ] Q4Q )  =  0.333  *  (A ( PS73 )  ♦  E.0#A(PS74)) 

A(?SG1*?4)  =  A  (  1  os i j  ?4 ) 

01  =  C  A ( 3 )  #  A(1«S1«G74) 

0?  =  A(30?MS74)  -  A(3R0]*S7a)  ♦  A(S?7MS7*)  -  A(SP77*S74) 
03  =  CP(3)  *  Op 
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C  C  C  C  V-  C  O'  C  O  C  O  c  o  o 


a  ( 10=.]  J  S74  )  =  '  >1  ♦  03 

fwf  F 


. TRAASVFRSF  oi.ANtS  ?  -  )00,.,... 

00  8?  v( V -  ]  •  mC *1  L 
JPFL  =  ( JY-1  )  *  17808 
jf.fM  =  <JY-1)  *  1«00 

.  .  .  .  .8  X  1TFP AT  ION.  .... 

assign  ni , ,r>Yw.S7T 

ASSIGN  f)?.  .OYf.873 

assign  ot. .r-Yr'.t'7T 
or  to  ma=i,t 

m  =•  JDF  L  ♦  SQ3is «  (VA-1  ) 
vc  =  J0h8  ♦  800* (CA-1 ) 

. SOFT  LATTICE  TRUNCATION . 

n  =  O.G  *  <  A  (71P8*C  )  ♦  A  (718t**M)  ) 

a  (8S84*NS  PI  )  =  0.333  *  (  A  (7188+MJ?!)  ♦A(718q*v;?i) 

1  ♦A(7190*v;pi) ) 

A(8GpC«.M)  =  0  .  ^  *  (  A  (  7?09  ♦!*)  ♦  A  (  7  ?  1  0  )  ) 

4  <718P.8;?.3)  =  A(AS«8*mi?3) 

.  ,.,.N A  I  A  F  X  LOOPS . 

n i  =  a  <<3Q43*m:-.73)  *  a  ( *S8F  +  m  t  S73 > 

up  =  A(ll?lT4M;^7T)-A(e;?77*M|A7T)*An0Sl3  +  8jA7T) 

1  -A (1053«*MtS73) 

03  =  ffxp  (f  ?7*nm;*73)  *  0? 

A ( ffbb+m ; )  =  n  ♦  03 

FRFF 


. FY  ITERATION . 

ASSIGN  HI  * .0YM.S74 
ASSIGN  f>?.  .Ovr  ,87  4 
ASSIGN  03  .  .  OYN  *  574 

no  40  *A=I,1 
«  =  JOF  l  ♦  8  Q  1h»(M«l) 
vc  s  JOFM  ♦  *0(>*(VA-1) 

.....SOFT  LATTICE  TRUNCATION . 

6f7ft8?*M)  =  0."  *  (  A  (R487  +  Ci)  ♦A  (  K4H8  *M  )  ) 

4 <7883**: PR)  =  0.3^3  *  ( A (84H7»m*9?) (P488*m»2?) 

V  ♦ A (H489*v;pp) j 

A  (  7  8  >•'<  G  +  N  )  =  0.333  *  <A(88no^M)  ♦  ?.0*A  (8510**>)  > 

A  (t  4H7*»*I?4)  =  A  (  7887+M  I  ?4  ) 

C  . NAT-'  fY  L 00°S . 

,i]  =  /*(7?<‘?*nji-.7a)  *  A(78A7»v»c74) 

0?  =  A(QPf?.M;->74)-A(R8  37*c*874)*A(ll?l?4H|874) 

\  -A(11?13^MIS74) 

0.3  =  OF  YH  (f?f,  +  «t»m74|  «  02 
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o  o  o  n  n  .->  n  n  o  n  n  r>  n  no  on 


A  (7P«7*M«874) 


01  ♦  03 


C 


. fy  Fr.vFi.oPE  computation . 

04  =  cavc*pps  (  a  <7*61  ♦Mifno )  ♦•*>  *  ru  > 

OS  =  VAPS ( 04  * PS ) 

RY  =  08  •  GT  .  A  (PC'1?*M|?4) 

40  A(M81?*M*?4)  =  Q8VCTRL  (08,MY*A  (PM?*M»?4)  ) 
FPFE 


. F  7  ITERATION . 

ASSIGN  Dl,.OYN.S99 
ASSIGN  0P..DYN.S99 
ASSIGN  03..0YN.8R9 

00  SO  MA= 1 «  3 
M  =  JDFL  ♦  8936* (NA-l } 

MM  =  JOFM  *  600* (MA-1 ) 

. F  A  I N  f  z  LOOPS . 

01  =  A ( PS37  .M  t  S9Q )  »  A (916?«.m;8qQ) 

02  =  A  ( 1  0C  1  3>^sPyR> -A  ( 1  091  +A  (  3901  6Q9) 

1  -A  (9M37  +  MSt'Q9) 

03  =  CFZF.  <601  ♦MM1 999)  *  02 
A (916?*M;Fqqj  =  P]  *03 

. F7  SOFT  TFm  WAVE  SOtPCF  CONOITION . 

IF < JY.OF.?.OP.hA,LF.?)OC  TO  47 
A  (  ?  1  0  34  *  SQCj  )  r  TFPV  ♦  A  {/>  1  0  34  : ‘•QQ  ) 

. F7  FNVFLOPE  COMPUTATION . 

47  04  =  VAPS (A (9737*M»?4) *04) 

RY  s  04.OT.A (Q7H7+V*?4) 

A  ( 9787*M  *  ?4  )  =  08VCTPL  (OA,bY*  A  («7F7*Mir>4)  ) 

04  =  PPVCMPPS  (  A  (91M*M*600)  ,4>l  *04) 

08  9  VAPS ( 04  I  OS ) 

BY  =  08. GT ,A (9F1?+m*?4) 

80  A (981P.M  *  ?4 )  =  OPVCTPL (08.BY* A (9P1?*M*?4)  ) 

. MX  ITERATION . 

DO  60  MA=1 *3 
M  =  JOFL  ♦  9936*(MA-1> 

. MAIN  HX  LOOPS . 

01  =  A  (  390  1  *99U) 

0?  s  A  (  1981 +M  *899)  .a  (  iQ?f,*M{S9i,  j  .*  {  3p?f,*M*SRa) 

1  -A (9i6p.M*SVV) 

A (3901 ♦M*8Q9)  =01  *0? 

. MX  FNVELOPF  computation . 

04  r  VAPS  (A  (64764.MJP4)  *04) 

BY  =  04. GT. A (49P6.MIP4) 

A (4S?6.M*?4)  =  08VCTRL (04»BY* A US?6*M*?4) > 

04  =  OPVCMPRS ( A (3900*M*600) »PX *04) 
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I 


o  “)  r>  o  r> 


QS  =  VARS  (fHJu*) 

HY  r  OS.PT.A (4SS1*M|?4) 
f-.O  A  (4BS1  ♦*  ?  ?4  )  =  (jAVCTPl,  (OS  »Py»  A  (  4551  ♦M*?‘4)  ) 

FPF  F 

r 

r  . my  itfwatiom..... 

ASSIGN  1*1  .  .DYN.SPP 
ASSIGN  03..0YN.598 
C 

00  7  0  »*  A  =  1  ,  0 
A  »  JHFL  *  SC.'lfv*  ( MA-1  ) 

r 

r  . c0FT  LATTICE  TPl'NCAT  ION . 

A  A  (  1  )  =  n.S  *  IM11117*M)«4(11)38*M)| 

A  A  (  ?  J  ??  )  =  0.333  *  (Ani]37*<;??)»A(l])3ft  +  rtI?2) 

1  ♦ A ( 1  ) 1 39  ; ,?2 ) ) 

a  A  (  ?4  )  =  0.333  *  IA(Ult'')*M)  ♦  ?.0*A< 11160+M) ) 

AAA(IJAOO)  =  GBVXF-NO  (  AA  (  1  J?4)  *hWl  AAA  (  1»500)  ) 

A  (  1  1  1  37*M  5?4)  =  OPVOPPS  (  A  (  10M3*Mtft00)  *RWf  A  (1 1137*Mf  ?4)  ) 

. ►‘Air  my  LOOPS..... 

0]  r  A ( 1  ns  1 

0?  =  A  (Qlf- S*V  5  S9H) -A  (  p1F?»m;SmF  ) +A  (  ) 

1  -A  (G5B8*M|59fl) 

A ( 1  OS] 3*M 5^98 )  =  01  ♦  0? 

A  (  1  051  ?♦*«!  BOO  )  =  OSVCTPl  (AAA  (  1  SGOO)  ,BWJ  A  (  1  0  5 1  2  ♦  M  I  6  0  0  )  ) 

A(10Sl?+H{*nC)  =  OPVCTPL (0.0. nX IA ( ] 05 1 ?♦* f 600 )) 

. MY  f-NVElOPf  COMPUTATION . 

0 4  =  VAPS ( A ( 1 1 0S7*M|?4) ;D4) 

B Y  =  04.GT.A  <  1  1  1<S?*m;?4) 

70  A  (  1  1  1  f>  ?♦**  :  24  )  =  OBVCTPL  <04. *Y  S  A  (  1  1  ]  t?4)  ) 

SPF  F 

. H7  ITEPATION . 

ASSIGN  01 . .PYN.S73 
ASSIGN  02..DYN.573 

00  Pfi  ma=1 , 3 
M  =  JOFL  ♦  SPlG*(MA-l) 

. SOFT  LATTICF  TPUNCATIOn . 

AA(?)  =  O.s  *  ( A (5P77*M) *A (SH7F+M) ) 

AA(3*?1)  =  0.333  *  (  A  (SF*77  +  *121  )  »A  (5H7A*M»?1  ) 

1  (S«79*M»21 )) 

A  A  f  ?  i  )  =  n.S  (A(SFQB.N)  ♦  a(‘RQt)*M)  ) 

AAA(lSFOf-)  =  i)HVXPMJ(AA(1|?4)  «  hW  I  AA  A  (  1  1 60  0  )  ) 

A(5«7f*M?&)  =  OAVCMPPSIA (S?S?*M;goO)  ♦BW1A(K876*M|?4)  ) 


r 

C  . "*A  ,  ■  M7  LOOPS 

01  =  A (S?77+^;S73) 
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r>o  n  o  rt 


n  ?  =  A  (  b9*F  j  S7  3  )  -  A  (  99?  ♦  w  •  *>  7  3  )  ♦  A  (1  99  1 ♦“ l 97  3 ) 

1  -A  (  19S?*M97  <) 

A.(9?77*m»973)  =  01  ♦  0? 

A  =■  Pt-VCTOl.  (  A  *  A  (  1  tf-nn)  «  «U.  t  4  (9?9l  +M  J  *,00  )  ) 

a o  a  (Spa i  ; *- on )  =  CftvcT^t.  ( ft .  ft  *  ->x  s  a  (4?s i  ♦  *«  jaoo )  ) 
fpff 

8?  COK'TINIIF 

IF  (MOAl  L.l  F  ,3?)G0  TO  94 


•  .  *  •  .TPAf  SVFkSF  PLAFf  NO.  101 
'A  =  Q<5  *  4C-** 


. F  X, 

A (PSAP+MJFPS) 
A  (9Q37*M  J9?9) 
A (7?37*MJ6?C) 
A  (9]  9?4-MU  ?9) 
A  (8937*l»';9?9) 
A <9A37*MIP?9) 


87  TPUOrATTOMS. 
=  a  <9937*m;F.?S) 
=  A  (7?37*M:9?9) 
=  A (9?h+w 19?S) 

=  A  (A937+o:9?S) 
=  A ( 9937 ♦  9  ; 9?9 ) 
=  a <3??**m: 9P9) 


hx  itfpation . 

ASSIGN  01 « .OY 0*099 
ASSIGN  0?,  .l.yw.bPP 
01  =  A  (3901  *mS999) 

0?  =  A  (19ci->m:499)_a(19?9*n»999)  *a  (  3P?6«>  »S99) 

1  —  A  (91f'?*(k'»‘at*9) 

A ( 390 1 ♦M  t  80Q)  =01  ♦  0? 

04  =  VA99 (A (4^7^*m;?4) ;p4 ) 

BY  r  04.67.  A 

A  (4SP9+WS  ?4  )  =  08VCT9L  (04.KY!  A  (4C,?A4.M5?4)  ) 

04  s  0PVOPPS(A  (3900*M;6nO)  .t-XID*) 

09  =  VARS (04*09) 

BY  =  O^.PT.A (4991 ♦Bt?4) 

A  (499]  ♦MPA  )  =  08V6T9L  (nR.MYJ  A  (49(^1 +^;?4)  ) 

FPF  F 

C  . H7  ITFPATION . 

ASSIGN  01«.OYN.B73 
ASSIGN  rj?,.0YN.S73 

A  A ( P )  =  ft. 9  *  ( A (C977.M) ♦ A  (SB7/  ♦  *)  ) 

A A ( 3 1 ? 1 )  =  0.333  *  (A(SK77*m;?] ) «A (Stt7H+M| ?] ) 

1  ♦ a (^H7«*^ : Pi ) ) 

A  A  ( ?4 )  =  O.S  *  (A(989fl+.v)  ♦  A(9B99+M)) 

AAA ( 1 *900 )  =  08VXPN0( AA  (  1  J?4)  ,OvhaAA ( 1 JhOO)  ) 

A  (5A79+M??4)  =  0BV6MPPS  (  A  (9?Sp  +  i-(*hOO  )  .BWJ  A  (tP7G*N;?A)  ) 

01  s  A  (9?77*n>»S73) 

0?  =  A  (999MM973  )  -  A  (99?  +  M  *>73)  ♦  A  ( 1  99]  ♦*!  973) 

1  -A  (  1  99?»w»c'73) 

A  ( c?77  ♦**  ?  97  3  )  =  01  ♦  OP 

A  (4P91  ♦MJftOO)  =  O'^vCTOf.  (  A  A  A  (  1  if  on  )  ,R**A  (9P91  «M(  900)  ) 

A  (S?91  ♦  MFOO)  =  O^VCTPL  (o.0«-.x;a  ( 9?Sl  ♦•* » 9  0  0  )  ) 

FPFF 


IS 


n  n 


<uu  t n  =  sfcofmcp) 

°c>TMT  1^0.  f> 

. FJFU'  FMVFLOPF  printout  POUTINf . 

no  inn  l=np<u  f«nf‘Ax.nhalf 

if  (fv.Fn.i  >00  T it  101 
ICO  CO^TIM'F 

IF  (M. FO. NUAX  )<•<>'  T n  101 
on  TO  1QP 

r  .....AT  aOP  J  70MT  AL  FYmmFTpy  PLANF..... 

1C]  IF  (A  .FO.M-  £X  )  Io0N=l 
°P  T  NT  10?i  M 

IP?  fc^MftT  <  1  M!  •*■?*  «?7HF  7  FNVFl.OPF  FOR  TI*t  STEP  =  .IH, 
1  //.OPX  «  lFMPLANF  7  =  P4*r\X  ,//.?X*  1HJ,/) 

Oil  L  F>'V  (!««]  .PC«  TPI/N) 

r 

JPTf|T  1  04  .  F' 

104  FORMAT  (  lHl  ,c.?X,?7FPX  FmVFLOPk  f OP  TI*F  STEP  =  ,IS. 
1  //«f  ?*  «  l^Mpt.AMF  7  =  ?4*0X,//.?X.lHj,/) 

CA|  L  FK‘V  (4C?F,,:<.7  7F  ■fF.IPHIsi) 
r 

PRINT  104.  M 

10*  fofmat (J hi ,f?x »?7phy  FNVFiOPfc.  fqp  TTME  STEP  =«Ib* 

1  //,*?>  .  l^RPLAFiF  7  =  ?A*PX.//,?X.1HJ./) 

CAI  l.  FMV  (  R  ?  ?  F  *  ?  .  7  7  F  +  6.  T  P  l  P  •• ) 
r 

C  . AT  VFPTirAl.  FYMVFTPY  PLAFjF..,.. 

PPIMT  1  OF .  M 

)  0*  FpPMAT  (]H1  ,t?x  ,?7HF  7  Fr-IVFLOPF  Fpp  tImF  STEP  =  .IF>, 
1  //.F?X  ,  17HPLANF  X  =  P4.P*DX»//»?X»1HJ,/) 

fail  fmw  (  ?ftf  .  wf.  IPijm) 

r 

PRINT  10F,  N 

10P  FORMAT ( ]R] *s?X*?7hm<  ENVELOPE  Fop  TIME  STEP  =  ,IS» 
1  //.F?X,17mPLAMF  x  =  P4.F*DX,//»?X» 1PJ./) 

CALI  FNV  |4C  41  «?.77F  ♦‘‘.TPIIN) 

c 

Pp  TNT  10P.  N 

IOC  FOPMAT  UF1  «F?X  .P7PFY  ENVELOPE  FOR  TImF  STEP  =  ,IS, 

1  //fF?X«]7HPLANF  X  =  24,S*0X.//*?X,lHjf/) 

C A l  L  FMV  (Ftyf.^c,  T P 1 1 N ) 

0 

1  GO  C.OK  T  I  V|  IF 

pee,  COATIncF 

PCI  Ta  =  SFCOmO (Cg) 

PPTp.'T  ISO.  TA 
STOP 

F»  0 

G',400I,TTF  F  f-MV  (I.  nr  A,SCA|.F  *  TP.  INCH) 

.  ■  i  >  I- VO  I  ON  L  («-  go Cy1h)  ,  I  p  (  ?S  00  )  .  FA  (  ?*> ) 

r  (  '  •  >/i  o  o  >  A 
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00  A  1  =  1 
A  NN ( I )  =  I 

on  i  ly=?.1oo 

Lor  =  t.oc*  ♦  rLY-i  >  *503* 

LOCI  =  1  «•  (LY-I)*?s 

ip  (Lnn  sps)  =  SCAI  *  *  Mton?1) 

1  ailoci?*)  =  o. 

IF ( IPUNCH.AQ.O) GO  TO  4 
no  ?  t  y=? • 1  on 
LOCI  =  P  ♦  (LY-1)*?S 

LOCII  =  tool  ♦  14 

?  WPTTF (P*3)  (IP(LL) .LL=LOCI*LOCII) 

3  FOPMAT(IOTO) 

4  DO  1  I Y=?« 1 00 
L Y Y  =  10?  -  LY 

LOCI  =  P  ♦  (LYY-1 ) *?S 
LOCIT  =  LOCI  ♦  19 

*>  POINT  fr.  LYY,  (IP(LL)  ♦Ll=l-0CI»10CIT) 
A  F0PMAT(1X,I3»^X,?0I6) 

POINT  7.  (I'-MLl.)  *LL  **>♦?*) 

7  FOPMAT (//. PX»?0IM 

9FTIJRM 

FNO 
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2.4  Problem  D  --  Task  2,  Case  2 
(Section  4.2  of  Volume  1) 


The  following  12  pages  list  the  computer  program  for  the  24  x  100  x  48 
cell  --  1800  time  step  run  of  Problem  D.  The  problem  solved  is  penetration 
of  a  12.8  cm  diameter,  28  cm  long,  missile  guidance  section  by  a  300  MHz 
plane  wave  at  axial  incidence,  for  the  case  of  the  interior  dielectric 
components,  metal  components,  and  wires  modeled. 
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o|o  o|o  o|o  u]u  u|o  «j|o  |  1  !  j  I  I  |  |  |  |o  u|  I  I  I  <j|<J 


PROGRAM  FDTD ( INPUT, OUTPUT  * TAPE  6  0= INPUT, TAPtb=T APE  8) 

RUN  TASK?—  STEADY  300  MHZ  PLANE  WAVE  I RRAOI ATI ON  OF  A 
12.8  CM  DIAMETER  MISSILE  GUIDANCE  SFCTION 
CASE  II-  INTFRIOR  METAL  AND  DIELECTRIC  COMPONENTS  MODELED 
INCIOENT  WAVE  COMPONENTS  ARE  EZ  AND  HX 
24  X  100  X  4*  CELL  CUBIC  SPACE  LATTICE  IS  USED 
UNIT  CELL  DIAMETER  *  OX  s  0.33  CM  =  WA VELEndTH/300 
EVEN  SYMMETRY  ABOUT  LATTICE  PLANE  X  =  24.5*OX  IS  ASSUMED 
SOFT  TFM  WAVE  SOUPCE  CONDITION  IS  USED  AT  PLANE  Y  =  3.0*DX 
SOFT  LATTICE  TRUNCATIONS  APE  USED 
PROGRAM  IS  OPTIMIZED  FOP  THE  CDC  STAR-100 


REAL  MUR,  MUZ 

DIMENSION  A(lSESSOO)  ,Z< 15500)  , AAA ( 1200) »&UE  (600) ♦ AA (48)  , 

1 _ DO  (24)  .DE  (4B)  *QF  (48)  *  CA (9 )  *  CB (9)  *EPS(9)  *SI6(9) 

COMMON  A. 

DESCRIPTOR  D1*02.03*D4«P5*P6*NFD*AI * A?*BU*bV .BW *BX ,BY *87 
BIT  BU.BBU ( 60  0 ) »BV » BBV (1200) .HW , BBW ( 1 200 ) »BX *8BX < 1200 > . 

1  BY  «BBY ( 24 ) *  B7  *  BBZ (48) _ 

ASSIGN  8U,PBU<1»600) 

ASSIGN  BV  *BPV (111200) _ 

ASSIGN  BW«BPW ( 1 1 1200) 

ASSIGN  BX  » BBX (111200) 

ASSIGN  BY,8BY<1124) 

ASSIGN  BZ»BBZ(H48) _ 

ASSIGN  D4,DD<1«24) 

ASSIGN  D5*0£ ( 1 148) _ 

ASSIGN  06.DF (1*48) 

T 1  =  SECOND ( CP ) 

PRINT  150*  Tl 

150  FORMAT  (E20. 5) _  _  _ 


. . I.  problem  parameters . . 

FPEQ  =  3 . OF 

DX  =  0.01/3.0 _ 

MPR  =  9 

DATA  EPS  /l.Ot  l.Ot  1.0*  5.5*  4,5*  5.3*  1.0*  1.0*  B.O/ 

DATA  SIG  /0.0*  3.7E*7*  0.025*  0.0024,  0.0008,  0.0*  6.6fc*6, 

1 _ 3.7E»7*  0.01/ _ 

NMAX  =  1800 


.....II.  BASIC  AND  DERIVED  CONSTANTS 

PI  «  3,14159265 _ 

MUZ  =  4.0  *  PI  *  1.0E-7 

EPSZ  =  8.854E-12 _ 

DT  =  DX  /  6 , OE  *8 

NHALF  =  0,5  /  FPFQ  /  DT _ 

P  *  07  /  2.0  /  EPSZ 

PA  a  OT»»2  /  DX»»2  /  MUZ  /  EPSZ _ 

PB  =  OT  /  OX  /  MUZ 


olo  I  jo  o  n 


C 


an  =  2.0  *  PI  *  FREQ  *  DT 
I PUN  =  0  _ 


CA(HQ)  =  0. _ 

CB ( 1 1  9)  =  0. 

b w  =  obvrko  ( i  .gstnm _ 

ex  =  0BVMK7 (24,25*RX) 

.  no  ?  i=i, RPR  . 

EAF  =  R  *  SI gTI  )  ''  EPS  ( I ) 

C A ( I )  =  (1,0  -  E AF )  /  (1,0  ♦  EAF) 
?  CB  ( I )  =  RA  /  EPS(I)  /  (1.0  ♦  E  AF ) 


. III.  LOAD  VECTOR  A. 

. ZFRO  INITIAL  FIELDS 

Z ( 1 l 15500 )  *  0. 

A(ltl5S00)  =  0. _ 

A ( 1550001 *15500)  =  0. 


.....TYPE  OF  MEDIUM . 

BEAD  (60,4.ElND  =  ?01  ,ERR=201  )  ( Z ( 1 ) , 1=1 , 600) 

4  FORMAT ( 75F 1.0) 

REA0(60.4.END=201 ,EPR=2P1 )  (Z  <  I ) ♦ I =50 0 1 , 560 0 ) 
RFAD(60,4,EMD=201 ,ERR=201 >  (Z(I)  ,1=2501.3100) 

DO  5  J=2.6 _ 

JOEL  =  (J-l)  *  15500 

5  A (JDEL+1* 15500)  =  Z(lllShOO) _ 


00  55  Js7 .14 _ 

READ (RO.A.EMD  =  201, ERR«?01)  ( Z  ( I ) , 1  =  1 ,600 ) 
REAO(6Q.4,E  ND  =  ?0 1,FpR=201)  (Z(I) ,1=5001,5600) 

RFAD(60.4.END=201,EPR=201)  (Z(I) .1=2501,3100) 

JDEL  =  (J-l)  »  15500 _ 

55  A  ( JDEL  ♦  1 1  1550  0)  =  Z(1US500> 


DO  6  JA  =  15.71 .7 

R£AD(60,4,END=?01.ERR=201)  (  Z  (  I )  , 1  =  1 ,600) 

RFAD(60,4,END=201,ERR=?01)  ( Z ( I ) ,1=5001.5600) 

PEAD(60,4,END=201 ,FRR=?01)  (Z(I) ,1=2501,3100) 

JDEL  =  (JA-1)  *  15500 

A(JDEL*1*15500)  =  Z(1*I5500) _ 

RE  AD ( 60 ,4 , END  =  ?0 1 ,ERP  =  ?01 )  <Z(1),I  =  1,60  0) 

READ(60,4,END=20I ,ERR=2Q1)  ( Z ( I ) , 1=50 0 1 , 560 0 ) 

DO  6  JR=l . 6 

JDEL  =  (JA.JB-l)  «  15500 _ 

6  A (JDEt +1*15500)  =  Z(l*15500) 

_ Q£.  L  J.=  q3j_LQ0 _ 

jrf.L  =  (J-l)  *  15500 

■I—&  ( J.D.I..I  tJL  11& SO 0  )  z.  2  {Li  L55MJ _ 

READ(60,4,END=201,ERR=201)  ( Z ( I ) , 1=1 ,600 ) 

PE  AD (60  *♦, END=20 1 ,ERR*?0 1 )  (Z(I) ,1=5001,5600) 

PFAD(60,4,E  MD  =  20 1 «ERR=201 )  (Z (I ) .1  =  2501,3100) 


00  8  J*75*77 

_ JOFI  a  (J-l  >  ♦  155QQ - 

8  A (JOEL* 1 1 15500 )  =  ZUI15500) 

C _ _ _ 

READ(60»4,END=201,ERR=201>  (Z ( I ) . 1=2501 *3100 ) 

_ J  »  78 _ _ _ 

JOEL  *  (J-l)  *  15500 

_ Af.inn  jJJLLSSflJLL  *  (155110) - 

C 

_ RF  AD  (  61)  ,  4  «  E  ND*2Q 1  §  ERR  a?  0  1  )  (Z(I) *T* 1.600) _ 

READ(60*4,END*201 *ERR=201 )  ( Z ( I ) *  1*500 1 .5600 ) 

_ J-g-Z.9 _ 

JOEL  *  (J-l)  •  15500 

_ A  f.lDFi  till  550Q1  *  7MI1550.Q) - 


C 

RFAn<An.4.FMn=2Q1  « F.RR*2D D — (Z  (  I  )  t  T  =  1  tfefl-Ol 


READ(60*4*EnD=201 .ERR *2 01) 
READ (60.4, END=201 ,ERR*201 ) 

(Z  (I) ,1*5001*5600) 

(Z(I) ,1=2501*3100) 

DO  9  J=80*64 

JOEL  *  (J-l)  *  15500 

C 

9  A (JDEL* 1 1 15500 )  =  Z(1115500) 

READ(60,4,END=201,ERR*201 > 

J  *  85 

(ZTI) ,1*2501*3100) 

JDEL  *  (J-l)  *  15500 

A (JDEL* 1 115500)  *  Z(1*15500) 

C 

READ(60*4,END=201 .ERR *201 > 

(Z(I), 1*1*600) 

REAO(60,4,ENO*H01,ERR*201> 
READ (60*4, END *201 «ERR*2 01) 

(1 (I ) *1*5001 *5b00) 
(Z(I)*I*25 01*3100) 

DO  10  J=86,87 

JDEL  *  (J-l)  *  15500 

C 

10  A  (JDEL*  1 115500)  =  ZUI15500) 

READ (60*4, END =201, ERR *2 01) 
READ(60*4,ENO=201,ERR*201  > 

(Z(I) *1*1*600) 

(Z(I) *1=5001*5600) 

READ (60  *4 ,END*201 «ERR*201> 
00  11  J*88*91 

(Z(I) ,1*2501*3100) 

JDEL  *  (j-l)  *  15500 


11  A(JDEL*1»15500)  a  ZUU5500) 


READ(60*4,END=201,ERR»201>  (Z ( I ) *  1*1 *800) 
READ(60*4*END*201*ERR=201)  (Z (1 ) ♦ 1*5001*5600) 

J  *  92 _ 

JDEL  =  (J-l)  *  15500 

A  ( JOEL*  1  (  15500 )  »  Z  ( 1  >  15500  >  _ ______ 


C  . PARTIAL  SYMMETRY  ABOUT  Z  *  24*0*UX 

DO  12  J*2* 100 

_ JOEL  «  (J-l)  •  15500 _ 

KMAX  *  49 

DO  1?  JA*1 *  3  _ 
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JR  =  JOEL  ♦  ( JA-1 ) *2500  ♦  1 

BUF(1I600)  «  A ( JB  >600) _ 

BU  =  RUF (11600) .GT.A.5 

BUFdlfeOO)  »  QBVCTRL(3.0«BUIBUF  (11600)  ) 

IF { JA.EG.3) KMAX=48 

A  ( JB+600I25)  =  BUF  ( 576 1?5 ) _ 

DO  1?  K«1 *  24 


_ KA  =  (K-l)  ♦  25 _ 

KR  =  (KMAX-K)  *  25 

1?  A(JB*KBt25)  =  BUF ( K A* 1 1 25 ) _ 

C 

C . . NON-SYHMETRIC  GROUND  WIRE 

DO  13  J=43,57 

_ JDEL  =  (J-l)  »  15500 _  _ 


A (JDEL+3199)  «  7.0 
A(J0EL+3224)  =  7.0 
13  A(JDEL*5699)  *  7.0 


00  1*  K*1R»2R 
KDEL  =  (K-l)«?5  ♦  24 
A(JDEL-»2500*KDEL)  *  7.0 
A (JDEL*5000*KUEL)  a  7.0 
14  A (JDEL +20500* KDEL )  *  7.0 


C  .....MEDIA  COEFFICIENTS 

_ ASSIGN  Alf.DYN.1200 _ 

ASSIGN  A2..DYN.1200 

_ ASSIGN  NFD, .OYN.120Q _ 

A 1  =  0. 

_ A?  =  0. _ _ 

DO  21  J=?,100 

_ JDEL  a  (J-l)  «  15500 _ 

DO  21  JA=1»3 

_ JB  =  (JA-1) *2500  ♦  1 _ 

JC  =  JB  ♦  1250 

_ NFD  x  A ( JDEL*JB< 1200  ) _ 

DO  20  JJ=1*MPB 

_ BV  «  NFO.EQ.JJ _ 

A]  =  08VCTRL(CA(JJ) .BV1A1) 

20  A?  =  QftVCTRL (CB ( JJ) .BVIA2) 

A (JDEL* JR 1 1200)  »  A1 

21  A { JDEL* JC 11200)  =  A2 

FREE  — ~ 


T?  =  S  ;:OND(CP) 

_ PR I M  150.  T2 _ 

C 

C_ _ _  _ IV,  TIME-STEPPING  LOOP 

00  2)0  Nsl.NMAX 

_ TEPM  =  SIN (FLOAT (N) »RD)  »  KB 

MCALL  =  3  ♦  I  FIX (FLOAT (N) /6.0) 
_ IF (MCALL. GT.33)MCALL«33 _ 
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O  Cl 


C 

C 


C 


C 


L 


C 


C 

c. 


. TRANSVERSE  PLANE  NO.  1 . 

. EX,  EZ  TRUNCATIONS . 

A(7S0H1?50)  =  A  ( 1 1 1250  ) _ 

A  ( 1 1 1  250 )  »  A  <23001  I  1250) 

A<10?om?5Q)  =  A(12S1I1ZS0) _ 

A (1251 1 1250)  =  A(25701U?50) 


. EY  ITERATION . 

ASSIGN  D1 , .DYN. 1174 _ 

ASSIGN  0?,.DYN.U74 

ASSIGN  D3 « .OYN. 1174 _ 

A(8B261  a  0.5  ♦  (A  (10076)  ~»A  (10077)  ) _ 

A(fl827l22)  =  0.333  *  ( A (1 0076 < 22) *A ( 1 0077 * 22 ) *A ( 1 007b 122 ) ) 

A  1 8840 )  =  0.333  *  (  A  (  10098) »2.0*A ( 10099)  ) _ 

A ( 1 0076 1 24 )  =  A (8651 124) 


A ( 10026 )  =  0.5  *  (A(10101)*A(10102)) 

A ( 1 0027 122 )  s  0.333  *  (A ( 10101 122) ♦  A ( 1010212?) ♦  A ( 10103122) ) 

A ( 1 004R )  =  0.333  *  (A(10123)*?.0*A(10124) ) 

AXLQXO-l  124  )  a  A(  1000H24) _ 

01  =  r.A  (3)  «  A  (  88S  1  111741 _ 

D?  =  A ( 1 1551 1 1 1 74) -A ( 1 15261 1 174) ,A ( 14201 ♦ 1 1 74 ) -A ( 14202 1 1 174 ) 

01  =  CB(3)  »  02 _ _ 

A ( 8851 11174)  =  01  ♦  03 

FREE _ 

. . TRANSVERSE  Pj.JttlES  2  -  UHL . . 

DO  82  JY=1,MCALL 

JOFL  =  (JY-1)  *  46500 _ 


EX  ITERATION. 
ASSIGN  Dl,. OYN. 1173 
ASSIGN  02, .OYN. 1173 
ASSIGN  03, .OYN. 1173 
C 

DO  30  MA*1,3 
M  s  JDEL  ♦  15500* (MA-1 ) 


SOFT  LATTICE  TRUNCATIONS 


A (23002*M)  = 

A (23003*N»?1 ) 

0.5  *  (A(242$2*M)«.A(2A253*M)  ) 

*  0.333  *  (AC24252*MJ21)4.A(242S3*M21) 

1 

A (230?4*N) 

s 

♦A(242S4»MI21) ) 
0.5  *  (A (24273*M)*A(24274*M) ) 

c 

A (2425?*H»23) 

=  A (?3027*M»23) 

A(2420?*N)  = 

A  (24203*^ 121) 

0.5  *  < A (2427?*M) ♦ A (24278+M) >  ' 

*  0.333  *  (A<242Y7*M|21),A(24278^M»21) 

1 

A(24224*H) 

s 

♦ A ( 24?79*M | 21 ) ) 
0.5  •  (A (24298*M) *A (24200) ) 

A (?4277*MI?3>  =  A <?4177*MI23) 

C 

C  . main  ex  loops . 

Pi  =  A<15b?7*H»1173)  *  A (23027*Mlll73> 

0?  s  A<29702**U173)-A<14202*P»1173)«-A<28352*MI1173> 

1 _ -A(28377»MI1 173) 

03  =  A<16777*F»1173)  *  D? 

_ A(?3027»HI1173)  =  01  »  03 _ 

C 

C  . EX  ENVELOPE  COMPUTATION..... 

04  =  VABS  (A  (23<S01*M»24)  »04) 

BY  =  04.GT.A(24301*MI24) 

30  A (?4301*M*24)  *  U8VCTRL  (04 , tY  » A  ( 2430 1  24 )  ) 

FPFF 


C . . EY  ITERATION 


ASSIGN  D1*.0YN.1174 
ASSIGN  Q2f.DYN.1174 
ASSIGN  D3,.0YN.1174 


00  40  MA=lt3 

M  =  JOEL  »  15S00»(MA-1)  _ _ 

C 

C . . SOFT  LATTICE  TRUNCATIONS . . 

A(243?6*M)  =  0.5  *  ( A ( 25576*M) ♦A(25577*M) ) 

_ A  ( 24327+M I 22 )  =  0.333  »  ( A ( 25576* M 1 22) ♦A(2b577*MI22) 

1  ♦A(25S78»M*22> ) 

_ A(?434Q*M)  s  0.333  »  (A(?S598*M)*2.0*A(25599+M>  ) _ 

A ( ?5576*M 1 24 )  *  A (24351*M124) 


A ( 25526+M )  =  0.5  *  <A<25601*M)+A(?5602*M) > 

A (25527+M122)  =  0.333  •  ( A ( 2580 1 *M I 22 ) »A ( 25602+M 1 22 ) 

1  ♦A(25603*MI22) ) 

A ( ?5549*M )  =  0.333  ♦  ( A <2b623»M) »2.0»A (25624+M)  ) _ 

A  (25601  *M24)  =  A  (2550 1  *M 1 24 ) 


C  . MAIN  EY  LOOPS . 

_ 01  =  A(180 26 11174)  «  A ( 2435 1 *M 111 74  ) _ 

0?  =  A(27051*N|1174)-A(?7026*NI1174)*A(297014MU174) 
1  -A(Z9702*MIU74) 

03  =  A  ( 19276  +  N ;  H74)  *  D2 

_ A  ( 24351 *M  < 1174)  g  01  ♦  03 _ 

c  “  . . "  ~ .  ' 

r . . FY  F NVELOPE  COMPUTATION . . 

04  =  VABS (A (24926*MI24) 104) 

BY  =  04.GT.A (25628»M|24) _ 

A  (?S6?f>->*»24)  =  QBVCTPL  (04*  BY  I  A  ( 25626  +  M 124 )  ) 


OS  =  08VCMPRS  (A  <24325*Mn  200)  »BXID5) 

06  =  V ABS ( 05 1  08 ) _ 

RZ  =  08.GT.A <2b651*M»48) 

40  A  (25551  •»M|4B)  =  QflVCTRL  (P6«  BZ I A  (25651  ♦MI46)  ) 


c 


c 

. EZ  ITERATION,.... 

ASSIGN  Dlt.0YN.1199 

ASSIGN  02t.DYN.H99 

ASSIGN  D3,.0YN.1199 

C 

00  50  HAal,3 

c 

M  a  JOEL  ♦  15500* (MA-1 ) 

c 

.....MAIN  EZ  LOOPS..... 

01  =  A(20501*M»1199)  *  A ( 257 0 1 *M *1 1 99 ) 

02  a  A(?8352*MtU99)-A(?ft351+Mlll99)*A(1152o*f4|U99) 
1  -A(27026*M(1199) 

D3  =  A  <21751 *MI 1 199)  *  D? 

A(25701*M|1199)  »  D1  ♦  03 _ 

C 

C 

. EZ  SOFT  TEM  WAVE  SOUKCE  CONDITION . 

IF<JY.GE.?.0R.MA.LF.2)G0  TO  47 

A  (56701  ( 1 199)  =  TERM  ♦  A (5670 1*1199) 

C 

C 

. E7  ENVELOPE  COMPUTATION . 

47  A ( P6926*M 1 24 )  =  0.5  *  ( A ( 26276*M » ?4 ) ♦ A ( 2630 1 *M 1 24 ) ) 
04  a  VABS (A (?6926*MI24) ID4) 

BY  *  04 «GT • A ( 26951 *M| 24 ) 

A ( 26951 *M 124 )  a  Q8VCTRL (04. BY < A ( 2695 1 *M * ?4 ) ) 

c 

_ 05  r  QBVOPRS(A<257QO*MH?QO)  »hxt05) 


06  *  VABS (D5I06) 

B7  =  06.GT.A(26976*M|4B) 


c_ 

50  A ( ?6976*M 148)  a  Q8VCTRL (06 t BZ * A (26976*M 1 48 ) ) 

C 

. HX  ITERATION,.... 

DO  60  MA= 1.3 

c_ 

M  a  JOEL  ♦  15500* (MA-1 ) 

c 

. MAIN  HX  LOOPS . 

_ D1  a  A ( 1 1526»M (1199) _ 

0?  a  A(8851»MI1199)-A(8826*M»1199)*A(10201*M|1199) 

_ 1 _ ■A(?5701*M1119Q1 _ 

c_ 

A ( 1 15?6*MI 1199)  a  01  ♦  02 

C 

. HX  ENVELOPE  COMPUTATION . 

_ A  ( 1 275 1 *M (24 )  a  Q.5  »  ( A ( 1 2 1 0 1 *M 1 24 ) + A (12126* M 1 2  4 J) 

04  a  VABS(A(12751*MI24) 104) 

_ BY  s  DA-.GT .At  1,2.7 7 6 f M.L2.4 ) _ 

A(12776*M(24)  =  Q8VCTRU04»BYtA<12776*MI24)  > 


D5  a  QBVCMPRS (A (115?5*M»1200> *HX(05> 

06  =  VABS  <05 1 06 )  _  _ 

BZ  »  D6.GT.A ( 12801*M|48> 

60  A (12801 ♦HI  48)  =  Q6VCTRL  ( 06 ♦ BZ ( A ( 1 ?80 1 *M 1 48) ) 
FREE 
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ASSIGN  02t.0YN.1198 


DO  70  MA= 1 ♦ 3 
M  =  JOEL  ♦  15500* (MA«I ) 


. SOFT  LATTICE  TRUNCATION . 

A  A ( 1 )  =  0.5  *  (A(29601*M) .A (29b02*M> ) 

AA(2I46)  =  0.333  »  (A(29601*HI46)*A(29602«MI46) _ 

1  ♦  A (29603*M|46> ) 

AA<4«)  =  0.5  *  (A (29647tM> tA (29648tM) ) 

AAA(1 11200)  =  Q8VXPN0  (AA  ( 1148)  ,BWI  AAA  ( 1  U200)  ) 

A (29601 ♦MI48)  =  Q8VCMPRS (A (28352.M* 1200) f 8«IA (29601*MI48) ) 


. MAIN  NY  LOOPS..... 

D 1  =  A (2R352.MI 119P) 

D?  a  A(25702*MIH98)-A(2570l*Mlll9fl)*A(23002*M|1198) 
1  -A (23027.MI 1 198) 

A  ( ?835?*M <1198)  a  01  ♦  02 


A (28351*MI1200)  =  08VCTRL <AAA(1 11200) »BWIA(28351.MI1200) ) 
A (28351 ♦Mil ?0  0)  =  Q8VCTRL ( 0 . 0 , BX I A ( 2835 1 1 1 *00 ) ) 


. ....hy  Envelope  computation . 

A ( 29576 +M 124)  a  0.5  *  < A ( 28926  +  M » 24 > *A ( 2895 1 tM 1 24 ) ) 
04  a  VABS (A (29576>MI24) ID4) 

hY  =  D4.GT.A(29651^MI24) _ 

70  A  (2965 1  *M  1 24 )  =  GI8VCTRL  (04  »  BY  I A  (2965 1  ♦M 1 24 )  ) 

FREE 


( 


. . HZ  ITERATION 

ASSIGN  Dlt.DYN.1173 
ASSIGN  02..0YN.1173 


QQ  ftp  MA=1«3 _ 

M  a  JOEL  ♦  15500* (MA-1 ) 


. SOFT  LATTICE  TRUNCATION,,... 

AA ( 2 )  =  0.5  »  (A(15427*M).A(154?8.M)  ) _ 

A A ( 3 1 45 )  =  0.333  *  ( A { 15427 »M!45) *A ( 15428*M|45) 

1 _ ♦A(15429^M|45)  ) _ 

A A ( 48 )  =  0.5  *  (A(15472^M)*A(15473*M) ) 

AAA  (111200)  =  Q8VXPND (AA(l|48)tBWIAAA(lll200)) _ 

A  ( 1  5426  *M  1 48 )  *  G8VCMPRS  ( A  (  14177  ♦**  11200)  *  BW  |A(1 5426  .M 148)  ) 


. MAIN  HZ  LOOPS..... 

01  a  A ( 14202*M| 1173) _ 

02  a  A(23027*MI1173)-A(7527*MI1173)*A(8851*MIH73) 

1 _ «A(8e52*M|1173) 

A(14?0?*Mni73)  =  01  ♦  0? 


46 


IO  OI  I  OI  I  oio  |o  OI  I  I  |nn 
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A(14176*M»1200)  =  08VCTRL (AAA (1 11200) tBWl A ( 14174*M»1200) ) 

A(14176*M 11200)  =  Q8VCTRL<0.0,exiA(14176*H»1200) ) 

. . HZ  envelope  COMPUTATION . . 

04  =  VABS(A(14776*M124)  104) 

by  «  Q4, GT .  A ( lb476»M|24) _ 

80  A { 1 5476  +  M I ?4 )  =  Q8VCTRL ( 04  * B Y 1  A ( 1 54 76*M » 24 )  ) 

FREE _ 

82  CONTINUE 

IF (MCALL.LE.32) 00  TO  94 _ 


. . TPANSVERSE  PLANE  NO,  101 

M  =  99  *  15800 


. EX,  EZ  TRUNCATIONS . 

A (2300l*Mt 1250)  =  A (1 55 0 1 U 2b0 ) 
A  < 1 550 1 ♦M 1 1 250 )  *  A(16751*M112bO> 
A  ( 16751  ♦Ml  1250  )  *  A(7501*MH2bP) 


A  (25701+m  1 250  )  =  A  ( 1 800 1  *M  tl  25Q  ) 
A ( 18001*M*1?50)  =  A (19251*Mll2b0) 
A  ( 1  9251  1 1 250 )  *  A  (10201+MH250) 


_ HX  ITERATION... .. _ 

ASSIGN  D1 » . DYN, 1 199 

ASSIGN  02 , , OYN « 1 1 99 _ 

OI  =  A ( 1 1526*Hl 1 199) 

D2  a  A  (8851»**11199)» A(88?fe*M|lI  99 )  »A(1020l*Mtll99) _ 

1  -A{25701^m;hq9) 

A ( I 1528+HI 1199)  =  01  ♦  D2 _ 

C 

_ A  (12751** 124)  s  0. 5  *  (A(12101»M124) »A(1212b»Ml?4)  ) _ 

04  =  V  A  6  S ( A ( 1 2  7 5 1 ♦ M 1 2 4 )  104) 

_ BY  =  D4,GT.A(12776*HI24) _ 

A ( 12776*^124)  *  Q8VCTRL <0A»BY1 A ( 1277b*M»24) ) 

C _ 

OS  =  Q8VCMPRS (A(11525*MI1?00) ,  H  X  t D5 ) 

_ Dfe  »  VAPS (05  I Ufe  ) _ 

HZ  =  OG.GT.A ( 12801,M|48) 

_ A  (  1  380  1  148)  s  Q8VCTRL (D6,BZ t A  (12B01*M4d)  ) _ 

FREE 

c _ 

C  . HZ  ITERATION..... 

_ ASSIGN  01, .DYN. 1173 _ 

ASSIGN  02, .OYN. 1173 

_ A  A  1 2 )  e  0.5  «  (A(15427,H)*A(lb428*M)) _ 

AA ( 3 1 45 )  =  0.J33  *  ( A (15427*M»45) *A ( 15428*m|45) 

1 _ *A(15429*M145)  ) _ 

A A ( 48 )  *  0.5  *  ( A ( 1547?*H> ,A { lb473,M) ) 

_ AAA  fl  1  1200)  =  Q8VXPND ( AA ( 1 1 48 )  , BW I  AAA ( 1  I  1 20 0 )  ) _ 

A(15426*M|48)  *  08 VCMPRS ( A ( 1 4 1 77*M fl 200 ) , BW I  A ( 1 5426,M J 48 ) ) 

r 
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D 1  =  A(14?0?*HJ1173) 

- Q?  =  A  (  230-27*  M  LI  1731-A(7527*H11173)  ♦A(fl851*Mtl  173) _ 

1  4A (88&2*M»11?3) 

_ a(14?q?»m;ii73)  =  D1  ♦  n? _ 

A  ( 14176*M51200)  =  Q8VCTPL  (AAA  ( 1  U200)  »HWI  A  (14176*rt|  1200)  ) 
_ A(14176»M-LL2PQ)  «  Q3VCTPL  (  0 . 0  ,  PX  >  A  ( 1 4 1 76*M  1  1 200  )  ) _ 

_ H4 — =_  VABS(AM4776  +  MI?4>  104) _ 

FY  =  D4.GT.A ( 1S476.MJ24) 

_ A-L15476  +  ML24)  =  OBV CTRL ( D4 . BY t A ( 1 5476+M I 24 )  } _ 

FREE 

94  T 3  =  SF CONI')  { CP ) 

PRINT  150.  T3 


....’.field  envelope  printout  routW 
00  100  l=nhalf.nmax .nhalf 
ifTnTeqTlIgo  To  ioi 
100  continue 

if<n.eq.nmaxT5o  to  101 
GO  TO  199 


. AT  HORIZONTAL  OBSERVATION  PLANE . 

101  IF (N.FO.NMAX) IPUN=1 
PP I  NT  102.  N 

10?  FORMAT (1H1.52X.27HEX  ENVELOPE  FOP  "TIME  STEP  =  ,I5. 
1  //.6?X, 15HPLANE  Z  *  24*DX ,//,2X, 1HJ,/) 

CALL  ENV (8801 .20.RC.IPUN) 


PRINT  1 03.  N 

103  FORMAT ( 1H] .52X.27HEY  ENVELOPE  FOR  TIME  STEP  =  ,I5, 
1  //.62X,  15HPLANE  Z  ^~24*t)X » //♦  2X » 1HJ. /) 

CALL  ENV ! 10126, 20. RC. IPUN) 


PRINT  104,  N 

104'TORMATTThTT52X,27HEZ  "ENVELOPE  FOR  TIME  STEP  *,15." 
1  //.62X.15HPLANE  Z  =  24*DX , //» 2X » 1HJ, / ) 

CALL"  ENV(11451,20,RC.IPUN) 


PRINT  105,  N 

105  FORMAT ( 1H1 .52X.27HHX  ENVELOPE  FOR  TIME  STEP  *,I5, 
1  //.62X.15HPLANE  Z  *  24*DX,//»2X» 1HJ,/) 

CALL  FNV (l?776,20,3.77£»fe. I PUN) _ 


PRINT  106,  N _ 

106  FOP MAT (1H1.52X, 27HHY  ENVELOPE  FOR  TIME  STEP  =,15. 

1 _ //.62X.15HPLANE  Z  3  24»DX,//,2X,  1HJ,/) 

CALL  FNV(14151,?0,3.77E*6»IPUN) 


PRINT  107,  N 

107  FORMAT (1H1.52X.27HHZ  ENVELOPE  FOR  TIME  STEP  =,I5, 
1  //,62X, 15 HP LANE  7  *  24*DX »//, 2X ♦ 1HJ, / ) 

CALL  FMV< 1 5476, 20, 3. 77E *6, I PUN) 
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_ .  , AT_VtRT  I  CAL  SYMMETRY  PLANE . 

PRINT  108.  N 

1QP  F PPM AT ( lHj  »5?X»?7HE7  ENVELOPE  FOR  TIME  STtP  =  ,15. 
I  //»6?X, 17HPLANE  X  =  2 4 ,5*DX ,//, 2X . 1 KJ , / ) 

CALL  FNV(11476«40,RC,IPUN ) 


PRINT  109,  N _ _ _ 

109  FORMAT ( 1H1»5?X,27HHX  ENVELOPE  FOP  TIME  STEP  =  ,I5, 

1 _ //«6?X,  17HPLANE  X  =  C4.5»DX,//,2X»lHJ,/) 

CALL  FNV(1?801»40.3.77E*6.IPUN) 

CUR  1  =  A (»9«15S00»1017B)  »  RC  »  73.33 _ 

CUP?  =  A  (76*15500*10161 )  *  RC  *  73.33 

PRINT  150,  CUR  1 _ 

PRINT  150*  CUR? 


PRINT  110*  N 

no  format <ihi»5?x«?7hey  envelope  for  time  step  =,ib. 
1  //»A?X* 17HPLANE  X  =  24.5*DX»//»2X»1Rj,/) 

CALL  ENV ( 10151 .40. RC, IPUN) 


IRQ  CONTINUE _ 

?00  CONTINUE 

?01  Ta  =  SECOND (CP ) 

PRINT  150  *  T4 

STOP _ 

FNO 

SUBROUTINE  ENV ( LOCA , NUM , SCALE ♦ I PUNCH! 
DIMENSION  A (1565500) , IP (4000 ) .NN ( 50 ) 
COMMON  A 

DO  1  I =TV50 
1  NN ( I )  a  I 


DO  2  L Y=? *100 

LOC  =  LOCA  ♦  ( L Y- 1 ) *  1 5500  ♦  4 

LOCI  =  1  ♦  (LY-1)»40 _ 

I P  <  LOC  1 1  NUM )  =  SCALE  *  A(LOCINUM)' 
2  A(LOCINUM)  =  0. 


IF ( IPUNCH.EQ.O) GO  TO  5 
DO  3  LY=2,100 
LOCI  =  1  *  <LY-1)*40 
LOCI  I  =  LOCI  -  1  ♦  NUM 

3  WRITE (8,4)  ( I P  <  LL ) *  LL*L0C I ♦  LOCI  I ) 

4  Format  (ToT6l 

5  LXM  =  NUM/?0 
DO  8  LX=1.LXM 

LXA  a  1  ♦  (LX-1 ) *20 
LA  =  LXA  ♦  4 

LZ  =  LA  ♦  19  " 
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DO  6  LY=?*100 

LYY  =  102  -  LY 

LOCI  =  LX A  ♦  (LYY-1 >  *40 

LOCII  =  LOCI  ♦  19 

6  PRINT  7*  LYY*  ( IP <LL) *LL=LOCI »LOCII ) 

7  FORMAT  {1X«I3*?>X*20I6) 

P  PRINT  9*  (NN(LL) »LLsLA»LZ) _ 

9  FORMAT (//*PX*20I6»//////) 


RFTURN 

END 


2.5  Problem  E  --  Task  5,  Case  1 
(Section  7.4.1  of  Volume  1) 

The  following  10  pages  list  the  computer  program  for  the  24  x  163  x  24 
cell  --  800  time  step  run  of  Problem  E.  The  problem  solved  is  a  hybrid  meth¬ 
od  of  moments/FD-TD  analysis  of  penetration  of  a  19.0  cm  diameter,  68.5  cm 
long,  open-ended  aluminum  cylinder  by  a  300  MHz  plane  wave  at  axial  incidence. 
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PROGRAM  FfTH  ( INPUT* OUTPUT* TARE  60  =  INPUT * TARF8»T APEfl ) _ 

r 

C  PliM  TASK6--  INTERFACING  THF  FL-TO  METHOD  WITH  THF  METHOD  OF 

C  MOMENTS 

C _ CASE  I-  STFADY  300  MHZ  TFM  IRRADIATION  OF  A  19.0  CM _ 

C  OIAMFTFH,’  68. s  CP  LONG*  OPEN-ENDED  ALUMINUM 

C  CYLINDER  _ _ _ _ _ 

C  AXIAL-INCIDENCE  plane  wave  with  COMPONENTS  F7  AND  HX 

C  MOOFLFD  USING  SCHELXUNOFF  EQUIVALENCE  THEOREM  WITH _ 

C  PROF.  DON  WILTON’S  M-O-m  APERTURE  CURRENT  DATA 

C _ ?_4_  X  163  X  ?4  CELL  CUBIC  SPACE  LATTICE  IS  USED _ 

C  UNIT  CFLl  riAMFTER  a  DX  =  0.5  CM  =  WAVELf N6TH/200 

C  EVEN  SYMMETRY  ABOUT  LATTICE  PLANES  X  =  ?4.5»DX  AND _ 

C  7  -  ?4.0*OX  IS  ASSUMED 

C  SOFT  LATTICF  TRUNCATIONS  ARt  USED  _ _ 

C  PROGRAM  IS  OPTIMIZED  FOP  THE  CDC  STAS-1 00 

C_ _ _ 

RFAL  UIR.  MUZ*  JXMAG,  JXPHA,  JXSRC,  JZMAG*  JZPHA,  JZSRC 

D I M  F  N  S I  ON  A (97 3504) .CEXR < 97  6  00) *CEYB( 97800) *CEZR (97800) . _ 

1  7 ( S936 )  »  A AA (600) »JXMA 6(600)  » JXPHA (6  00) *  JXSRC (600 ) * 

?.  JZMAO  (6  00J  .  J  Z  P  H  A  (600  )  *JZSRC(600)  *A  A  (gS)  *DO(?4)  » 

3  DF(?4) ,CA(9) ,CH(9) *EPS(3) *SIG(3) 

_ COMMON  A.CEXP*CFVh*CEZ8 _ 

DESCRIPTOR  01 . D? , 03 , D4 , D5 ,NFD , A  1 , A? , R V ,BW ,BX . RY , RQ 

BIT  ev.HPR (600) *  PW  * RRW (600) *BX*BrtX (600) *BY*PHY(?4) * _ 

1  RO  »BHO  (600 ) 

_  ASSIGN  RO . RPO  (1  16  0  0)  _ _ 

ASSIGN  PV  »RRR ( 1  1600 ) 

_ ASSIGN  PWfFPWtl >600) _ 

ASSIGN  HX  « PRX ( 11600) 

_ASS_IGN  BY  «PPY  ( 1  1  ?4  ) _  _ _ 

ASSIGN  n4.DD(ll?4) 

_ ASS  LGN  D5_»  OE.il «  ?4 ) _ 

T 1  a  SFCONO(CP) 

_ PPJLNI _ L5JQjl_L1 _ 1 _ 

ISO  FORMAT  (FP0.5) 

C  . I.  PROPLFM  PARAMETERS..... 

.  ..  _ FPJECL  _ _ _ _ _ _ _ _ _ _ 

DX  a  0.005 

_ MPP  =  3 _ 

DATA  EPS/  1.0,  1.0,  1.0  / 

_ DATA  S.IP/  0,0*  3. 7F  *7  *  0.01  / _ _ 

NM A  X  a  800 

c _ _ 

C  . II.  BASIC  AND  DERIVED  CONSTANTS..... 

_ PI  ~  3,H.IS92G5 _ 

MU7  =  4.0  *  PI  *  1.0E-7 

_ FJ>iZ_=_£^f<54F-_L2_ _ 

DT  a  nx  /  6.0E*P 

_ NHAIF  a  0.5  /  FREQ  /  DT _ 

R  a  DT  /  ?.0  /  EPSZ 


j 
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R  A  = 

nT«*?  /  DX**2  /  MU7  /  EPSZ 

RR  = 

DT  /  OX  /  MUZ 

PC  3 

1  .  0  /  RB _ 

RO  = 

?.0  *  PI  *  FRFR  *  DT 

1PUJL 

*  J) _ 

c 


Ctniq)  =  (L. 

CP(IIO)  =  0. 


RW  =  QPVMK  0  ( 1  «  25 *  BW ) 
RX  =  OPVMKZ (24.2SIBX) 


DO  ?  1*1, MPR 

EAF  *  P  *  510(1)  /  EPS ( I ) 

CA ( I )  *  (1.0-EAF)  /  ( 1 «  0*EAF ) 

C 

2 

CP(I)  =  RA  /  EPS ( I )  /  ( 1 • 0  *E  AF ) 

c 

. TIT.  1  DAD  UPrrnR 

A . 

C 

.....ZERO  INITIAL  FIELDS . 

Z ( 1 1593*)  =  0. 

A  ( 1  1 593* )  =  ft. 

C 

A(967569I5936)  =  0. 

c 

.....TYPE  OF  MEDIUM... 

RFAD(60.4.END=?01 ,ERR*?oi  ) 

(Z ( I 1 . T*1 ,600) 

♦ 

FORMAT  (75F1.0) 
READ(60,4,ENO=201,ERR=201> 

(Z(I)  .1*2601.3200) 

READ < 60.4, END*?01,ERR*201 ) 
DO  5  J*2  » 1 3 

(7(1) ,1*1301.1900) 

5 

JDEL  a  (J-I)  *  5936 

A  (JDFL*  1  15936)  =  Z(1I5936> 

DO  *  J*15?,163 

_ JDFL  =  (J-l)  »  593* _ 

C 

6 

A  (JDEL*1  *5936)  =  ZHI5936) 

RFAD(60.4,FND*?01 ,ERR=201 ) 

READ(60.4.END*201 .FRR=?ft] ) 

(2(1  ).  1*1.600) 

(Z(I) .1*2601,3200) 

READ(60,4»END*201,ERR=201 ) 
J3Q  -Z.  J=14  t.I50 _ 

(7(1) ,1*1301.1900) 

JDFL  =  (J-l)  *  593* 

_ 7  A  ( JDEL+1 15936)  =  ZAH332M 

C 


READ(60,4.END*?01«ERR*2ftl)  ( Z 1 1 ) , I  *  1 . 600 ) 


REA0(60»4,END*201 ,ERR  =  201 ) 
_ READ (60.4. ENDsPOl «  ERRsPft  1  ) 

<Z(I>, 1*2601,3200) 
(Z(I) ,1=1301,1900) 

J  =  151 

_ JDFJ _ *  (J-l)  »  5936 

c 

A  (JDEL*1 15936)  =  Z(l*5936) 

c 

. APFRTURF  FXOITATI ON . 

RFAD ( AO. 1 ft.FNDsPftl  .FQRsPftl ) 

( JXMAfi ( n .1*1 .6001 

RFAD(60*11»FND*201.ERR*?01) 
RFAD (60.1  ft. FND*?01 .FRRs?ftl ) 

(JXPHA(I) ,1*1,600) 
(JZMAfilll .1*1,600) 

RFAD (6ft. 1 1 ,FNO*201 ,ERR*201 ) 

(JZPHA(I)  ,1*1,6(10) 

53 


n  r> 


11  FORMAT  (25F3.0) 


C _ 

C  . MEDIA  coefficients 

_ _ ASSIGN!  Alt._DYN.t6QJ) _ 

ASSIGN  A?,.OYN.60Q 

_ ASSIGN  NF Qt .QM. 6_Q Q _ 

A 1  =  0. 

_ A2-  s.-Q-. _ 

no  ?3  J  =  ? « 1 63 

_ JDEL  =  ( J-U  »  SiE36 _ 

JOFM  =  (J-l)  *  600 

C _  _ _ _ -  .. 

NFO  a  A ( JOF L ♦ 1  *600) 

_ DO  20 _  J  J=1  _ 

BV  =  NFD.FO.JJ 

AJ  a  QBVCTPL (CA ( JJ) ,RVt A1 ) _ 

20  A?  =  OBVCTPL (CH ( JJ) ,RV»A?) 

_ A  ( JDF.L*  I  1 6 0  01  E  A 1 _ 

CF XR  (  JOEM* 1 1600)  a  A? 

C _ _ 

NFD  a  A(JDtL>1301 1600) 

00  ?1  JJ=1,MP«  _ 

BV  =  NFD.FO.JJ 

A1  =  QBVCTRUCA  (JJ)  .RV1A1  ) _ 

?1  A?  a  QRVCTPL(CB(jJ, ,BV!A?) 

_ A_1  JHEL-t  130116001  a  A.l 

CEYB < JOFM.l |600)  a  A? 

£ .  _  _ _ 

NFD  a  A <JDFL*?601 1600) 

_ 00_2?  JJ=1,MPP _ 

BV  a  NFD.FO.JJ 

_ A I  a  QRVCTPL(CA(JJ) ,BV1A1) 

2?  A?  a  QRVCTPL (CH ( JJ) ,BV1 A?) 

_ A ( JDEL* 260 1 1600)  a  A1 _ 

?3  CF?B  < jnpM*l 1600)  a  A2 

_  FPFF  _  _ 

C 

_ T?  a  SFCOND(CP) _ 

PRINT  1 50 «  T? 


. IV.  time-stepping  LOOP . 

DO  20  0  NaljNM/ix _ ___ _ _ _ _ _ _____ 

C 

C - - 1_. . COMPUTE  aperture  source  currents . 

ASSIGN  DU.DYN.600 

..ASSIGN  0?,.DYN,600  _  _ _ 

ASSIGN  D3..DYN.600 

£ _ _  _ 

0 1  a  FLOAT(N)  -  JXPHAUI600) 

_ BP  a  Dl.LT.0.0 _ 

01  a  OBVCTRL(O.O.PQIDl) 

_ D_2  a  01  *  JFLD. _ _ _ 
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do  I  '  o  cj  M  o  |  1  o  j  o  u  !  loo!  !  oo  oo 


. MAIN  EX  LOOPS...,. 

963*M  I  573  )  *  A/658fl*M» 


D?  =  A<11?13*M*573)-A<5?77*M»573>*A<10513*MI573) 


03  s  CFXR  <6?7*MM573)  *  D? 
AifiSfiB  «513U=_D1>_03 _ 


uaaa  111:1  aKfiin  mwunn 


IF  ( M  ,F0  •  1  ?*5936 )  A  < 656?*M I 60 0 > = JXSRC ( 1» 600 >  ♦A  ( 656?*M  1 60 0  ) 

_3ILJCI1NT1NJL1E _ _ _ _ _ _ 

FPFF 


,F Y  ITERATION, 


ASSIGN  D?».DYN.574 

_  DO  4G  itA=I-*3- _ 

M  s  JPFL  ♦  5936* ( MA—1 ) 


A(7P6?*M)  *  0.5  *  <  A  ( 8487 *M) ♦ A ( 84PS  +  M ) ) 


♦A (fl489*MI2?) ) 


A(P487*MJ?4)  =  A (78B7*M|?4) 


. MAIN  FY  LOOPS . 

D 1  =  A  (7?6?*m;574)  »  A  ( 7fl87»M  1574 ) _ _ 

0?  s  A ( 986? $574 ) — A (9837*M 1 574 ) ♦ A ( 11 21 574 ) 

_ 1 _ -A  ( 11213*MI574) 

03  =  CEYF*  <6?6*MK*574)  *  0? 

A  (7R«7*m;574)  =  Dl  ♦  03  _ 

.  « , • • F  Y  ENVELOPE  COMPUTATION . . 

04  s  ORVCMPRS ( A (7861*Mtf 00) »BX 104) 

_ 05  a  V APS (P4 I  05  ) _ 

BY  s  05.GT. A (B51?*m;24) 

40  A<P5l?*M|?4j  =  08VCTPL (D5»BYI A (851?*Mt?4)  ) _ 

FREE 


. EZ  ITERATION, 

ASSIGN  Plt.DYN .599 
ASSIGN  0? « ♦ OYN • 599 
ASSIGN  03*. OYN. 599 

00  50  *A=1,3 
m  =  JOEL  ♦  5936*<MA-1) 
s  JOEM  ♦  600*(MA-l) 


c 

C  .....MAIN  F7  LOOPS..... 

01  =  A<R537*Ml~599)  *  A  ( 9 1 6?  *M  1599)" 

D?  a  A  (10513*  *1599) -A  (  1  051?*M|599)  «_A  <3901*M  *S_99_)_ 

1  -  A  (  98  37  *M 1 5  99  ) 

_ P3  =  CE?S(601*mhI599)  »  02 _ 

A ( 9 16?*M I 59Q )  a  01  ♦  03 
C 

C  . J7  APERTURE  SOURCE  CONDITION . 

IF (M.FQ. 12*5936)  A ( 916?*M 1600 )  =  JZSRC ( 1 t 600 ) ♦ A ( 91 6?»M » 60 0 ) 
C 

£ . . E  7  ENVELOPE  COMPUTATION . . 

04  =  VARS (A (9737*M|?4) 104) 

_ HY  a  04 .GT.A  (9787*MI £4) _ _ _ 


A (9787  I ?4 )  =  Q8VCTRL (D4.HYIA (97B7*MI?4) ) 

04  =  Q8VCMPRS  <A(9161*MI600).6XID4) _ 

OS  a  VABS (04105) 


c 

BY  a  P5.bT.A(9612*M<Z4) _ 

SO  A ( 981 2*M | ?4 )  a  QPVCTRL (0S.BY1A (9fll?*M124) ) 

c 

. HX  ITERATION..... 

DO  60  MAa 1,3 

c_ 

M  a  JOEL  ♦  5936* (MA-1) 

c 

. MAIN  HX  LOOPS . 

01  =_A (39QL*M»599) 

0?  a  A(1951*MI599)*»A(1 9?6*M  1599  )  ♦  A  <  3?26*M  1599) 

1  -A (9162+MI599) 

c_ 

A ( 390 1  *M  1 599 )  a  D1  ♦  0? 

C 

. HX  ENVELOPE  COMPUTATION . 

04  a  VABS(A(4476*M124) IDA) _ 

BY  a  D4.GT.A <4526*M124) 

_ A  ( 4526. M I 24 )  a  Q8  VCTPL ( D4 . BY  1  A ( 4526 ♦M | 24 )  ) 

D4  a  QRVCMPPS <A <3900*MI600) .HXI04) 

05  g  VAflSfOAlPS) _ 

RY  a  05.GT.A (4551*M»24) 

GO  A (4551.MIP4 )  a  08 VCTRL (05  »HYIA(4551*M124)  ) 

C 

FRFE 

C 

.....HY  ITERATION,..,, 

ASSIGN  D1 . .DYN.59B 

C 

ASSIGN  DP..0YN.598 

DO  70  MAa 1.3 

M  a  JQFL  ♦  5936* (MA-1 ) 

C 

£_ 

.....SOFT  LATTICE  TRUNCATION . 

AA  ( 1 )  a  0.5  *  (A<1I137-»M)*A(11138*M)  ) 

_ AA(2J?2)  =  0.333  *  (A(11137*MIP?)*A(1113B4M|??» 

1  *A (11139*MI?2) ) 

_ A_A_{ 24 )  =  0,333  *  (AJ11159*M)  ♦  2. 0»A  ( 1 1 1  6Q»*)  ) 

AAA (11600)  a  «8VXPNQ(AA(1»24) . faW 1  A A A ( 1 1 60 0 ) ) 
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on  I  o  n 


C 

A  (10512  ♦  M«ft0  0>  =  QBVCTRL  (AAA(1<600)  ♦ BW t A ( 1 051 2*H < 600 )  ) 
A(1051?*MI600>  =  QHVCTRL(0.0«BX|A(10512*MI600)  ) 


. HY  ENVELOPE  COMPUTATION . 

_ QA  a  VAB.S  ( A  LLLflflXtflL 2 411.04  ) _ 

BY  a  04.GT.A ( 11 16?*M»24) 

70  A ( 1 1 1 6?*M I ?4 )  =  QBVCTRL (D4«BYtA (1 n6?*H<2») ) 
FRFF 


c 

. HZ  ITERATION..... 

ASSTfiN  D1..DYN.S73 

c _ 

ASSIGN  D2. .0YN.573 

DO  80  MA= 1 .3 

M  a  .inFI  ♦  5936*  <MA-1  1 

c 

c. 

_ -SOFT  1  ATT  T  CF  TRUNCATION . . . . 

AA ( 2)  a  0.5  *  (A(5877*M)*A(5B78*M) ) 

AA/3t?J)  a  0^.  333  *  ( A  ( 5fl77*M  S  21 )  ♦  A  ( 5ft7B^M  1 21 ) 

1  ♦A(5B79*M»21) ) 

AA ( 24 )  =,0.5  ♦ _ CA(5B9fl»Ml-  ♦  A-L589.9*itU _ 

AAA ( 1 <600 )  =  QflVXPND<AA<ll24) . BW * AA A ( 1 1 60 0 ) ) 


BO  A  <5?51*M|600)  =  Qfl VCTRL ( 0 ,0 f BX I  A ( 5251 ♦MlbOO ) ) 

—  FJ?E£ _ _ _ 

«?  continue 

IF  f  MfiAl  1  -I  F.  «i3»  fiO  TQ  <54  _ 


C 

c  ...TRANSVERSE  PLANE  NO.  164..... 

M  s  16?  *  5936 


. .  .  . .EX  »  EZ  TRUNCATIONS. 
M  656?»M  t  fe?5 )  =  A (5937 »M  <  625 ) 
A  (~5937*M|6?5)  =  A (7237*M»625) 
A (7?37»M<6?5)  =  A(6?6*H<625) 

A ( 9 1 6  ? ♦ M 1 6 ?  5 )  =  A (BS37*MI6?5) 
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rvr>|  |  I  i  r»  ir> 


A  (R537.»M*F?5)  =  A {9R37*Mj62b) 
A<9837*M|8?5)  =  A <32?8fM|6?5) 


. HX  ITERATION . 

ASSIGN  Dlt.nyN .5 9 9 _ _ _ 

ASSIGN  D?. .DYN.599 

_0I  _*_A1 350 1  ♦M  IS  991 _ .. _ _ _ _ 

DP  a  A ( 1951+MI599) -A ( 19?6*M*S99) *A (3??6*M»S99) 

A _ "A  ( Q1 62*M I  599  ) 

A ( 390 1  ♦M  | 599 )  s  01  ♦  02 

0&  a  VARS  (A  (4476*Mt?4)  104) _ _ _ _ _ 

BY  a  04.GT.A (45?6*M|24) 

A  t4526»Ml?4)  a  QflVCTBL  t 04 «BV I A (4526+M I ?4 )  )  _ 

D4  a  OAVCMPRSIA  (  3900  *M  I  f>0  0  >  «6XID4) 

D5  a  VAPS (04  I  PS ) _ 

BY  a  OS •  GT ,  A (4551 ♦PI24) 

Aj_455_L*_P L? ILL.  =  Q8VCTRL (05«BYI A (4551 *M|24)  ) 

FREE 


HZ  ITERATION.,... 

ASSIGN  0lt.0YN.573 _ 

ASSIGN  D?« .DYN.573 

AA (?)  a  0.5  ♦  (A(5877*M)*A(5678*M) ) _ 

AA(3I?1)  *  0.333  *  { A (5877*M»?1 ) *A (5878*M»21 ) 

A _ _ _ ♦A(5P79*HI?1  )  ) _ 

A A ( ?4 )  =  0.5  *  (A(SR9R*M)  ♦  A(5B99*M)) 

—  AAA 1 1  IftOO)  «  QP-VXRND LAAHA2±LshMiAAA ( LLfeMli _ 

A(5R76*M|?4)  -  Q8VCMPRS  (A  (5?5?*M»600)  »BW»A  (5«76-»Ml24)  ) 

01  =  A (5277 «M i 573) _ 

D?  a  A(65P8*MI573)-A(65?*M$573)*A(1951«-M|573) 

_ 1 _ _ -AUS52*Hl5IU _ 

A  (5277*M|573)  =  01  ♦  0? 

A  <5251 +MI600)  =  QP VCTRL ( A A A ( 1 1 60 0 ) t BW I  A (5?51 ♦  Mlf.QQ )  ) _ 

A  (5251+MI600)  =  QRVCTRL (O.OtbXI A (5?51*MI800> ) 

_  FREE _ 

94  T3  a  SECOND ( CP ) 

PPINT  150*  T3 _ 

C 

C . . FIELD  ENVELOPE  PRINTOUT  ROUTINE . . 

00  100  LcNHALftNMAX.NHALF 

_ 1F_LM»EQ»L)  QQ_TO_  1QJ _ 

100  CONTINIJF 

_ IF(N.EO.NMAX)GO  TO  .101 _ 

GO  TO  199 


C  . AT  H0RI70NTAL  SYMMETRY  PLANE . 

_ UL1 _ I£.t.N«EQ«NMAXl_I-P  UNrJ _ 

PPINT  10?»  N 

_ 10?  FORMAT  ( 1H1  «5?X«27HEZ  F-NVELOPE  FOR  TIME  STF.P  atlS. 

1  //»(S?X»15HPLANE  Z  =  ?4*DX»//t?X» 1HJ,/) 

_ CAI  I  FNV  (3851  «RC.  IP1JN) _ 


PRINT  104*  N  _  .  _ 

104  F0RMAT(1H1,5?X*27NHX  ENVELOPE  FOR  TIME  STEP  **I5* 

—  _  1 _  7/*6?X*15HPLANE  Z  =  ?4«Qx,//,2X*lHJ«/) _ 

CALL  FNV(4526«376,7,I PUN) 

£ _ 

PRINT  105,  N 

1(15  FORMAT (1H1*52X*27HHY  ENVELOPE  FOR  TIME  STEP  «,I5, 
1  //*62X,15MPLANE  Z  =  24*0X ,//*2X * 1HJ ,/ ) 

CALL  ENY  ( 5 22fe.il7A. TjIPJJNJ _ _ 

C 

£ _ ...,.AT  VERTICAL  SYMMETRY  PLANE . . 

PPINT  106,  N 

106  FORMAT  flhl *5?X * 27HF.Z  ENVELOPE  FOR  TIME  STEP  »«I5* 
1  // * 62X *  1 7HRL ANE  X  *  24.5*DX*//*2X* lHJ* /) 

_ _ CALLENV  l3676tRCjJPUNJ _ 

C 

_ PPINT  108*  N _ 

10fl  FORMATHHl  ,5?X,27HHX  ENVELOPE  FOR  TIME  STEP  **I5* 
1  //,6?X,17HPLANE  X  =  24 ,5»DX * //*2X * lHJ,/) 

CALL  ENV( 4551*376, 7* I PUN) 

C _ _  _ 

PRINT  1  OR ♦  N 

109  FORMAT ( 1H1 «5?X*27HEY  ENVELOPE  FOR  TIME  ST^P  «%15> 
1  //* 6?X » 1 7 HP LANE  X  «  ?4,5*0X»//*2X*1HJ*/) 

..  CALL  ENV(2576*RCUPJiRl _ 

C 

199  CQNT  lMlE _ _ _ 

200  CONTINUE 

?(U  T 4  =  SECOND  (CPI _ 

PPINT  150*  T  4 

STOP _ _ _ 

END 

SUBROUTINE  ENV (LOCA,SCALF,IPUNCH> 

DIMENSION  A(973504) ,IP{4100) *NN(25) 

_ COMMON  _A_  _ _ 

C  . ~  . . 

_ 00  8  1=1*25 _ 

«  NN ( I )  s  I 

c _ _ _ 

00  1  LY*2, If 3 

_  LOC  =  LOCA  ♦  (LV-1)»S936 _ _ 

LOCI  =  1  ♦  (LY-1)*?5 
_ IPILOCM25)  =  SCALE  »  AU0CI25) 

1  A(L0CI?5)  =  0. 

c _ _ _ 

IF ( I PUNCH ,FO . 0 ) GO  TO  4 

_ 00  2  L Ys? ,  1  63 _ _ _ _ 

LOCI  *  5  ♦  <LY-1)*25 
_ LOCI  I  =  LOCI  ♦  19 

2  *RTTE(8»3>  (IP(LL) *LLsLOCI*LOCII) 

3  FORMAT  (1 016) _ 
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_ 4  DO  5  LY  =  2»lf>3 

LYY  a  165  -  LY 

LOCI  =  5  ♦_  tLYY-lL»25 _ 

LOCH  a  LOCI  ♦  19 

_ 5l  PR1NX  i.Yy^,  (LP1LL 1  •LL=LOCI  .LOCIJJ 

6  F OPM AT (lXtI3«SX«?076) 

_ PRINT  J  ♦JMNlLU  f  LL=5  *  J?  4 ) _ 

7  FORMAT (//«HX.20I6) 

C _ 

RFTUPN 

_ END.  .  ....  .  _ 
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2.6  Problem  F  --  Task  5,  Case  2 
(Section  7.4.2  of  Volume  1) 

The  following  11  pages  list  the  computer  program  for  the  24  x  163  x  48 
cell  --  800  time  step  run  of  Problem  F.  The  problem  solved  is  a  hybrid 
method  of  moments/FD-TD  analysis  of  penetration  of  a  19.0  cm  diameter,  68.5 
cm  long,  open-ended  aluminum  cylinder  by  a  300  MHz  plane  wave  at  45°  incid¬ 
ence  relative  to  the  cylinder  axis  and  TM  polarized. 


62 


on  no 


PPOGR  AM  Full  (ii'iPuT,UuTPUT»l4r-'Fr(iiH,PUT»TAPtr-TAPri> 


<: 

f  rtff  TASNb--  I  N  IF *F  *»C  I  No  1  ri  -Tn  mFTFiOu  mT'-  T  nf  "KT<"Ui, 

C  ()T  MtiKhNlb 

C  CaSF  II —  S  T  r  a  i ;  Y  3t>  0  *HZ  Tf.-i'l-LAKlZtii  PLfli«t-»«vt-  i  ^k*n  I  aT  I  GN 

C  (F  "  1R.0  Cm  1'l^f.lTh*  6P.S  CM  LiMot  fiHrM-F.Mir.iy 

C  CONOUCT  Il.G  CYLInUF.p 

C  OHL  I  our  -  I  NC  I  OF'nCF  PLANE  WAVc.  kiln  CUMPUNF-W  T  S  fcZ  A  ,U  HA 

C  MuOhLtD 

C  USING  SCMF  LKUNOr  F  F  OU I  v  a(_  r  r-CF  ThHi«tM  run  kIlTI^/ 

C  GLISSUN  m-u-p  A  P  FI  P  T  i.l  •*  F  COr.PtNT  i  i  A  I  A 

C  ?4  X  163  X  4 F  CtLL  Cun  1 C  SPACt  L A T T 1  Lt  IS  UStO 

C  UNIT  CELL  DIAmf.TF.P  =  OX  =  i),r»o  Cm  =  wAVt  LtM-TH/POO 

C  F.  V  F  N  SYMMF.ThY  ahC'UT  LATTH-F.  pLaNF.  X  =  ?A.S*i/A  IS  ASSli^cU 

C  SOFT  LATTICF.  TRUNCATIONS  Ant  OScL 

C  PROGRAM  IS  OPTImIZmi  F  Or  Tut  Ll  (  STap-IOo 

C 

pfAL  PUR*MUZ»  Jaman,  JXPhA»uaSpC  .  J/P  A(-  ,  jZppa  *  oZSRL 
U  I  PENS  I  ON  A  (  ?S420OO  )  *  Z  (  lSSOU)  ♦  A  A  A  (  1?U0  )  •  HUF  (hfli)  )  .  afi  (4r>)  , 

1  00(64)  ,  UF  (4H)  ,0F  (40)  *(  A  (4)  ,CH  (H)  *F  r-S  (  3)  *SlG  (  J) 

01  PENS  I  ON  JXMAo  (  1  POO)  »  JXPPA  (  1EC  0}  ,  JXSRC  (  lP.Ui:  )  , 

1  j/m ah ( 1 60  0 )  »  JZpp A  ( 1 6 C 0 )  »JZSRC(16uo) 

COMMON  A 

OF  SCR  IP  TOP  D1  ♦OZ,»03,04,0S»0n»r  F  o,  A  1  ,  A6,HU  «f>V  »*W  «nX  ,  hy  ,p,Z 
UFSCPIPTOP  HO 

HIT  HU«  Hbli  ( 60<O  ,H\/fnnV(l?Oli)  *i  I1Z0U)  •oA«HftX  (1600)  » 

I  HY.HbY  (?4)  ,p2,-<-iZ  (4H) 

HIT  HO,  FFO(lZllO) 

ASSIGN  HOfHPUd  ?nOO) 

ASSIGN  HV  ,HH\/  ( 1  i  I?00> 

ASSIGN  b*»RH« < 1  *  1600) 

ASSIGN  bX  ,FhX  (  1  J  1  6U0) 

ASS  I G Hi  HY  ,i-r  Y  ( 1  5?4) 

ASSIGN  H  7  ,  HP  /  (  i  j4t() 

assign  f(mfwj(  i ;  i^oo> 

ASSIGN  04,f.()(  1  5?4) 

ASSIGN  OS, OF.  (1UH) 

ASSIGN  06, OF (1446) 

T I  =  SFCONO(CH) 

PRINT  ISO,  T 1 
ISO  FORMAT  (  F  ?(:  ,  S  ) 

..,,,1,  PPUPLFM  PAHAMtTciPS  .  ,  ,  ,  , 

F'RFO  =  3, CF  *b 
OX  =  0 . 0 0 H 
MPR  =  .3 

1'aTA  fps  /  1  , 0  ,  1.(1,  1,0/ 

uaT A  S Is  /0,0«  3,7c *7,  0.01/ 

NMAX  =  HOC 

. II.  HASIC  AnO  UFHiVFl;  constants . 

PI  =  3.141Sc«?r'S 
MUZ  =  4.0  *  Pi  4  1.0t-7 
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r>  r>  no 


r_DHZ  -  *  •  F  ~  i  'i 
i/l  =  Ox  /  f- .  0 1  ♦  ri 
I'-htLF  =  t'.b  /  FPtu  /  OT 
«  =  1 1 T  /  ?.()  /  t  HK/ 

=  i)T**r  /  HA****  /  'Jin?  /  LPSi 

=  I > T  /  !  X  /  I'M!/ 

«f  =  1  .  (I /I-  H 

r*  I  =  ?.('•  *  HI  *  F  H  F  O  o  LT 

i  r  l  •  N  -  ( 

C 

c,.  (  ]  ;m  =  (, . 

CM  UM  =  o. 

~t  V  =  ft  h  V  t'  F  0  (  1  »  K  *3  ?  n  '*  ) 

ox  =  ( ?*  Mb  ;r  x ) 

ur  ?  i  = 

r.  A  F  =  D  *  S I  iJ  (  i  )  /  f  *  S  ( I  ) 

LA  (  I  )  =  (  1 .0  -  F.  AF  )  /  (1.0  ♦  t  Al-  ) 

CM  I  )  =  HA  /  r.HS(l)  /  (l.o  +  cAF) 

C 

. HI.  LOAD  VtCTO*  A . 

. ztou  initial  fulus . 

ZddbSOO)  =  o. 

A  <  1  I  ISMiO  )  =  it  . 

MPbir’iSbOldb'iOO)  =  0. 

. TYHc  OK  Mti)IU**.  •  •  .  • 

*F  AD  (  00 .4  « t.  P>.r>  =  r'G  1  «fc  df<=?0  1 )  (  Z  (  1  )  *I=l*bOG) 

4  F  0  U  ft  t  (7HF  l  #(t) 

"t-<>u(O0«4«r:M;  =  e0  1.FKP  =  2'fil>  ( Z  (  1  )  »I  =  St)01«bbQu) 

*  t  A  D  (  Ml  .4  «  t.lv0=Z  0  1  «  KPh=?(i  1  )  (  Z  ( 1  )  »I=Kb01»31U0) 

DC  b  J  =  £  .  13 
jr.Fl  =  (o-l)  *  1  r»  S  0  0 

b  A  (JOt  L*  H  l'xbOO)  =  Z(ldbbOu) 
o(.  t  J=lb?»  lb-1 
JOFL  =  ( o  —  1 )  *  lbbOO 
o  a  (oOKI  *1  ;  loSOO)  =  ZddbbOM 

HtAD(Afl.4»fcAo-<?01»r. =  0  1 )  ( Z  (l)»I=ltbOO) 

Hf  AO  (AO  »4  ,tM)  =  ^0  l.tPP  =  <?Oi)  ( Z  (1  )  »  I=b001*bb0t) 

ri  F  <•  u  (  A  0 . 4  *  F  4 1  >  -  C  (j  1  »FWK  =  ZOl  )  (Z  (I  )  »  I=Zb01»31UO) 
oO  7  , J  =  1  4  •  1  Fu i 
JnrL  =  (o-l)  *  lbbin) 

7  4  ( jl)F  .(.♦  1  5  IKSf/fi  )  =  2  (  1  •  1  bU'  U  0  ) 

HFAO(nO**«tM>  =  KOl«fc.KPsi*(il)  (Z  (1  )  *  I  =  1 »  bO  U ) 

*-  c  A i '  (  A.)  ,4  .!-  M‘  =  r.O  1  »F  *P=2fi  1 )  (  Z  ( 1 )  »  I=b(>01  *bbl)<; ) 

*F  AO  (  ao  .4  ,t  f.O  =  KU  1  .r.ww  =  ?01  )  (Z  (  1  )  .  I  =  ?bOl»31uu) 

J  =  Ml 

•  IDF  L  =  (J-l)  *  lbSOO 

**  (oOF  l  ♦  1  ;  lSbOl) )  =  7(lilSAU0) 

C  . aFD'IuaF  KxCITaTIun..... 

r-  r  A!)  (4(1  .  1  *- O'M  =  ?0  1  .hKD  =  ?{i  1  )  (OAPAO(I)  .1  —  ltlKOO) 
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PF  6  0  •  1  0  *F  N\-  =  C  0  1  tf  PP  =  2U  1  )  (OAHhti  (  I  )  «  I  =  1  t  1  2  f; ) 

r<  F  A  H  (  M)  «  1  0  .  F.  f  I  u  -  2  0  1  »  c!  P  u  =  2  0  1  )  (  J  /  'i**  (i  (  1  )  «l  =  l»lfi»0) 

*F  AO  (60  •  lo»t  NG  =  2G  1  .t-(«  =  ?Ol  )  {  wi  HnA  (  I  )  *  1  =  1  t  IrVi)) 

10  RiMHAT(l^^.?»/.nr^.?) 

1A  FORMAT <?0F 3.0) 

C  .  .  *  •  .  S  YN  "if  T  H  Y  AmOuT  /  =  ''A>n*oA...«. 

GO  12  J=?*l<~3 

JGr  L  =  ( J-  1  )  *  1  oo<j  ii 

f\yA  a  -  a? 

i  >  0  1?  J  A  -  1  «  3 

jo  =  JGFL  ♦  (ua-1)*2oO<>  ♦  1 

-h;f  <  l ;  auo )  =  «(Jiu.no) 

IF  (  JA.FO.  J)  Kma/=a,< 
a  ( jH  +  efio; 20)  =  -iuf(o7o;?o) 

1*0  1?  K  =  l«>4 
f\  A  =  (  K  -  1  )  -»  2  0 

Kl-  =  (KMAX-K)  *  20 

1?  A  (  Jf.»+KK»?S)  =  -i  *JF  (2a*  1  ;?o  ) 

c 

c 

C  . Of  0  I  A  COF  F'F  ICItMS . 

hFSII'F  A  1  ♦  .  uY1'  .  120  0 
AO?  I  Go  A2. .GYN. 1200 
ASSIGN'  N  F  i »  »  »  0  Y  *  i «  1  2  (Ml 
A 1  =  n. 

A?  =  0. 

GO  21  J  =  2  »  1  G  3 

JGF  L  =  ( J-l  )  *  loooO 

DO  21  JA=1,3 

JH  =  ( JA-1 ) *2000  ♦  1 

JC  =  JR  ♦  1 2S  0 

NFf*  =  A  (  JGF  L*  Jo  S  1200) 

U(i  20  JJ=1*mPk 
r  V  =  nFIj.FG.  JJ 
A  1  =  GMVCTHL  (Ca.  (JJ)  »oV5Al  ) 

?0  A  2  =  OHVCTpL  (C->  (  JJ>  .pV;A2> 

A ( JDKL* JR s 1200 )  =  A] 

21  A { JOF  L  ♦  JC  •  1200)  =  62 
KOF  t 
C 

T  2  =  SF" COfio  (  C.a  ) 

PRINT  1  b  0  »  T2 
C 

C  . IV.  T I Mh-STtPP 1 NG  U0(  H . 

I/O  200  N=1»GMAA 
ASSIGF  D  1  .  ,('Y  ii  1  2  0  0 
ASSIGN  U?»  .DY'i.  1  2;mi 
ASSIGN  G 3* .GYM.  1  20o 
C 

oi  =  fluaT(n)  -  jxPha ( l ; 1200 > 

of)  =  0 1  .1  T  .  0  .  H 

i*l  =  GAVCToL  (  o  .  0  »  ru  5  ')  1  ) 
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o?  =  ni  *  «ij 

i.  ?  =  v  s  I  n  ( 1 1  ? ;  i.  o ) 

J>bMC  (  1  ;  lr(MJ)  =  jXl'fi(i(l!ROU)  *  Crt(l)  *  1*4 

c 

i»l  =  FLOAT  (  g  )  -  JZPha  ( 1 1 1*00) 

=  lil.LT.n.,) 

i *  1  =  ( ■  V C T " L  (ii.O.GOtt)l) 

op  =  *  l  * 

u?  =  vsiMiy;'*j) 

J/SPL ( 1 5 IfuI )  -  J7*- A<,  ( 1 ;  1*00)  *  Crf(l)  *  D  4 

F  F  1 1 

c 

m  C  A  L  L  =  7  ♦  IF  I*  (F’LOaT  ( f\i  )/*->•  <j  ) 

IM^CAI  L.t>T  ,s*>  MCALL=r>4 
C 

C  . TrtAOSvtHSt  PLANfc  'MU  •  1..... 

C  . t.X,  1 2  TPUnCaT  IONS . 

A (7S01 5 12N0)  =  A(l*l?r»«) 

A  u  ;  1  l* )  =  A  (  'i.  30  0  1  5  1  2S  0  ) 

A  (  1  0?01  5  12S0)  =  A  (  1Xt>1  5  12M* ) 

A  (  1?<=>1  5  12S0)  =  A  (?S701  5  l?bu) 

r 

C  . tY  iIFwaTIOn..... 

ASSIGN  L  I ,  . L Y f  .1  17a 
ASSIGN  '*?,  ,l;YO.  1  1  74 
ASSIGN  0  4,  .GY*'  .  11  74 
L 

A(H8?M  =  0  .  b  *  <a<1007G)*A(1uu77>) 

A(Fh?7  ???)  =  0.333  *  (A  (10(*7c>522)*A  (10077522)  *A  (10078522)) 

A(884A)  =  0.344  *  (  A  (  1  i)  OOf> )  *2 • 0*A  (1  00S*S* )  ) 

«  ( 1  0(*7s  5  f4  )  =  A(orsl524) 

c 

A(ion?G)  =  o.s  *  (  a  ( loion  ♦  a  ( loio?)  > 

a  (10027522)  =  0.333  A  (  A  (  10  101522)  *A  (  101025  22)  *A  ( 10103522)  ) 
a(1004u)  =  0.443  *  (A(10123)*2.0*A(10124)  ) 

4(10101524)  =  *  <  100  OU  24) 


c 

Ol  =  C A  (3 )  *  A ( oFSl ? 1  174) 

OP  =  A ( 1 1SS1 ; 11 74) -A u is?h5 i 174 ) ♦ A (14201  5  11 74)-A (  142021 1174) 
u?  =  0;<3)  *  Of 
a (AHSl 51174)  -  01  ♦  03 

fall 

c 

r  . IaaiSv/FkSF  planf.s  p  -  1*3 . 

on  f?  jy=i,mc«ll 

JI  FL  =  (UY-1)  *  4OS00 
C  .....ex  i T tKhT I  ON  .  .  .  .  . 


ASSIGN  01..0YU, 1173 
ASSIGN  I»2».i)Y«.1174 
ASSIGN  1*3.  .OYO.  1  174 
C 

On  30  NA=1,3 
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non  no  noon  on  n  no 


A  (?4?S2«-»- ?23)  =  a  (23i'>27*m;23) 

A(?4?0?+M)  =  0.5  *  (A  (242  7  7  ♦v')  ♦A(^4f’7M*C')  ) 

M?4203»m*21)  =  0.333  *  ( A  ( 242  7  7*f«t21)*A(  2*2  7  » 2 1 ) 

1  ♦  A  (2427W«'M;21  )  ) 

A(i>42?4*M)  =  0  .  t>  *  (  A  <?42*rt  ♦  •*  )  ♦  A  (  A42SS  )  ) 

a  (  P4277  +  M  ;  ?3)  =  A  (?4l77*P523) 

. MAIN  t  X  LOOPS . 

01  =  A( 1  S  S  2  7  ♦  ^  5 1173)  *  A (23027*m; H 73) 

0?  =  A  (  2W02*w»  1 1  73) -A  (14202*F  ♦  1 173)  ♦  4  <2h3b2+Mt  1 1  73) 

1  -A (2b3  7  7*m; 1  ]  7  <) 

u 3  =  A  ( lb777*">«  1 1  73)  «  02 
A (?30?7*P5 1 173)  =  01  ♦  03 

. JX  APFHTUPF  SOllHCt.  CONDITION . 

If  <M.Ln.i2*is^oO)  a  (2iooifi-i;i2oo)  =  jxspc  ( i ;  i2oo)  ♦ 

1  a  ( 23o o  1  ♦(« ;  l  ?o o  ) 

. FX  cNvtLOPt  COMPUTATION . 

04  =  VAbS  (  A  <23601  ♦M*24)  JI/4) 

t-  Y  =  (14 .07  .  A  (  2430  1  *t»t  24  ) 

30  A  (  ?430  1  ♦M  J  24  )  =  QnvCTkL  (04  ,  r>r  ;  H  (2430  1  +m  *  2*+ )  ) 

FPf  t 

. tv  IT tK AT  ION . 

ASSIGN  01 . .OYN. 1 174 
ASSIGN  02, .OYN. 1174 
ASSIGN  03».UY'M.1174 

lift  40  M  A  —  1 .3 
M  =  JOtL  ♦  lbbOO* {fc-A-1 ) 


SOFT  LATTICF  TPUNCA T I ONS . 

A ( 24326 )  =  O.b  *  (A(2Ss7F  +  m) ♦A(2SS77*m) ) 

A  (  24327 ♦M  5  2?)  =  0.333  *  ( A ( 2nS 7p.m{ 22 ) ♦A ( 2b^7 7*P 5 22 ) 

1  ♦ A ( 2r>t>7M.p; 22) ) 

m(2434<3»P)  =  0.333  «  <  A  (  2b3'*«"M-  )  ♦?.  0*  A  (  2S3'»4*M)  ) 

A ( 2SS  76 *F  *  24 )  =  A ( 243S1 ♦M 1 24 ) 

C 

A(?S5?6*N)  =  O.S  *  ( A ( 2Sn  0 1 ♦ M ) ♦  A ( ?S  NO  2*  M  )  ) 

A  (  2SS27  ♦  P  J  2? )  =  0.333  *  (  A  (  2->"0  i  ♦  M  J  22  )  ♦  A  (  2b<-  0  2* '  S  2  2  ) 

1  .A  (2bt-o3*F  •  22  )  ) 

A(2S44S*M)  =  0.333  *  (A  <2bb23'*f  )  ♦?.0*A(2bo24*M)  ) 

A (2SG01 ♦PS?4)  =  A (2SS01 *Mi2A) 

r 


i 
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on  r.  o  o 


C  . .  t-Y  LOOMS..... 

Ul  =  A  ( 1 17«4 )  *  A  1  174) 

■j?  -  A  (  ? 7 r>s  l  ;  1 1 7** )  —  A  ( *  1 1  7a  )  ♦  A  ( ?•*  / 0 1  *m;  1 1  74 ) 

1  -A  (2R702*mU  174) 

03  =  A  (  1V27N.'-'  i  1  1  74  )  *  02 

A  (?43sl  *m;  1  174)  =  i>1  ♦  f  3 

C 

C  . r  Y  r.NVt  LOOK  COMMUTATION . 

JL  =  VAt’S(*--(?442N  +  M$?4)?f,'4) 

r> Y  =  l>4  .OT  .  A  (  i'Shi'fltK  i  24  ) 

a (2Gn2n*m;24)  =  w«VCTML  (04,oY;a  <2S6264M*24)  ) 

c 

os  =  ohvcmk-ms  (  a  (?43?s*m;  i20o)  ♦t<x;i>s) 

Oh  =  VAnSfl'SSoo) 

N  ?  =  UN  .07  .  A  (  2">MS  1  ♦  **  »  4tl  ) 

40  A  (  1  ♦  ♦  4  h  )  =  Oh  V  CTRL  (  ON  »  o/.  *  4  (  2S0S1  4M  ?  4n  )  ) 

THE 

. t.7  ITERATION..... 

nSS  IG'j  0  1.  .t-Y  ..1144 
ASSIGf  02.  .OYN.l  WR 
ASSIGN  U„i.  .UYn.  l  1V4 

00  SO  MA  =  1  «  3 
v  =  jOFL  *  1 n  s  0  0  * (nA-1 ) 

. NAlf'  t2  LOOMS..... 

l’  1  =  A  (  ?()S()  1  ♦<*,  J  1  1  vi>4  )  *  A  ( ?S7o  1  *m;  1  1  44) 

O?  =  A  (?h3S/?^E  J  llvM)-A  (?A3M*Mi  11MV)  4A(  11S204M*  11^4) 

1  -A(2702t>*M;ilOM) 

1-3  =  A  (?17nl*M«  llvw)  *  0? 

A  (?S701*«; 1 1M4)  =  oi  ♦  03 

C 

C  . Jt  AMtuTuMf  S 0 U M L t  CONDITION. .... 

IF  ( v  .  M  0  .  l?*lSt>00)  A  «?S701*M}  1200)  =  JZSRC ( 1 ! 1200 )  ♦ 

1  A  (2S7014MU200) 

C 

0  . t  7  t  .VtLOME  COMMUTATION . 

A  (?hOe>h,*h‘l?h)  =  O.S  *  (  A  (2627N.Mi<>4)  .A  (2o30l4M»?4)  ) 
04  =  VAhS  (  A  (?N426*m;?4)  *C,4) 
hY  =  04. UT .A (20VS1 *N|24) 

a (2N4S1 *m;24)  =  OEVCTRL (04,SY}A (26VS1 4NI?4) ) 

C 

.)4  a  OhVCnMks  (o  (?S700mm;  12(»U)  .rXJDs) 
t-f  =  V  AMS  ( i-S  i  oo  ) 

«7  =  ON.GT .A ( EOV764MJ4W) 

So  A  (2NV7N.m;4m)  =  OsVCT»L  (l-b,B/ »A  (2FR764M;4o)  ) 

C 

C  .....MX  ITERATION..... 

uO  NO  ma=  1.3 
m  =  JOFL  ♦  1nS(<0«  (*A-1  ) 

C 
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no  on  o  on  onooo  n  on 


C  . MA  Il‘  HX  LOOMS . 

'  01  =  A  ( 1 1 5?H ; ] 1 44 ) 

,  n?  =  A  (  AA51  *v  ;  1  lvs) -A  { Kh?r>  +  *{  ]  144  )  ♦  A  ( 1020  1  ♦.*>  S  11 

1  -a  { 2t>7u  1  *f  *  1 1 *-** ) 

a  ( 1  1S?h*h;  1144)  =  i,i  ♦  ip 

. HX  tNVcLOPb  COMMUTATION.,... 

A  {  l?7b  1  ♦M  5  ?4  )  =  0.5  *  (  A  {  littl  i)  1  +  P  t  0*4  )*  (s  { 1  f  l  f  )  ) 

04  =  VAHS(A(12751  *m»24)  ;(>4) 

5Y  s  04 . G  T  .  A  (  1  2  7  75  ♦  '4 ;  P4  ) 

A  (1277H*m;?4)  =  UMVCTKL  (04,-»Y  «a  (1?775*m;c4)  ) 

OS  =  QpVOHPS  <  A  (  1  ]525*M  (UN)  »rx  J  OH) 

OH  =  VAt'S  ( 05  J  uo  ) 

H 2  S  D6.GT.A ( 1PB01HM14M) 

HO  A(l?ft01+M«4H)  =  •jSVC7ML(r)ti.’^Z*A(l?^01*>"*4r)) 

FREE 

. h  Y  I  TEH4TIUN.  .  .  .  . 

ASSIGN  Dl  ♦  .0Y!M.  1  lHfc 
ASSIGN  0?«.()YN.114h 

0(i  70  MA=1»3 
M  =  Jflfc'L  ♦  1550(i»  <«A-1  ) 

. SOFT  LATTICE  TRUNCATION . 

AA(1)  =  C  •  5  *  <  A  <24601  «>)  *A  (2Vt>02*M)) 

AM2I4H1  =  0.3j3  *  (  A  <?4m0  1  +■*  i  AH  )  ♦  a  ( 24602+N  ii>h  ) 

1  +  A  (P4H03*  4J4fe)  ) 

AA  (  48  )  =  0.6  *  <  A  (24H47*.*)  *A  (/•Rt>4h*M)  ) 

AAA  (111  20  0  )  =  wbVXMNl)  (  AA  (  I  ?  Ac)  ;  AAA  (  1  i  l^uL  )  ) 

A  ( 2S60  1  ♦m  5  4  A  )  =  QhVCMMHS  (  A  <26352*m;  I?00)  .CWiA  1  ♦M54a)  ) 

. MAIN  MY  LOOPS . 

ul  =  A (2M352+M ?1 1 HA ) 

Op  s  A  (2570?*1"  ♦  1 1  ) -A  ( 2570  1 +m«  1 1  HA  )  ♦  A  (  23002  +  m  { 1 1  H*' ) 

1  -a<23u«?7*m;11hm) 

a  (?83S?*M;  1  ISM)  =  01  ♦  u? 

a  (?83S1  ♦MU?0(>)  =  onYCTKL  ( AAA  u  ;  1?00)  A  <P.i3si  ♦ -i  1200)  ) 
A(?H351*M»lp00)  =  OAVCTHL  (  0 . 0  »t  X  ?  M28361  ♦  M»  1  200)  ) 

. HY  hNVtLOPt  COMMUTATION . 

A <?957H*M }24)  =  0.->  *  ( A ( 2ah2o*m; ?4) ♦ ft ( prtVDi *MS ) ) 

04  =  VANS  (A  (PHS7H  +  M24)  U'4) 

AY  =  fl4.GT  .  A  (  2^65  1  *M;  p4  ) 

70  A  (?4Hbl+MI?4)  =  OHVTTPL  (04.  <Y  5M2SH5)  ♦!"<£*)  ) 

F  PFE 

. HZ  ITERATION..... 

ASSIGN  01  ,  .OYN. I  173 
ASS  IGN  02, •  0  Y4  • 1 1 73 
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DO  8  0  *»A=1»3 

(•;  =  jdfl  ♦  is^uo* (f a-i > 

c 

C  . SOFT  LATTICfc  TRUNCATION..... 

A  A<2)  =  0.S  *  (  A  (  1S427*M)  *A  ( lb42G*M)  ) 

AA(3}  4b)  =  0.3J3  «  (A  ( lb427*«M4S)  +  a  (  lb42«*M.4b> 

1  »A  (  lb424»MJ4b)  ) 

A/i(48)  =  O.b  *  (  A  ( 1  b472*M  )  «•  A  (  lt>473*M)  ) 

AfAfDlifcli)  =  U8VXPN0  (  AA  (  1  «4rt)  ibN  i  AAA  (  1  *  1*00  >  ) 

A  (  1S4?6*M?48)  =  U8VCKPPS  (  A  (  14177»M;  1200)  *h«»  SA(  1‘j42«5*MJ4B)  ) 

c 

C  . MIN  nZ  LOOPS..... 

i » 1  =  A  (  1 4  2  0  2  ♦  »  1 1  7  3  ) 

0?  =  A(23027*M;il73)-A<7b2:7*M1173)  ♦  A  <  beS  1  ♦*!  1 1  73 ) 

1  -A (H8b2+Mi 1173) 

A ( I4?0?*p;  1  173)  =  01  ♦  n? 

C 

A  (  1  41  7b*“  i  1  200  )  =  UAVCTPL  (AAA  (  1  ;  1200)  »ota?A  (  l4l7b  +  PU200)  ) 

A  U4178  +  8J  1200)  =  (vbVCTPL  (O.Ovt-x;  A  { I4l7t>«-I*«  1*00)  ) 


C 

C  . HZ  KNVtLOPe  COMPUTATION 

04  =  VAbS  (A  (  1477b*M!24)  ;04) 
r.Y  =  04  .GT  .  A  (  l-D47b>M!  24) 


ho  A(lS47b*M524)  =  Ob VCT PL < 04 . n Y ♦ A ( 1 b4 It ♦« 1 24 ) ) 

fpf  e 

h2  CONTINUE. 


IF (MCALL.Lh.b3) bO  TO  04 
C 

C  . TPANbVEbSe  PLAM:  nu.  Ib4 

m  =  18?  *  lbbOu 
C 

C  . to  tZ  TRUNCATIONS.  .  .  .  . 


A (23001 *m; 12b0)  =  A ( ISbo 1 ♦M  » I2b0 ) 

A ( 1 S50 1 ♦M  *  1 2b0 )  =  A ( lb7bl*MJ icbO) 
a  ( l(S7Sl*MU?bO)  =  A <7S01*M; lcbO) 

C 

A <?S701 ♦M? 1?SU)  =  A ( 1 rt o  0 1 i 1 2b  0 ) 

A(1h001*m?12S0)  =  A ( lR2bl*M» l^bO) 

A  ( 1 1 *m  *  1 ?5u )  =  A  <  10?01*MI  lirbO) 

C 

C  .....bX  ITlRaTIOM . 

a SSI  ON  01 , .OYN. 1 144 
ASSIGN  u2* .DYN. 1144 
01  =  A ( 1 lu?b+M; 1194) 

U?  =  A  <88bl*b?  1  149) -A  <MH2t>*M«l  199)  *A  (  1020l*M  1 1  49  > 

1  -A ( 2a70 1 ♦M  *  1 199) 

h  (  1  lb?G*M;  1144)  =  U1  ♦  02 

C 

a  (  1P7S1  ♦*524)  =  (,.b  *  (  A  <  12I01*m;?4)  *A{  1212b  +  M;24)  ) 

1)4  =  V  AbS  (A(l?7bi+Mt?4)  L04) 
rY  =  D  4 . G  T .A ( 127  /  b ♦ m ; ?  4 ) 

A ( 1 2778+85 ?4)  a  ubVCTWL (04»rtY;A (127764MJ24) ) 

C 

70 


n  n  n  n  o  n  no 


05  =  OPVCmRRS  (A  (1  1S?5*m;  1^00)  ,<-X«05) 

Ob  =  VAHS(os;ub) 

b7  =  06  »  GT  .A  (  1  ?b0  1  ;  4tt  ) 

A<12R01*M?4b)  =  (JSVCTRL  <1)6,  1  a  ( 12«0  J  ) 

FRFt 


. hZ  1  T  K  W  A  T ION . 

AS$1C-N  i>l,.0YN.1173 
ASSIGN  02  ,  • UYN .1173 

AA(?)  =  0.5  *  (  A  <  154?7*M)  +a  (  ) 

AA(3»45)  =  (1.333  *  (A < 15427+m ,45) ,A < 1542h*M 5«5) 

1  ♦  A(1542'*«,M!4h)) 

A  A  (  4b  )  =  0.5  *  ( A < 1 64 72*m ) ♦ A ( 15473,n ) ) 

AAA(i;i?00)  =  UMVXPNO(AA  (ll'fb)  .Mu'SaAA  (i;i200)  ) 

A  <  15426+M  *46)  =  U8VCMPRS  (  A  <  1*1  77**:  121)0)  ♦»»  »  M  1S42**k  }4t>)  ) 

01  =  A ( 14202*Mi 1173) 

0?  =  A  (  23027*iU  1 1  73  ) -A  {  7S27  ♦(»  *  1  1  73)  *A  (HbSl  ♦  •  ,  1  i  73  ) 

1  -A(bMb2*i'«1173) 

A  (  1420?*MS  1  173)  =  1<1  ♦  0? 

A  (  141764.*;  1200)  =  G»<4VCT*L<AAA<m?00)*t»W**(l4l7h*Mil?00)) 

A  (1417b*N?l?00)  =  OHVCTRL  (0.0  *lrX  »A  ( 1  4l  7b*M«  1200  )  ) 

04  =  VAt-S  (A  (  lA776*b524)  *u4) 

b  Y  =  D4.GT . A < 15476+M5 24) 

A  ( 15476**  (24)  =  uHVCTbL<04,oY,MlS47fc*M»24)  ) 

FREE 

94  T3  =  SECOND <C2) 

PRINT  150,  T3 

. F'lELO  tNVELOPt  PbIMOuT  POUTINt . 

DO  100  L  =  l<bAl_F  ,Nmax  .nhalf 
IF (N.EO.L)GO  TO  101 

100  CONTINUE 
IF (N.FQ.NMAX)GO  TO  101 
GO  TO  199 

. AT  HORIZONTAL  ObSbK VAT  1  ON  PlaNF . 

101  IF (N.EQ.NMAX) IPUN=1 
I PUN  =  0 
PRINT  102,  N 

10?  FORMAT  (  lHl,5?X,?7Ht  X  EnVFLGPc.  FOR  TI*E  SI  tP  =  .I5. 

1  //,b?A, 15HRLANE  7  =  ?4*UX .//,2X, lnj,/ ) 

CALL  FNV  (bF.Ol  ,20,«C,  IPUN) 

C 

PRINT  103,  N 

103  format  ( ini  ,s?x,27mey  envelope.  for  time  step  =,is. 

1  //«b?x,  IsRRLaNE  7  =  24<»Ux,//,2X,  lFlu,/) 

CALL  FNV(1012b,?0,PC,lPUF ) 

C 

PRINT  104,  N 

104  FOkMAT  ( 1 H 1  ,5?x,27hEZ  ENVELOPE  F(.w  TIME  STEP  =,Ib. 

1  //,b?X, 15HRLANE  Z  =  24*dx,//,2a, lno,/) 
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CM.l  FNV  (1  1aS1,P»»»kC,  IPuN) 


OkINT  10b»  N 

10b  E  OPmAT  (  lh  1  ,b?A  *  27mhx  fNVELOht  E  Ok  TlMt  bTtP  =  ,IS. 
1  // •  b?  a  *  1  bt-ikL 2  =  t4*i>X  » // ,  PA  »  1  fl O  »  /  ) 

CALL  FnV  <  lP77o,PO,37b.7,lPUN) 

►'PINT  10b»  N 

1  0  4-  f  (ikMAT  (  lh  1  ,S?X,?7hhY  ENVELOPE  E  Ok  T  I  Kit  STcP  =  ,Ibi 
1  /  /  »  E>?X  »  lbftPLANfc  2  =  ,  //  »?X»  1  rt  j  »  /  ) 

CALL  fcNV(l*»lbl«20,37b.7,lPUN> 

10  7*  N 

1 1>7  M)(4»*AT  (1H1  ,c?A«27HhZ  ENVELOPt  K)K  TlMt  STcP  =  ,Ib 
1  // , hPX  »  lbhPLMNt  7  =  P44UX»//»2X»lriU»/) 

CALL  ENV  ( lb47b»?0»37fi.7*lPUN) 

. mT  VtRTICAL  SYKMtTKY  PL  ANt  »  .  .  ,  • 

HUNT  10n,  N 

ink  EOPmaT  (  lnl  ,b?A,P7KtZ  f-.MVf.LOkt  h  Ok  Tl^t  bTtP  =  ,Ib 
1  //,b2A,  17htPLANt  X  =  P4.b*OX»//,2X»lt~J»/) 

CALL  EOV ( 1 l47o,40,kC« 1PON) 

PklNT  109,  N 

lfi<4  fOP^AT  (  lhl  ,S?A*P7mhX  fNVtLOht  E  OK  TlMt  SIE.K  =  ,IS 
1  //*h?A «  17hPl.ANt  X  =  P4.b*DX»//»?X*ihJ»/) 

CAl  L  E  N  V  ( 1  ?P0  1  ,40,37b»7,I PUN ) 

PPINT  110,  N 

110  EOPMAT ( 1H1 ,S?X,^7Hfc Y  tNVELUkt  E  OP  TlMt  STtP  =,I5 
1  // «  6? X  »  17hkL  ANEI  X  =  £'4.b*OX,//,?X,lhj,/) 

call  FOV(l0lM,40,4C*IP0M 

i9«  continue: 

P 00  continue 

P01  T  4  =  SFCOru(Ck) 

PPTNT  lbO,  T4 
bT  OP 
tr  P 

bOKKOOTlNfc  f MV (LOCA,MJM,SOALt» IkUnCH) 

O  I*F.  NS  I  ON  A  (?b4^000)  ,  IP  (MOdO)  ,NMbO) 

COvmQN  A 

00  1  1=1, bO 
1  Mill  =  1 

IPdJMOOO)  =  o 

no  ?  L  Y  =  ? , 1 h3 

Lf C  =  LOLA  ♦  (LY-l)*lbSOO  ♦  4 
LOCI  =  1  ♦  (LY-11440 

IP(LOCUNOm)  =  SCALt  *  A(LOC»NOM) 

?  A  ( L  0C  S  E  UP  )  =  0  . 


c 

IF  <  IPUNCh  .(■<,.  U  )  C-0  TO  4 
00  3  LY=?»l63 
LOCI  si*  <LY-n*40 
L0CI1  =  LOCI  -  1  ♦  NUP 

3  wPITE(fl.4>  (IP(LL)  «LL=LOCl»LuClI) 

4  FORMAT (1016) 

C 

<=■  lxp  =  num/po 
00  0  LX=1»LXM 
LX  A  =  1  ♦  <LX-l)*?G 

LA  =  LX A  ♦  4 
L7  =  LA  ♦  19 
00  b  LY=?tlr3 
LYY  =  164  -  LY 

LOCI  =  LX  m  ♦  (LYY-1 ) *40 
LOCII  =  LOCI  ♦  19 

b  POINT  7  »  LYY.  (IP(LL)  «LL*LOCi. LUCII) 
7  FOPMAT  (1X*I3*3X.<?0I6) 
n  POINT  9.  (f  f4(LL)  «LL=LA,L2) 

9  FORMAT  //////> 

C 

RETURN 

EMJ 
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2.7  Problem  G  --  Task  5,  Case  3 
(Section  7.4.3  of  Volume  1) 

The  following  13  pages  list  the  computer  program  for  the  24  x  100  x  48 
cell  --  1800  time  step  run  of  Problem  G.  The  problem  solved  is  a  hybrid 
method  of  moments/FD-TD  analysis  of  penetration  of  a  12.8  cm  diameter,  28  cm 
long,  missile  guidance  section  by  a  300  MHz  plane  wave  at  axial  incidence, 
for  the  case  of  the  interior  dielectric  components,  metal  components,  and 
wires  modeled. 


n  o 


PROGRAM  FDTD  ( 1 NPUT , OUTPUT » T APE b0= I NPUT * T APt 8  = T APE  « > 


C 

C  PUN  TASKS —  INTERFACING  THE  PD-TD  METHOD  wITh  ThE  NET HOf) 

C  OP  MOMENTS 

C  CASE  III-  STEADY  300  MHZ  PLANE  RAVE  IRRADIATION  OP  A 

C  12. «  CM  DIAMETER  MISSILE  GUIDANCE.  SECTION 

C  INTERIOR  METAL  AND  OIELECTPIC  COMPONENTS  MODELED 

C  AXIAL-INCIDENCE  PLANE  WAVE  WITH  COMPONENTS  tZ  AND  HX  MODELED 

C  USING  SCUtLKUNOFF  EQU I  V ALt  NCE  THEOREM  EUR  WILTON/ 

C  GLISSON  M-O-M  APERTURE  CURRENT  DATA 

C  24  X  100  X  48  CELL  CUBIC  SRACt  LATTICE  IS  USED 

C  UNIT  CELL  DIAMETER  =  DX  =  0.33  CM  =  wAVELtNGTh/300 

C  EVEN  SYMMETRY  ABOUT  LATTICE  PLANE  X  =  24,S*DX  IS  ASSUMED 

C  SOFT  LATTICE  TRUNCATIONS  ARE  USED 

C  PROGRAM  IS  OPTIMIZED  FOR  ThE  CDC  STAR-100 

C 

REAL  MUR*MUZ» JXMAG. JXPHA. JXSrC* JYMAGtUYPHA* JYSRC* JZMAG* 

1  JZPH A  »  JZ  SR  C 

dimension  a  <  is65soo  >  ,z  <  issoo  ,aaa  azooi  »huf  <t>oo)»AA(48), 

1  DD(24) ,D£(48) ,DP (48) ,CA(R) ,CB<9) »£PS(R) »SIG(9> 

DIMENSION  JXM AG (  1200)  ,JXPHA(1200>  ,JXSRC(12U0)  , 

1  JYMAG ( 1200) »JYPHA ( 1200) » JYSRC < 1200) , 

2  JZMAG  (1200)  ,J7PHA ( 1/00) ,JZSRC (1200) 

COMMON  A 

Of SCRIPTOB  D1 »U2,D3»D4»D5»D6»NPD,A1,A2»HU,BV,Bw»BX»BY»8Z 
DESCRIPTOR  BQ 

BIT  BU, BbU ( 60 0 ) ,bV  « BBV (  1200) * BW ,HBW ( 1 20 0 ) ,ax ,HHX ( 12  00 ) * 

1  B Y , bb Y ( 24 )  , B7 »  BBZ ( 48 ) 

bIT  BQ,  BBQ (1200) 

ASSIGN  BU,BBU<1»600) 

ASSIGN  BV , BBV (111200) 

ASSIGN  HW , BBW ( 1 M  20  0 ) 

ASSIGN  fax  ,  BBX  ( 1  1  1200) 

ASSIGN  BY  ,  PiBY  (1  J24) 

ASSIGN  BZ, RHZ(l?4a) 

ASSIGN  BG,PBG(1U200> 

ASSIGN  04,DD(1»24) 

ASSIGN  D5,DE(lJ4b) 

ASSIGN  D6 , OF ( 1  I  48 ) 

T 1  =  SECOND ( CP ) 

PRINT  150,  T 1 
150  FORMAT (F20. 5) 

. I.  problem  parameters . 

FREQ  =  3.0E*8 
DX  =  0.01/3.0 

MPR  =  Q 

OATA  EPS  /1.0*  1.0,  1.0,  5.5,  4.B,  5.3,  1.0,  1.0,  8.0/ 

OATA  SIG  /0.0«  3. 7E  *7,  0.025,  0,0024,  0.0008,  0.0,  8.8E*b, 

1  3.7E*7,  0.01/ 

NMAX  =  1800 
C 
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on  non 


C  . II.  nASIC  AND  DENlVtD  CONSTANTS.... 

PI  =  3.1415928b 
MUZ  =  4.0  *  PI  *  1.0E-7 
EPSZ  =  8  .  8b4£  -12 
DT  =  OX  /  6 . 0t  *8 
NHALF  =  0.5  /  EREU  /  OT 
P  =  DT  /  2.0  /  EPSZ 
HA  =  DT**?  /  UA**2  /  MUZ  /  tRSZ 
88  =  DT  /  OX  /  MUZ 
«C  =  l.G/PP 

80  =  2.0  *  PI  *  FRED  *  DT 
I  PUN  =  0 

C A  < 15  9)  =  0. 

Cb ( l ; 9 )  =  o. 

riv»  *  G8VMK0  (  1  ,?s;pw) 

8 X  =  (58VMKZ  (?4.2S?hX) 

UO  ?  1=1. MPR 

EAF  =  R  *  SIG ( I )  /  EPS  ( I ) 

C A ( I )  =  (1.0  -  EAF)  /  (1.0  ♦  tAF) 

?  CP  (  I )  =  PA  /  E8S(I)  /  (1.0  ♦  EAF) 

. III.  LOAD  VECTOR  A . 

.....ZERO  INITIAL  FIELDS . 

Z  ( 1  U5500)  =  0. 

A ( 1 515500)  =  0. 

A ( 1550001 5 15500 )  =  0. 

. TYPt  OF  MEDIUM..... 

WEAD(60.4.E  NO =20 1. ERR®? 01)  (Z ( I > . 1  =  1 .600 ) 

4  FORMAT (7SF1 .0) 

RFAD<60.4,END=201,ERR=201)  (Z(I> , 1=5001 .5600) 

READ (60.4, END  =  ?01 .ERR=201 )  (Z  ( I ) .1=2501*3100) 
DO  5  J=2.14 
JDEL  =  (J-l)  *  15500 

5  A ( JDEL* 1 ) 1 550  0 )  =  Z  ( 1 5 15500 ) 

DO  8  JA=15»71 .7 

READ(60.4.END=201 .ERR=?01 )  ( Z ( I ) , 1= 1 . 60 0 ) 

READ(60.4,END=201.EHR=201)  ( Z  < I ) .1=5001.5000) 

8FAD(h0.A,END=201 .ERR=201 >  (Z ( I ) . 1=2501*3100) 
JOEL  =  (JA-1)  *  15500 
A <JDEL*1 ;  15500)  =  Z ( 1  1 1 5500 ) 

READ(60,4,ENO  =  201 .ERR  =  201 )  (Z ( I ) .  1  =  1 .600 ) 
RFAD(60.4 .END  =  2  0 1 .ERR  =  201 )  ( Z ( I ) *  1  =  5001*560  0) 

00  6  JR  =  1  *  6 

JOEL  =  ( JA ♦ JH- 1 )  *  15500 

6  A  ( JDEL*  1)  15500)  =  ZUUbbOO) 

DO  7  J=93 ♦ 1 0 0 
JDEL  =  (J-l)  *  15500 

7  A  (JOEL*l  » 15500)  =  ZUU5500) 

RE  AD (60«4,END  =  ?01 .ERR  =  ?01 )  ( Z ( 1 ) *  I  =  1  *  60 0 ) 
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n  o 


RE  AO  <60,4  ,ENO=Z01 ,ERR=?Ol  i  <Z  ( 1  )  ,1=5001,5600) 
»EAD(60,4,END=201 ,ERR=201 >  (Z ( 1 ) * 1=2501 *3100) 
00  8  J=75,77 
JDEL  =  (J-l)  *  15500 

8  A  ( JDEL*1  115500)  =  ZUU5500) 

C 

READ(60,4,END=201 ,ERP=201 )  <Z (I) «I=?b 01*310 0) 
J  =  78 

JDEL  =  (J-l)  *  15500 
A(JDEL*1  U5500)  =  Z(1U5500> 

C 

READ (60*4  ,  EnD  =  2  0 1,ERP=?01)  (  Z  (  1  )  *1  =  1*80  0) 

READ(60*4,END=201 ,ENrt=?01 )  (Z ( I ) * 1=5001 ,560  0 ) 
J  =  79 

JDEL  =  (J-l)  *  15500 
A  ( JDEL  *1(  1  550  0  )  =  Z(lUbbOO) 

C 

READ(60»4,EnD=201,ERR=201>  ( 2  ( I ) ♦ I  =  1 ,600 > 
READ(60.4,END=201  , ERR =201  )  <Z  (1)  .  1*5001  ,5600) 
READ(60«4,END=201*tPR=201)  (Z ( 1 ) .1=2501,3100) 
DO  9  J=60  «  84 
JOEL  =  (J-l)  *  15500 

9  A  (JOEL+1*  15500)  =  ZUU5500) 

C 

READ(60»4,END=201,ERR=?01 )  (Z (I) ,1=2501*3100) 
J  =  85 

JDEL  =  (J-l)  *  15500 
A (JDEL*11 15500)  =  Z  < 1 1 15500 ) 

C 

HEAD(60,4,ENO=201 ,ERR=201 )  (Z ( I ) , 1=1 ,600) 

READ  (80,4, END  =  201  ,ERR=201  )  (Z(I)  ,  1  =  50  01  ,5t>0  0  ) 

READ (60»4fEND=201 »£RR=201 )  (Z( 1) ,1=2501,3100) 

00  10  J  =  86 ,87 
JDEL  =  (J-l)  *  15500 

10  A (JOEL* 1 1 1550 0 )  =  Z  ( 1  *  1 S500 ) 

C 

READ (60,4,END  =  201 ,ERR  =  201 )  (Z ( 1 )  ,  1  =  1 ,600) 

READ (60,4, END  =  2 01 ,ERR  =  201 )  ( Z  ( 1 )  ,  I =50 0 1 , 56u 0 ) 

READ(60,4,ENO=201 ,ERR=201 )  (Z ( 1 ) ,1=2501,3100) 
DO  11  J=B8 ,91 
JDEL  =  (J-l)  *  15500 

11  A  (JDEL+1M5500)  =  Z(l?15500> 

REm 0(60,*,, END  =  201, ERR =201 )  (Z  ( I  )  ,  1  =  1 ,600) 

READ (60,4,END  =  201 ,ERR  =  201 )  (Z(I)  ,1=5001,560  0) 
J  =  92 

JDEL  =  (J-l)  *  15500 
A  (JDEL*1 1 15500)  =  ZUU5500) 

.... .APERTURE  EXCITATION . 

RFAD (60,15, ENU=201 «ERR=201 ) (JXRAG(I) ,1=1,625) 
READ(60,16,ENU=?01 ,ERR=201 ) (JXRHA (I ) ,1=1,625) 
RFAD(60*15,ENU=?01*ERR=?01) (JYRAG(I) ,1=1, 62b) 


i 


77 


1 


I 
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RF  AD(60»16,ENU=201.ERR=201>  (JYPHA(I)  , 1=1 ♦ 62b) 
RE AO < 60 ♦lb.FNU=201«ERR=20l)  <JZMAG(I  )  ♦  1  =  1  » 62b) 
READ  (60  ♦  16«END=20  1  »  ERR  =  2G  1 )  (JZPHA(I)  ♦  1  =  1  ♦  62b) 
lb  FORMAT ( 12F6.2«/«13F6.2) 

16  FOPMAT(2bF3.0) 

C 

C  . PARTIAL  SYMMETRY  AdOuT  Z  =  24.0*UX. 

DO  1?  J  =  2 »  10  0 
JOEL  =  (J-l)  *  lbbOO 
KMAX  =  4B 
DO  1?  JA= 1 «  3 

JB  =  JDEL  ♦  (JA-l)*?bOO  ♦  1 
riUF(U600)  =  A  (  JB  16  00) 

PD  =  BUF ( 11600) .GT.6.5 

BUF  (  1  5600)  =  JrtVCTPL  (  3.0»BU»OUF  {  H600)  ) 

IF  ( JA.EGl.3)  KMAX=48 
4(JB*600?2b)  =  BUF (576125) 

DO  12  K=1.24 
KA  =  (K-l)  *  2b 
KB  =  (KMAX-K)  *  25 
1?  A(JB*KB;?5)  =  dUF (K A* 1 J25) 

C 

c  . NON-SYMMETRIC  GROUND  . . 

DO  13  J=43,S7 
JDEL  =  <J“1)  *  15S00 
A ( JDEL>  3 1 99 )  =  7.0 
A ( JDEL ♦ 3224 )  =  7.0 

13  A ( JDEL  *6699 )  =  7.0 
C 

DO  14  K  =  1 8  »  28 
KDFL  =  (K-l  ) *2b  ♦  24 
A  ( JDEL*2bGO*KOtL)  =  7.0 
A  (jDEL*5000«-KUtL)  =  7.0 

14  A  (JDEL*20B00*KUEL)  =  7.0 
C 

C  . APERTURE  EXCITATION . 

DO  17  K  =2 ♦ 24 
KA  =  (K-1)*2B  ♦  1 
KP  =  (49-K ) *?b  ♦  1 

JXMAG(KB;?5)  =  -JXMAG <KA»2b) 

JXPHA(K6l2b)  =  JXPHA ( KA 1 ?b) 

J YMAG (Kb  I  25 )  =  - JYMAG (K A  1 2b ) 

JYPHA(KB;?5)  =  JYPH A ( K  A  I  2b ) 

KA  =  KA  -  2b 

JZMAG ( KB  1 ?b )  =  JZM AG { K A  1 2b ) 

17' JZPHA (KBJ 2b)  =  JZPHA(KA»2b) 

C 

C  . MEDIA  COEFFICIENTS . 

ASSIGN  A1 ,.DYN. 1200 
ASSIGN  A2».DYN.1200 
ASSIGN  NFD..DYN.1200 
A 1  =  0. 

A  2  =  0. 
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DO  21  J=?,100 

JDt  L  =  (J-l)  *  15500 

DO  21  JA=1,3 

JB  =  ( JA- 1 ) *250  0  ♦  1 

JC  a  JP  ♦  1250 

NFD  =  A  ( JOEL  ♦  JH 1 1 20  0 ) 

DO  20  JJ=1»BPH 
BV  =  NFD.fcG.JJ 
A  1  r  QflVCTRL (LA ( JJ)  ,HV»A1 ) 

20  A2  =  Q8VCTRL (Cb < JJ) ,BV» A2) 

A( JOEL* OB  11200)  =  A  1 

21  A ( JOEL ♦ JC i  1 20  0 )  =  A2 
FREE 

T ?  =  SECOND (CH) 

PRINT  150,  T2 

. IV.  TlMt-STEPPINb  LOOP . 

DO  200  N=1,NNAX 
ASSIGN  D1..DYN.1200 
ASSIGN  D2..0YN.1200 
ASSIGN  03  «  « OYN . 1200 

D1  =  FLOAT ( N)  -  JXPHA ( 1  I  1200) 

BO  =  D1.LT.0.0 

D1  =  Q8VCTPL(0.0,BQ;D1) 

02  =  D1  *  PO 
03  s  VSIMU2JJ3) 

JXSRC (1 .1200)  =  JXPAG(l;t200)  *  CB  ( 1  )  *  03 

01  =  FLOAT(N)  -  JYPHA (111  200 ) 

HO  =  D1.LT.0.0 
D 1  =  OflVCTHL<0.0,HQ|L)l) 

02  =  D1  *  HD 
03  =  VSIN(L)?|D3) 

JYSRC(U1200)  =  JYMAG(1J1200)  *  CB  ( 1 )  *  03 

01  =  FLOAT(N)  -  JZPHA (111200) 

BO  =  Dl.LT.0.0 

01  =  O0VCTRL(O.O,RGIO1) 

02  =  01  *  PO 
03  =  VSIN(D?»D3) 

JZSRCUHZOO)  =  JZB  AG  (111200)  *  CP  (  1 )  *  03 
FREE 

MCALL  =  33 


. TRANSVERSE  PLANt  NO.  1 

. ex,  t:  truncations.... 

A ( 750 1  I  1250 )  =  A ( 1 ; 1250) 

A (111250)  =  A (23001 J 1 250 ) 

A ( 1 0201  * 1250)  =  A ( 1 25 1 1 1 250 ) 

A ( 1 25 1 1  1 250 )  =  A ( 25701 f 1250 ) 
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. EY  ITERATION . 

ASSIGN  01 « .OYN.l 174 
ASSIGN  02. .OYN. 1174 
ASSIGN  03*. OYN. 11 74 

A (8826 )  =  O.S  *  < A ( 10078) *A ( 10177) ) 

A ( 8827 5 22 )  =  0.333  *  ( A < 1 00  To  1 22 ) + A ( 1 0077 * 22 ) ♦ A ( 1 0078 1 22 ) ) 

A { 8849 )  =  0.333  *  < A ( 1 0  OSH ) +  2 . 0*A  ( 1 0099 ) ) 

A(1007f>;24)  =  A  (  8hS  1124) 

A  (  1  0026 )  =  0.5  *  (A  UOlOl )  *A  (  10102)  ) 

A  ( 1 0027.22)  =  0.333  *  ( A ( 1 0 1 0 1 5 22 ) ♦ A ( 1 0 1 021 22 ) ♦ A ( 1 0 1 03 5 22 ) ) 

A ( 1 0049 )  =  0.333  *  (A(10l23)*2.0*A(10124n 
A (  1  01  0  1  524)  =  A  (  10001  524) 

L)1  =  CA  (3)  *  A  (8851  5  1174) 

02  =  A(llbSi;il74)-A(llb?651174)4A(l420l51174)-A(142025117*) 
03  =  CR ( 3 )  *  02 
A (885 151174)  =01+03 
FREE 


. TRANSVERSE  PLANES  2  -  100 . 

00  82  JY= 1 « mc ALL 
JOEL  =  <  J Y— 1 )  *  4bS00 

. tX  ITERATION . 

ASSIGN  01  * . OYN. 1173 
ASSIGN  D2..DYN.1173 
ASSIGN  03». OYN. 1173 

00  30  MA=1,3 
M  =  JOEL  ♦  1 5S0  0* ( mA-  1 ) 

. SOFT  LATTICE  TRUNCmT 1 ONS . 

A(23002*M)  =  0.5  *  ( A ( 242&2+M ) ♦ A ( 242S3+M ) ) 

A (23003+R521 )  =  0.333  *  < A ( 24252+M 5 2 1 ) ♦ A ( 24253+M 5 21 ) 

1  ♦ A ( 24254  +  R  5  21 ) ) 

A(?30?4*R)  =  0.5  »  (A(?4?73+w)+A(24274+M) ) 

A (24252+N523)  =  A < 230 27 +M 5 23 ) 

A(24?02+R)  =  0.5  *  (A(24277+N) +A(24278+M) ) 

A (24203*M5?1 )  =  0.333  *  < A  ( 24277+M 5 2 1 ) ♦ A ( 24278+R 5 21 ) 

1  +A (24279+M521 ) ) 

A ( 24224 ♦ f* )  =  O.S  *  (A  (24298  +  R) +A (24299)  ) 

A ( 24277  +  M  5  23)  =  A ( 24 1 77 +M 5 23 > 

. RAIN  EX  . . 

01  =  A ( 15527+M5 1173)  *  A  ( 23 02 7 +  M 5 1 1 73 ) 

02  =  A (29702+R. 1 173) -A ( 14202+R5 1173) +A(2835?+M5 1173) 

1  -A(26377+M51173) 

03  =  A ( 1G777+M5 1 173)  »  0? 

A (23027 5 1 173)  =  01  ♦  03 
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. OX  APERTURE  SOURCE  CONDITION . 

IF(M.EQ. 12*15500)  A ( 2300 1 *M » 1 200 )  =  OXSRC ( 1 » 1200)  ♦ 
1  A  (230  01  ♦M  1 20 0  > 

. EX  ENVELOPE  COMPUTATION . 

04  =  VABS(A(23b01>M»24UD4) 

BY  a  D4.GT.A (24301 *MJ24) 

30  A  ( ?430 1  ♦P I  ?4 )  =  QEVCTRL(D4»t<Y»A<24301>M»24)  ) 

FREE 


.....EY  ITERATION 
ASSIGN  Dlt.0YN.1174 
ASSIGN  D2..0YN.1174 
ASSIGN  03t.0YN.H74 


DO  40  MA=1.3 
M  a  JOEL  ♦  15500*(MA-1) 


A 

A 

I 

A 

A 


<2432b*M) 

<24327*M 

(2434Q*M) 
(25 576*M» 


SOFT  LATTICE  TRUNCATIONS . 

a  0.5  *  (A ( 25576+M ) ♦  A ( 25577 ♦m ) ) 

22)  =  0.333  *  (A (25b76.M;22) (  2:>57 7*M  J  22 ) 

♦  A  (2557K.M22)  ) 

a  0.333  *  (A <2559tt*M) .2.0*A (255*5.M) ) 

24)  =  A (24351 .M I  24) 


A(25526*N)  = 
A ( 25527 ♦M 122 

1 


0.5  *  (A (25b01*M> «-A (25602*M) ) 

=  0.333  »  (A(25b01*MJ22) ♦  A  ( 25b02*M  22 ) 

♦  A  (25b03*M.22)  ) 


A<?554R*M)  a  0.333  *  ( A (2bb23*M) *2. 0*A (25b24.R) ) 
A ( 25601 1 24 )  =  A ( 25501 .M  $  24 ) 


.....FAIN  EY  LOOPS.,... 

D 1  a  A(lfl02b*M;il74)  *  A (24351 1 174) 

02  a  A(27051*MU174)-A(?702b.M»1174)*A(29m*M»ll74) 
1  -A(29702*M»1174) 

03  =  A(1927b*M»1174)  *  02 
A (24351*M1 1 174)  a  01  ♦  03 

. JY  APERTURE  SOURCE  CONDITION . 

IF  (M.EO.b3*15bOO)  A  (2432b*Ml  120  0)  =  JYbRC ( 1 « 1 200 )  ♦ 

1  A  (2432o.m| 1200) 

C 

C  . EY  ENVELOPE  COMPUTATION . 

04  a  VABS < A (?49?6.M J ?4 ) $D4) 

BY  =  04.GT.A (25b2b»M;24) 

A ( 25626  124)  =  Q8VCTRL (D4.bY  >A (25b2b*MJ24)  ) 

C 

05  a  Q0VCMPRS(A(24325*m;12OO) t BX  J  05) 

Ob  a  VABS (05 1 Ub ) 

BZ  =  D6.GT.A (25651.MI4B) 

40  A ( 25651 4M I 4fi )  =  QbVCTRL ( D6 t WZ f A ( 2565 1 +M 1 4d ) ) 

FREE 
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. EZ  ITERATION . 

ASSIGN  01  *  .OYN. 1 199 
ASSIGN  D2..DYN.1199 
ASSIGN  03, ,0YN. 1 199 

00  SO  MA=1,3 
M  =  JOEL  ♦  15500* (MA-1 ) 

. MAIN  EZ  LOOPS..... 

D 1  =  A(20S01*MU199)  *  A ( 2570 1 *M  1 1 199 ) 

D?  =  A(?ft35?*M*1199)-A(283Sl*NJ1199)*A(ll 526+M I 1199) 
1  -A(27026*M»1199) 

03  =  A (21751*m? 1199)  *  02 
A ( 257  0 1 ♦N 1 1 199 )  =  01  ♦  03 

. JZ  APERTURE  SOURCE  CONDITION . 

IF (M.EQ. 12*15500)  A<25701*M1200>  a  JZSRC(1I1200)  ♦ 

1  A (2570l*M| 1200) 


. EZ  ENVELOPE  COMPUTATION . 

A ( 26926*M  1 24 )  *  0.5  *  <A(26276*MI24) ♦ A ( 2630 1 *M| 24 ) ) 
04  =  VAbS (A (26926*MI24) *04) 

BY  a  04 .GT . A ( 2695 1 ♦M l 24 ) 

A ( 26951 ♦  M  5  24 )  =  QbVCTRL (D4 * «Y t A ( 26951 *M 124 ) ) 

05  =  08VCMPPS (A (25700*M»1200) *«x;D5) 

06  a  VAbS  (05  }  1)6  ) 

HZ  a  06.GT. A (26976>M»48) 

50  A (26976*M{48)  =  G8VCTRL (D6f oZ *  A ( 26976*M ;4« ) ) 

. HX  ITERATION..... 

00  60  M A= 1 . 3 

M  a  JOEL  ♦  15500MMA-1) 

. MAIN  HX  LOOPS . 

D 1  a  A ( 1 1526*M? 1199) 

02  a  A(8851*m;1199)-A(882G*m|1199>»A(10201*M»1199) 

1  -A(25701*MU199) 

A  <  1  15?6-»MU199)  a  oi  ♦  02 

. HX  ENVELOPE  COMPUTATION..... 

a ( 1 2751 ♦Ml 24 )  =  0.5  *  (A(12101+MI?4)*A<12126*M|24) ) 

04  a  VABS (A ( 12751*MI24) »D4) 

HY  =  D4.GT.A ( 12776*9124) 

A (12776*m;24)  =  Q8VCTRL ( D4 ♦ bY *  A ( 1 2776*M? 24 ) ) 

C 

05  =  08VCMPRS<A <115?5*MI1200> *bX;D5) 

06  a  VABS (05.06) 

HZ  a  D6.GT. A ( 12801*M|48) 

60  A ( 1 280 1  *M 1 48 )  =  Q8VCTPL<D6*o2)A(12801*Ml4B) ) 

FPEE 
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C  . HY  ITERATION.,... 

ASSIGN  D1..DYN.1198 
ASSIGN  0?,.t)YN.H% 

00  70  MA= 1*3 
M  =  JDFl  ♦  lSbOG*(MA-l) 

. SOE T  LATTICE  TRUNCATION . 

A  A  (  1 )  =  O.S  «  (A  (  ?9f>0  1  ♦  M  )  ♦A  ( 296  ()2*M )  ) 

AA { 2  146 1  =  0.333  *  ( A ( 2960 1 *M 1 46 ) ♦ A (296 02*M 1 46 ) 

1  *A  (?Vt>03*M;4b)  ) 

AA<48)  =  O.b  *  (A  (29647*8  )*M(*'9648.M)  ) 

A  A  A  (  1  >1200)  =  U8VXRN0I  AA  (  1  ;4«)  ,bw  ;  AAA  (  1  I  120(>)  ) 

A (29601 ♦Ml  48)  =  QMVCMRRS < A ( 2«3b?*Mt  1200) *** ? A <29601 *Ml4R)  ) 

.....MAIN  RY  LOOPS..... 

01  =  A<283b2*M*1196) 

02  =  A  <2b70?*Ml  1 198) -A  <2S701*m;  1 19ft  )  *a  (23<j08*m;  1 196) 

1  -A  (23l>27*M{  119b  ) 

A  (?83S2*M1 1 196  )  =  01  ♦  02 

A  <283M*MI  1?00)  =  Q8VCTRL  (AAA  (1  *1200)  .  b«i  5  A  (  3S 1  ?  1  ?  00  )  ) 

A  (26351*8  11200)  =  Q8VCTRL  <  0 . 0  *  6  X  *  A  (  283b  1  *t*  *  1 200  )  ) 

.....RY  ENVELOPE  COMPUTATION . 

A(?9576*m»2A)  =  O.S  *  <A(2R92«->*MJ?4)*A(2rt>yi>l*m?4)) 

04  =  VA6S (A <29b76*M*24) *Q4) 

HY  s  04  »GT .A ( 2965 l ♦M 1 24 ) 

70  A  ( 29651  ♦H(?4)  =  08VCTRL(f>4»nY*A(296bl*M*24)  ) 

EPEE 

. HZ  ITERATION . 

ASSIGN  D1..0YN.1173 
ASSIGN  02 » . DYN. 1173 

00  80  MA=1 ,3 
M  =  JOEL  ♦  1SB00* (MA-l ) 

.....SOFT  LATTICE  TRUNCATION . 

A A ( 2 )  =  O.b  *  <A<lb427«-M)*A(lb42R*M)  ) 

AA(3*45)  =  0,333  *  < A ( lb427*M*4b) .A ( 1542B*M*4b> 

1  *A ( lb429*M*4b)  ) 

AA ( 48 )  =  O.S  *  (A(lS472*M)4A<lb473*M) , 

AAA(1?1200)  =  08 V  XPNO { A  A  < 1 ?  48 ) «  hw  *  A  A A (1 *  1 20  0 )  ) 

A  <15426*MJ4fi)  =  Q8VCMPRS(A(14177*MU200)  » b* *  a ( 1 b42b*M * 48 )  ) 

.....MAIN  HZ  LOOPS..... 

01  =  A(14202»m; 1173) 

02  =  A <23027*M» 1 173) -A <7b27*M* 1 173) *A (BHbl*"* 1 173) 

1  -A  <68b2*M  1173) 

A  (14202*H» 1173)  *  01  ♦  DP 
C 

A  (1417fc*HU200)  *  08VCTRL  (AAA  ( 1 J  1200)  .BR  I  A  (  1417h*m»  1200)  ) 
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A  <1417f»*M*l?00)  =  GHVCTPl  <0.0*hX  J  A  ( 1417b«-M;  1200)  ) 

. HZ  ENVELOPE  COMPUTATION . 

04  =  VAbS (A ( 14/7b>M|?4) JD4) 
rY  =  04.0T.A  (  1  b4  7b  ♦*t  ;  ?*  ) 

Hn  A < l547b*M5?4)  =  (JHVCTWL  <04,nY  *  A  (  lb47b^M»24)  > 

f-  PFfc 

HP  CONTINUE 

IF (rCALL.Lt .32) GO  TO  94 

. TRANSVERSE  PLANt  NO.  101 . 

m  =  99  *  15500 


. EX»  tz  TRUNCATIONS . 

A  <?3001*P*.  1250)  =  A  ( 1  550 1  *M  *  1  2b  0 ) 
A  U5501*m;1?S0)  =  A (16751*Mil2S0) 
A ( 16751*F U250)  =  A <7501*M; 12S0) 


A  (  257  01  ♦  F"  *  1  ?5  0  ) 
A ( 1 P001 ♦H  » 1250) 
A (19?S1*H? 1P5U) 


A (1H001*H« 1250) 
A  <  19251  ♦<*)  1250) 
A ( 10201*Mtl2b0) 


. HX  ITERATION . 

ASSIGN  U1 % .DYN. 1 199 
ASSIGN  02  *  •  OYN  *11  99 
01  =  A< llb2b*H» 1199) 

0?  s  A  (  8b5 1  ♦M  •  I  199)  -  a  (Hb?b*b5  1 199)  *  A  {  10201*Ml  11  99 ) 

1  -A(25701  +  h»hq9) 

A(1 15?G*m; 1199)  =  01  ♦  0? 

A ( 12751 *MI24)  =  0.5  *  (A ( 1210 1 ♦  Mi  24) ♦  A ( 1212b*MI ?4 ) ) 
04  =  VAtS(A(12751*w»?4) ;D4) 
riY  =  04. GT. A ( 1277G*m;24) 

A ( 1 2776 ♦  M • ?4 )  =  Q»VCT«l <D4.»Y ) A ( l?776*Mi24) ) 

i>5  =  OflVCMPRS  <A  (1 15?5*MJ  1200)  tbXIOS) 

Ob  =  VAhS(055Ub) 

07  =  Ob.GT.A ( 12801+MS45) 

A  (  12801O54B)  =  QrVCTRL(nfo*oZ*A(l?bOl4M*4on 
FREE 


HZ  ITERATION..... 

ASSIGN  01  ,  . l>Y<v.  1 173 
ASSIGN  U?»  .OYN.  1 173 

A  A ( ? )  =  0.5  *  ( A ( 15427*M) +A { lb42H*M) ) 

AA(3*45)  =  0.330  *  (A < 15427**;4S) ♦ A ( 1 5428 ♦M » 4 5 ) 

1  ♦  A ( 15429+HJ45) ) 

A  A ( 48 )  =  0.5  *  ( A ( 1547?*P> ♦ A < 15473«M)  ) 

AAA(i;i?00)  =  UMVXPNO  <  AA  (  H4h>  *PwiAAA  ( 1  J  1200)  ) 

A ( 1 5425 ♦P 148)  =  ObVCMPPS { A ( l4l77*M? 1200) » A ( lb42b*M|4P) ) 

C 

01  =  A(l4?G?**f 1173) 


04 


no  r>  o 


02  =  A <23027*M11 173>-A (7527*m»1 173)  *A  <8«51*M  1 173) 

1  -A (B852*m; 1 173) 

A  ( 14?0?*M»1 173)  =  D1  ♦  02 

A  ( 1  4 1 76  11200)  =  GHVCTRL  ( AA  A  ( 1  1 1 20  0  )  ,HW  I  A  (  1 4  1 7b*M  U200  >  ) 

A ( 141 76*M  1 1 200 )  =  08VCTRL < 0. 0 »HX *  A < 14l7b*M » 1 20  0 ) ) 

04  =  VA8S(A(14776*MI24) ID4) 

HY  =  D4.GT.A( 15476«.M|24) 

A  ( 1 5476+M  I  24 )  =  QWVCTRL  (D4,bY|A (1647b4M»24) ) 

FREE 

94  T 3  =  SECOND ( CP ) 

PRINT  150 ♦  T3 

. FIELD  ENVELOPE  PRINTOUT  ROUTINE . 

00  100  L  =  NHALf’  »NMAX»NHALF 
IF (N.EO.L  >  GO  TO  101 
100  CONTINUE 

IF(N.EQ.NMAX)GO  TO  101 
60  TO  199 


AT  HORIZONTAL  OBSERVATION  PLANE . 

101  IF (N.EO.NMAX) IHUN=1 
I PUN  =  0 
PRINT  102,  N 

102  FORMAT  (1H1*S2X*27HEX  ENVELOPE  E OR  TIME  ST eP  =  ,I5, 

1  //, 62X, 15HPLANE  Z  =  24*DX , // ,2X , lHU, / ) 

CALL  ENV(fl«01»20»RC»IPUN) 

C 

PRINT  103,  N 

103  FORMAT (1H1,52X»27HEY  ENVELOPE  E OR  TIME  STEP  =,I5, 

1  //,62X,15HPLANE  Z  =  24*DX  ,//,2X,  lrtj,/) 

CALL  ENV ( 10126, 20, RC* I PUN  > 

C 

PRINT  104,  N 

104  FORMAT ( 1H1 ,52X,27hEZ  ENVELOPE  E OR  TIME  STEP  =,I5, 

1  //,62X, 15HPLANE  Z  =  24*0X,//»2X, 1HJ,/) 

CALL  ENV(11451,20,RC*IPUN) 

C 

PRINT  105,  N 

105  FORMAT ( 1H1 ,5?X,27HHX  ENVELOPE  E OR  TIME  STEP  =,I5, 

1  //,62X» 15HPLANE  Z  =  24*UX , // » 2X , 1HJ , / ) 

CALL  ENV(1277b, 20, 376,7, IRON) 

C 

PRINT  106,  N 

106  FORMAT ( 1H 1 , S2X , 27HHY  ENVELOPE  E OR  TIME  STEP  =,I5, 

1  //,62X, 15HPLANE  Z  =  24*DX,//,2X*1MJ,/) 

CALL  ENV ( 14151 t20, 376. 7, IRON) 

C 

PRINT  107,  N 

107  F0RMAT(1H1,52X,27HHZ  ENVELOPE  FOR  TIME  STEP  =,IS, 

1  //,62X,15HPLANE  Z  =  24*0X,//,2X* lHj,/) 

CALL  ENV(1547b, 20,376. 7, IPUN) 

C 
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C  . AT  VtRTICAL  SYMMETRY  PLANE . 

PRINT  1 0 A  »  N 

108  F0PMAT(1H1*52X»?7HEZ  tNVELOPt  F  OR  TIME  STtP  »*I5 

1  //*6?X» 17 HP LANE  X  =  24.S*DX*//»2X»lRJ*/> 

CALL  ENV ( 1 147&*40*RC« IPUN) 

C 

PRINT  109*  N 

109  FORMAT ( 1 H  1  *  52X  ,  27FihX  ENVELOPt  FOR  TIME  STEP  =  ,I5 

1  //»#»2X»  17HPLANE  X  =  24.®>*DX*//*2X*lFtJ,/) 

CALL  ENV < 12801 *40*376.7, IPUN) 

CUR  1  =  A<49*15500*10178)  *  PC  *  73.33 
CUP2  =  A(7H*lSS00*-101bl )  *  *C  *  73.33 
PRINT  150*  CUPl 
PRINT  150*  CU«2 
C 

PRINT  110*  N 

110  FORMAT  < 1H1*5?X*27HEY  ENVELOPt  FOP  TIMt  STtP  *,I5 

1  //*62X, 17HPLANE  X  =  24.5*DX *//»2X* 1HJ,/) 

CALL  ENV(101S1*40.RC*IPUN) 

C 

19Q  CONTINUE 

200  CONTINUE 

201  T4  =  SECOND (CM) 

PRINT  150*  T  4 
STOP 

END 

SUBROUTINE  tN V < LOCA ,NUM * SCALE ♦ IPUNCH ) 

DIMENSION  A(lbbS5001  *  IP  <40001  »NNt50> 

COMMON  A 
C 

00  1  1=1,50 

1  NN(I)  =  I 
C 

00  2  LY=2*100 

LOC  =  LOCA  ♦  (LY-1 ) *15500  ♦  4 
LOCI  =  1  ♦  <LY-1)*40 
IP(LOCIINUM)  =  SCALE  *  A(LOCINUM) 

2  A(LOC*NUM)  *  0. 

C 

I F ( IPUNCh .EO. 0 ) GO  TO  5 
DO  3  LYS2» 1 00 
LOCI  =  1  ♦  <LY-1)*40 
LOCI I  =  LOCI  -  1  ♦  NUM 

3  MR I TE (8*4)  (IP(LL) ,LL=L0CI*L0C1I) 

4  FORMAT (1016) 

C 

5  LXM  =  NUM/20 
00  8  LX=1,LXM 

LX A  =  1  ♦  <LX-1)*20 
LA  =  LX A  ♦  4 
LZ  =  LA  ♦  19 
00  6  LY=?*100 
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LYY  =  10?  -  LY 

LOCI  =  LX A  ♦  (LYY-1 ) *40 

LOCI  I  =  LOCI  ♦  19 

6  PRINT  7 *  LYY.  <  IP < LL ) . LL  =  LOC I . LOC 1 1 ) 

7  FORMAT  <1X.I3.^X»20I6) 

8  PRINT  9.  (NN(LL) .LL=LA.LZ) 

9  FORMAT (//.«X. 2016. //////) 

RETURN 

END 


3.0  DATA  CARD  FORMAT 


3.1  Type  of  Medium  at  Each  Lattice  Cell 

For  most  realistic  interaction  geometries,  data  cards  are  needed  to 
specify  the  type  of  medium  at  each  location  of  an  electric  field  component. 
Up  to  9  distinct  media  can  be  efficiently  specified  within  a  lattice  using 
the  following  format. 

Using  the  Fortran  statement  4  FORMAT  (75F1.0)  a  medium-type  integer 
1,  2,...,  9  can  be  assigned  to  the  600  locations  of  an  electric  field  com¬ 
ponent  in  one  plane  j  =  constant  with  only  8  data  cards.  The  600  locations 


are  ordered  as  shown  in  Figure  1. 

Data  Card 

Within  the  75F1.0  format,  1 

Assigns  Type  Integer  t 
Consecutive  Locations 

1 

1  -  75. 

2 

76-150 

3 

151  -  225 

4 

226  -  300 

5 

301  -  375 

6 

376  -  450 

7 

451  -  525 

8 

526  -  600 

In  all  data  cards,  column  25,  column  50,  and  columns  75  -  80  are  left  blank. 
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3.2  Equivalent  Electric  Current  Aperture  Excitation 

For  the  case  of  the  hybrid  method  of  moments/FD-TD  analysis,  data  cards 

are  needed  to  specify  the  Schelkunoff  equivalent  electric  current  excitation 

at  each  point  in  the  aperture  plane  marked  by  the  location  of  an  electric 

field  component.  The  electric  current  excitation  is  given  in  amperes/meter 

and  is  normalized  (for  purposes  of  achieving  a  wide-range  print  out)  to  an 

4 

incident  wave  electric  field  of  10  volts/meter.  In  addition,  the  electric 
current  excitation  is  characterized  at  each  point  by  an  integer  number  repre¬ 
senting  the  number  of  time  steps  of  delay  of  the  current  relative  to  the 
earliest  appearance  of  current  at  any  point  within  the  aperture  plane. 

Using  the  Fortran  statement  15  FORMAT  (12F6.2 ,/ ,13F6.2) ,  a  current- 
excitation  magnitude  can  be  assigned  to  the  600  locations  of  an  electric 
field  component  in  one  plane  j  =  constant  with  only  48  data  cards.  The  600 
locations  are  ordered  as  shown  in  Figure  1.  With  this  format,  we  have: 

Assigns  Current  Magnitude 
Data  Cards  to  Consecutive  Locations 


1 

and 

2 

1 

-  25 

3 

and 

4 

26 

-  50 

5 

and 

6 

51 

-  75 

47 

and 

48 

576 

-  600 

In  all  data  cards,  the  quantity  0.0  is  assigned  to  electric  field  component 
locations  either  at  the  boundary  of  the  aperture  or  outside  of  the  aperture. 
In  all  even-numbered  data  cards,  the  quantity  0.0  is  assigned  to  the  thir¬ 
teenth  entry  (columns  73-78). 

Using  the  Fortran  statement  16  FORMAT  (25F3.0),  a  current-excitation 
delay  can  be  assigned  to  the  600  locations  of  an  electric  field  component  in 
one  plane  j  =  constant  with  only  24  data  cards.  The  600  locations  are  order¬ 
ed  as  shown  in  Figure  1.  With  this  format,  we  have: 
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Assigns  Current  Delay 

Data  Card  to  Consecutive  Locations 


1 

1 

-  25 

2 

26 

-  50 

3 

51 

-  75 

24 

576 

-  600 

In  all  data  cards,  the  quantity  0  is  assigned  to  electric  field  component 
locations  either  at  the  boundary  of  the  aperture  or  outside  of  the  aperture. 
Further,  in  all  cards,  the  quantity  0  is  assigned  to  the  twenty-fifth  entry 
(columns  73-75). 
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